%产生峰值为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));
Y=fft(y,512);
F =10*f*[0:256]/512;
fp=2*sqrt(Y.*conj(Y));％幅度谱
xp=angle(Y); ％相位谱
gl=abs(Y).^2; %功率谱
magif=ifft2(abs(f2));%幅度重构
pha=angle(f2);%取相位
phaif=ifft2(exp(j*pha));%相位重构
数据是x(i)，共N个点，采样频率是fsample
扩展资料 ：
信号源有很多种，包括正弦波信号源，函数发生器、脉冲发生器、扫描发生器、任意波形发生器、合成信号源等。一般来讲任意波形发生器，是一种特殊的信号源，综合具有其它信号源波形生成能力，因而适合各种仿真实验的需要。
传统都认为信号源主要给被测电路提供所需要的已知信号（各种波形），然后用其它仪表测量感兴趣的参数。可见信号源在电子实验和测试处理中，并不测量任何参数而是根据使用者的要求，仿真各种测试信号，提供给被测电路，以达到测试的需要。
参考资料来源：百度百科-三角波信号
%产生峰值为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;
