作业帮 > 综合 > 作业

求matlab周期三角波信号频谱分析的代码,能画出三角波信号、幅度谱和相位谱.

来源:学生作业帮 编辑:大师作文网作业帮 分类:综合作业 时间:2024/09/20 11:54:22
求matlab周期三角波信号频谱分析的代码,能画出三角波信号、幅度谱和相位谱.
求matlab周期三角波信号频谱分析的代码,能画出三角波信号、幅度谱和相位谱.
%产生峰值为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;