matlab.rar_2ASK matlab_2ASK频谱_2ask_ask

%不同采样频率对实验结果的影响 fs1=1.0e+4;%设定采样频率 fs2=2.0e+4; fs3=3.0e+4; N=1024;%取的样本点 n=0:N-1; t1=n/fs1;%采样 t2=n/fs2; t3=n/fs3; f1=1.0e+3;%设定正弦信号频率 %生成不同采样频率正弦信号 x11=sin(2*pi*f1*t1); x12=sin(2*pi*f1*t2); x13=sin(2*pi*f1*t3); %生成不同采样频率下的输入信号（将高斯白噪声加入正弦信号） xa=awgn(x11,10,'measured');%加入高斯白噪声 xb=awgn(x12,10,'measured'); xc=awgn(x13,10,'measured'); %**************************************% %输入信号时域波形 figure(22);%时域波形 subplot(3,1,1); plot(t1,xa,'k');legend('不同采样频率下的时域波形f=10k'); grid on; hold on; subplot(3,1,2); plot(t2,xb,'r'); legend('不同采样频率下的时域波形f=20k'); grid on; hold on; subplot(3,1,3); plot(t3,xc); legend('不同采样频率下的时域波形f=30k'); grid on; hold on; %*************************************% %进行FFT变换并做频谱图 ya=fft(xa,N);%进行fft变换 yb=fft(xb,N); yc=fft(xc,N); magya=abs(ya);%求幅值 magyb=abs(yb); magyc=abs(yc); fya=(0:length(ya)-1)*2*pi/length(ya);%进行对应的频率转换 fyb=(0:length(yb)-1)*2*pi/length(yb); fyc=(0:length(yc)-1)*2*pi/length(yc); figure(23);%做频谱图 subplot(3,1,1); plot(fya,magya,'k'); axis([0 2*pi 0 600]); xlabel('频率(Hz)'); ylabel('幅值'); title('不同采样频率下输入信号的频谱图f=10k'); grid on; hold on; subplot(3,1,2); plot(fyb,magyb,'r'); axis([0 2*pi 0 600]); xlabel('频率(Hz)'); ylabel('幅值'); title('不同采样频率下输入信号的频谱图f=20k'); grid on; hold on; subplot(3,1,3); plot(fyc,magyc,''); axis([0 2*pi 0 600]); xlabel('频率(Hz)'); ylabel('幅值'); title('不同采样频率下输入信号的频谱图f=30k'); grid on; hold on;