matlab作业:信号y = 5*sin(pi*20*t)+3*cos(2*pi*50*t)
clc;clear;
fs=1000;%采样频率
N=2000;%采样数据点
n=0:N-1;t=n/fs;%时间序列
y_source=5*sin(pi*20*t)+3*cos(2*pi*50*t);%信号
y_whitenoise=rand(1,N);%均匀分布白噪声
% y_whitenoise=randn(1,N);%正态分布白噪声
y_mix=y_source+y_whitenoise;%混合信号
y_fft=fft(y_mix,N)*2/N;
mag=abs(y_fft);
f=n(1,1:N/2)*fs/N;
subplot(2,2,1);plot(t,y_source);title('源信号');xlabel('时间/s');ylabel('振幅');
subplot(2,2,2);plot(t,y_whitenoise);title('白噪声');xlabel('时间/s');ylabel('振幅');
subplot(2,2,3);plot(t,y_mix);title('混合信号');xlabel('时间/s');ylabel('振幅');
subplot(2,2,4);plot(f,mag(1,1:N/2));title('傅里叶变换');xlabel('频率/Hz');ylabel('振幅');