MATLAB實(shí)驗(yàn)二-周期信號(hào)的頻譜

實(shí)驗(yàn)二周期信號(hào)的頻譜clearalln_max=[3214181];%計(jì)算N次N=length(n_max);t=-1.1:.002:1.1;omega_0=2*pi;%基波頻率fork=1:Nn=[];n=[1:n_max(k)];b_n=2./(pi*n).*(-1).^(1+n);%計(jì)算傅里葉級(jí)數(shù)anx=b_n*sin(omega_0*n'*t);%計(jì)算前幾項(xiàng)的和subplot(N,1,k),plot(t,x,'linewidth',2);%在N幅圖中的第k子圖畫(huà)波形axis([-1.11.1-1.51.5]);line([-1.11.1],[00],'color','r');%畫(huà)直線，表示橫軸，線為紅色line([00],[-1.51.5],'color','r');%畫(huà)直線，表示縱軸，線為紅色bt=strcat('最高諧波次數(shù)=',num2str(n_max(k)));%字符串連接title(bt);%在N幅圖中的第k子圖上寫(xiě)標(biāo)題endclearalln_max=[3214181];A=1;N=length(n_max);%計(jì)算N次t=-1.1:.002:1.1;omega_0=2*pi;%基波頻率figure;fork=1:Nn=[];n=[1:n_max(k)];%n=1,3,5等b=(pi*n).^2;a=(sin(pi*n/2)).^2;a_n=4./b.*a;%計(jì)算傅里葉級(jí)數(shù)anx=A*a_n*cos(omega_0*n'*t);%計(jì)算前幾項(xiàng)的和subplot(N,1,k),plot(t,x,'linewidth',2);%在N幅圖中的第k子圖畫(huà)波形axis([-1.11.1-1.51.5]);line([-1.11.1],[00],'color','r');%畫(huà)直線，表示橫軸，線為紅色line([00],[-1.51.5],'color','r');%畫(huà)直線，表示縱軸，線為紅色bt=strcat('最高諧波次數(shù)=',num2str(n_max(k)))；%字符串連接title(bt);%在N幅圖中的第k子圖上寫(xiě)標(biāo)題endclearalltau=0.5;t=[2510];tau_T=tau./t;%占空比omega_0=2*pi./tn0=0;n1=50;N=length(t);fork=1:Nn=[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([n0n1Fn_min-0.01Fn_max+0.01]);line([n0n1],[00],'color','r');%畫(huà)直線，表示橫軸，線為紅色title(strcat('幅度頻譜：頻譜密度=0.5,周期=',num2str(t(k))));%在圖上寫(xiě)標(biāo)題endclearalltau=[0.20.51];t=5;tau_T=tau./t;%占空比omega_0=2*pi/tn0=0;n1=50;N=length(tau);n=[n0:omega_0:n1];fork=1:NF_n=tau_T(k)*Sa(tau(k)*n/2);Fn_max=max(F_n);Fn_min=min(F_n);subplot(N,1,k),stem(n,F_n,'.');%畫(huà)幅度頻譜axis([n0n1Fn_min-0.01Fn_max+0.01]);line([n0n1