Question

用matlab绘出周期锯齿波信号的时域波形及频谱图，T=2pai

Answer 1

%产峰值1三角波析其0~63谐波幅值谱相位谱 clf; Fs =128; %采频率 T = 1/Fs; % 采周期 N = 128; % 采点数 t = (0:N-1)*T; % 间单位：S x=zeros(N); for n=0:N-1 b=fix((n)/(N/4)); m=n+1; A=1/(N/4); if b==0 x(m)=A*n; elseif b==1||b==2 x(m)=A*(N/2-n); elseif b==3 x(m)=A*(n-N); end; end; n=0:N-1; subplot(3,1,1) plot(t,x); xlabel('间/S'); ylabel('振幅');title('域波形');grid on; y=fft(x,N); %信号进行快速Fourier变换 mag=abs(y)*2/N; %求取Fourier变换振幅;*2/N转变真实幅值 f=n*Fs/N; subplot(3,1,2) plot(f(1:N/2),mag(1:N/2)); %绘Nyquist频率前随频率变化振幅 xlabel('频率/Hz'); ylabel('振幅');title('幅值谱');grid on; p=mod(angle(y)*180/pi,360); subplot(3,1,3) plot(f(1:N/2),p(1:N/2)); %绘Nyquist频率前随频率变化相位 xlabel('频率/Hz'); ylabel('振幅');title('相位谱');grid on;