武神裝攻 发表于 2010-11-25 08:10:03

请板上板友帮实现集总for回圈的程序实现

不好意思 恳请帮助 谢谢
如下图所示公式
http://addsnaps.com/img/original/26142
小弟假设一下值是
%%%%%%%%%%%
clear all
clc
t = 0:0.005:0.1;
w = 2*pi*60
Vm = ;   %是个a*1大的矩阵
θn = ;
for ot = 1:length(t)
   for n = 1:length(θn)
          m = n ;
               ?
               ?
               ?
   end
end
也就是第一个Vt(1;1)的值
=10*sin(1*w*0+0) + 5*sin(2*w*0+30) + 8*sin(3*w*0+60) + ............+ 2*sin(7*w*0+30)
第二个Vt(1;2)的值
=10*sin(1*w*0.005+0) + 5*sin(2*w*0.005+30) + 8*sin(3*w*0.005+60) + ............+ 2*sin(7*w*0.005+30)
第三个Vt(1;3)的值
=10*sin(1*w*0.01+0) + 5*sin(2*w*0.01+30) + 8*sin(3*w*0.01+60) + ............+ 2*sin(7*w*0.01+30)
第四个Vt(1;4)的值
=10*sin(1*w*0.015+0) + 5*sin(2*w*0.015+30) + 8*sin(3*w*0.015+60) + ............+ 2*sin(7*w*0.015+30)
第五个Vt(1;5)的值
=10*sin(1*w*0.02+0) + 5*sin(2*w*0.02+30) + 8*sin(3*w*0.02+60) + ............+ 2*sin(7*w*0.02+30)
...............依此类推
主要用途是FFT反转
将频谱转回波型
小弟看到此公式是这样想 还是说有错请指教
在此先谢谢帮实现 拜托 谢谢

qibbxxt 发表于 2010-11-25 08:53:21

clear all
clc
t = 0:0.005:0.1;
w = 2*pi*60;
Vm = ;   %是个a*1大的矩阵
theta = ;
n=1:length(Vm);
f=@(t)sum(Vm.*sin(n'.*w.*t+theta));
g=arrayfun(f,t);
页: [1]
查看完整版本: 请板上板友帮实现集总for回圈的程序实现