设有三个不同频率的正弦信号,频率分别为f1=100hz,f2=200hz,f3=3800hz.设有三个不同频率的正弦信号,频率分别为f1=100hz,f2=200hz,f3=3800hz.现在用抽样频率f3=3800hz对这三个正弦信号进行抽样,用MATLAB命令

设有三个不同频率的正弦信号,频率分别为f1=100hz,f2=200hz,f3=3800hz.设有三个不同频率的正弦信号,频率分别为f1=100hz,f2=200hz,f3=3800hz.现在用抽样频率f3=3800hz对这三个正弦信号进行抽样,用MATLAB命令
设有三个不同频率的正弦信号,频率分别为f1=100hz,f2=200hz,f3=3800hz.
设有三个不同频率的正弦信号,频率分别为f1=100hz,f2=200hz,f3=3800hz
.现在用抽样频率f3=3800hz对这三个正弦信号进行抽样,用MATLAB命令画出各抽样信号的波形及频谱,并分析频率混叠现象.

设有三个不同频率的正弦信号,频率分别为f1=100hz,f2=200hz,f3=3800hz.设有三个不同频率的正弦信号,频率分别为f1=100hz,f2=200hz,f3=3800hz.现在用抽样频率f3=3800hz对这三个正弦信号进行抽样,用MATLAB命令
clf;
fs=3800;N=38;   %采样频率和数据点数
n=0:N-1;t=n/fs;   %时间序列
x1=sin(2*pi*100*t); 
x2=sin(2*pi*200*t); 
x3=sin(2*pi*3800*t); 
y1=fft(x1,N);    %对信号进行快速Fourier变换
mag1=abs(y1);     %求得Fourier变换后的振幅
y2=fft(x2,N);    %对信号进行快速Fourier变换
mag2=abs(y2);     %求得Fourier变换后的振幅
y3=fft(x3,N);    %对信号进行快速Fourier变换
mag3=abs(y3);     %求得Fourier变换后的振幅
f=n*fs/N;    %频率序列
subplot(3,1,1),plot(f(1:N/2),mag1(1:N/2)); %绘出Nyquist频率之前随频率变化的振幅
xlabel('频率/Hz');
ylabel('振幅');title('f=100Hz');grid on;
subplot(3,1,2),plot(f(1:N/2),mag2(1:N/2)); %绘出Nyquist频率之前随频率变化的振幅
xlabel('频率/Hz');
ylabel('振幅');title('f=200Hz');grid on;
subplot(3,1,3),plot(f(1:N/2),mag3(1:N/2)); %绘出Nyquist频率之前随频率变化的振幅
xlabel('频率/Hz');
ylabel('振幅');title('f=3800Hz');grid on;
fft频谱图如下：