gaosi.m

十几个数字信号处理实验课的程序

n=0:15;
x1=exp(-1*(n-8).^2/2);
x2=exp(-1*(n-8).^2/4);
x3=exp(-1*(n-8).^2/8);
figure
subplot(3,1,1);stem(x1);title('信号图形(p=8,q=2)');
subplot(3,1,2);stem(x2);title('信号图形(p=8,q=4)');
subplot(3,1,3);stem(x3);title('信号图形(p=8,q=8)');
figure
subplot(3,1,1);stem(abs(fft(x1)));title('信号频谱(p=8,q=2)');
subplot(3,1,2);stem(abs(fft(x2)));title('信号频谱(p=8,q=4)');
subplot(3,1,3);stem(abs(fft(x3)));title('信号频谱(p=8,q=8)');
x1=exp(-1*(n-8).^2/8);
x2=exp(-1*(n-13).^2/8);
x3=exp(-1*(n-14).^2/8);
figure
subplot(3,1,1);stem(x1);title('信号图形(p=8,q=8)');
subplot(3,1,2);stem(x2);title('信号图形(p=13,q=8)');
subplot(3,1,3);stem(x3);title('信号图形(p=14,q=8)');
figure
subplot(3,1,1);stem(abs(fft(x1)));title('信号频谱(p=8,q=8)');
subplot(3,1,2);stem(abs(fft(x2)));title('信号频谱(p=13,q=8)');
subplot(3,1,3);stem(abs(fft(x3)));title('信号频谱(p=14,q=8)');