用matlab实现连续时间周期信号(方波信号，三角波信号）的频谱分析

Posted 2023-03-06

调整信号的周期和宽度，观察信号谱线的特性。 方波信号：周期T=1s，T1=0.5s，A=±1. 三角波信号：周期T=1s，A=1. 求代码！

参考技术A %产生峰值为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;