[求助] 编码 下面函数 画频谱
 

菜鸟遇到了，下面的函数，声音信号的函数模型
N
s（t）=∑cos(t)*cos(i*0.01*pi*t)
i=1

其中：N=6,
t=0：0.01：15；
怎么画出 函数的频谱？
谢谢大侠的帮助

帮顶，同问。。。

先计算时域信号，再通过fft函数获得幅频序列即可

[QUOTE=watcher;9017]先计算时域信号，再通过fft函数获得幅频序列即可...[/QUOTE]
N=6;
n=0;
for t=0:0.01:15;
for k=1:N;
a=cos(t);
y=cos(t*pi*0.1*k);
m=a*y;
end
n=m+n;
hold on;
plot(t,n);
end
上面是我写的时域信号，可是时域信号不是一个函数，而是个点---冲击，所以傅立叶变换后就是一常数....
这是老师给的毕业设计里面语音信号的函数模型。
恳请高手指点

[QUOTE=young266;9034]N=6;
n=0;
for t=0:0.01:15;
for k=1:N;
a=cos(t);
y=cos(t*pi*0.1*k);
m=a*y;
end
n=m+n;
hold on;
pl...[/QUOTE]

外循环没必要，如下程序可生成时域信号序列
N=6;
n=0;
[COLOR="Red"]t=0:0.01:15;[/COLOR]
for k=1:N;
a=cos(t);
y=cos(t*pi*0.1*k);
[COLOR="Red"]m=a.*y;[/COLOR]
[COLOR="red"]n=m+n;
end[/COLOR]
plot(t,n);

太感谢了 呵呵