木的人回答

自己做吧:(

function [t,xt]=MySin(A0,f0,P0,N,fs)

%   MySin产生一个指定幅值和频率的正弦信号.

%   MySin(A0, f0)返回一个时间t数组和一个幅值xt数组.

%   A0：幅值, f0：频率, P0：相位, dt：时间间隔, N：数据点数,fs:采样率

%   例如要生成幅值为5，频率为10Hz，相位为0的正弦信号，如此输入便可：[t,xt]=MySin(5,10,0,2048,1024)

%   当前版本：1.0

%   作  者：richardpan,Matlab

%   完成日期：2009年03月04日

if nargin<5

error('输入变量太少!');

end

if nargin >5

error('输入变量太多!');

end

%A0=5;                                                                                              %幅值为5

%f0=10;                                                                                             %频率为10HZ

%fs=1024;                                                                                            %采样率

t=0:1/fs:1;                                                                                         %采样点个数为length(t)

xt=A0*sin(2*pi*f0*t+P0);                                                                            %根据公式x(t)=A*sin(2*pi*f*t+P)生成指定频率和幅值的正弦信号

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

A0=5;

f0=20;

N=2048;

fs=1024;

P0=pi/2;                                                                   %因为三角波的傅里叶展开式是余弦函数，欲用MySin叠加生成，则需将初相位设为90°

b=f0;                                                                      %将基频赋给b

a=1:2:1024;

PanSin=[];                                                                 %正弦信号叠加数组

for n=1:30                                                                 %n为叠加个数，设最大30

f0=a(n)*b;                                                              %各个谐波频率

[t,xt]=MySin(A0,f0,P0,N,fs);                                            %调用[t,xt]=MySin(A0,f0,P0)函数

PanSin(1:length(t),n)=1/(a(n)^2)*xt                                      %进行叠加

end                                                                        %30次叠加完成后，PanSin数组中是30列长度为length(t)的数据

b=PanSin';                                                                 %求转置

y2=A0/2+sum(b(1:2,:));                                                     %取前2行叠加

y5=A0/2+sum(b(1:5,:));                                                     %取前5行叠加

y10=A0/2+sum(b(1:10,:));                                                   %取前10行叠加

y20=A0/2+sum(b(1:20,:));                                                   %取前20行叠加

y30=A0/2+sum(PanSin');                                                     %取前30行叠加

plot(t(1:103),y2(1:103),'-r');                                             %绘图，只显示2个周期

hold on

plot(t(1:103),y5(1:103),'-g');

plot(t(1:103),y10(1:103),'-y');

plot(t(1:103),y20(1:103),'-m');

plot(t(1:103),y30(1:103),'-b');

title(['用正弦信号构造三角波信号'],'fontsize',14);

xlabel('t(ms)','fontsize',14);                                             %设置x轴标题为t(ms)

ylabel('x(t)','fontsize',14);                                              %设置y轴标题为x(t)

axis tight;

legend('2个正弦信号','5个正弦信号','10个正弦信号','20个正弦信号','30个正弦信号');