MATLAB实验二 周期信号的频谱

实验二 周期信号的频谱

clear all

n_max=[3 21 41 81]; %计算N次 N=length(n_max); t=-1.1:.002:1.1;

omega_0=2*pi; %基波频率

for k=1:N n=[];

n=[1:n_max(k)];

b_n=2./(pi*n).*(-1).^(1+n); %计算傅里叶级数an x=b_n*sin(omega_0*n'*t); %计算前几项的和

subplot(N, 1,k),plot(t,x,'linewidth',2); %在N幅图中的第k子图画波形 axis([-1.1 1.1 -1.5 1.5]);

line([-1.1 1.1],[0 0],'color','r'); %画直线，表示横轴，线为红色 line([0 0],[-1.5 1.5],'color','r'); %画直线，表示纵轴，线为红色 bt=strcat('最高谐波次数=',num2str(n_max(k))); %字符串连接

title(bt); %在N幅图中的第k子图上写标题 end

clear all

n_max=[3 21 41 81]; A=1;

N=length(n_max); %计算N次 t=-1.1:.002:1.1;

omega_0=2*pi; %基波频率 figure; for k=1:N n=[];

n=[1:n_max(k)]; %n=1,3,5等

b=(pi*n).^2;

a=(sin(pi*n/2)).^2;

a_n=4./b.*a; %计算傅里叶级数an x=A*a_n*cos(omega_0*n'*t); %计算前几项的和

subplot(N,1,k),plot(t,x,'linewidth',2); %在N幅图中的第k子图画波形 axis([-1.1 1.1 -1.5 1.5]);

line([-1.1 1.1],[0 0],'color','r'); %画直线，表示横轴，线为红色 line([0 0],[-1.5 1.5],'color','r'); %画直线，表示纵轴，线为红色 bt=strcat('最高谐波次数=',num2str(n_max(k)))； %字符串连接

title(bt); %在N幅图中的第k子图上写标题 end

clear all

tau=0.5;t=[2 5 10];

tau_T=tau./t; %占空比

omega_0 =2*pi./t n0=0;n1=50; N=length(t); for k=1:N

n=[n0:omega_0(k):n1]; F_n=tau_T(k)*Sa(tau*n/2) Fn_max=max(F_n); Fn_min=min(F_n);

subplot(N,1,k),stem(n,F_n,'.');

axis([n0 n1 Fn_min-0.01 Fn_max+0.01]);

line([n0 n1],[0 0],'color','r'); %画直线，表示横轴，线为红色 title(strcat('幅度频谱：频谱密度=0.5,周期=',num2str(t(k)))); %在图上写标题 end