已知信号 x(t) =cos(t)试画出采样频率分别为 1200Hz、 800Hz 和 600Hz 时， 信号采样值的波形及其频 谱 ,matlab

下面是在Matlab中画出采样频率分别为1200Hz、800Hz和600Hz时，信号采样值的波形及其频谱的代码和结果。

% 采样频率
fs1 = 1200; % 采样频率为1200Hz
fs2 = 800; % 采样频率为800Hz
fs3 = 600; % 采样频率为600Hz

% 时间序列
t = 0:1/(2*pi*fs1):2*pi; % 采样频率为1200Hz时的时间序列
t2 = 0:1/(2*pi*fs2):2*pi; % 采样频率为800Hz时的时间序列
t3 = 0:1/(2*pi*fs3):2*pi; % 采样频率为600Hz时的时间序列

% 信号
x = cos(t); % 原始信号

% 采样
xn1 = cos(2*pi*fs1*t); % 采样频率为1200Hz时的采样序列
xn2 = cos(2*pi*fs2*t2); % 采样频率为800Hz时的采样序列
xn3 = cos(2*pi*fs3*t3); % 采样频率为600Hz时的采样序列

% 绘制波形
subplot(3,1,1)
plot(t,x)
hold on
stem(t,xn1)
title('采样频率为1200Hz时的信号波形')
legend('原始信号','采样信号')
xlabel('时间/s')
ylabel('幅值')

subplot(3,1,2)
plot(t,x)
hold on
stem(t2,xn2)
title('采样频率为800Hz时的信号波形')
legend('原始信号','采样信号')
xlabel('时间/s')
ylabel('幅值')

subplot(3,1,3)
plot(t,x)
hold on
stem(t3,xn3)
title('采样频率为600Hz时的信号波形')
legend('原始信号','采样信号')
xlabel('时间/s')
ylabel('幅值')

% 绘制频谱
N = length(t);
f = -fs1/2:fs1/N:fs1/2-fs1/N; % 采样频率为1200Hz时的频率序列
f2 = -fs2/2:fs2/N:fs2/2-fs2/N; % 采样频率为800Hz时的频率序列
f3 = -fs3/2:fs3/N:fs3/2-fs3/N; % 采样频率为600Hz时的频率序列

X = fftshift(fft(x));
Xn1 = fftshift(fft(xn1));
Xn2 = fftshift(fft(xn2));
Xn3 = fftshift(fft(xn3));

figure
subplot(3,1,1)
plot(f,abs(X))
hold on
stem(f,abs(Xn1))
title('采样频率为1200Hz时的信号频谱')
legend('原始信号','采样信号')
xlabel('频率/Hz')
ylabel('幅值')

subplot(3,1,2)
plot(f,abs(X))
hold on
stem(f2,abs(Xn2))
title('采样频率为800Hz时的信号频谱')
legend('原始信号','采样信号')
xlabel('频率/Hz')
ylabel('幅值')

subplot(3,1,3)
plot(f,abs(X))
hold on
stem(f3,abs(Xn3))
title('采样频率为600Hz时的信号频谱')
legend('原始信号','采样信号')
xlabel('频率/Hz')
ylabel('幅值')

下面是在Matlab中绘制出采样频率分别为1200Hz、800Hz和600Hz时，信号采样值的波形及其频谱的结果。