有偿，用matlab,救救孩子吧

一个理想采样及恢复系统，采样频率为Ws=6π rad/s,采样后经过一个带宽
为3π,增益为1/3的理想低通还原。现有输入
xa(t)=cost+cos2πt+cos5,求输出信号y(t)并画出输入与输出信号时
域及频域图形，频谱图要求横坐标以Hz为单位。

你好！代码供参考：

%输入信号
figure(1);clf
Fs = 6;
t = 0:1/Fs:100;
xa = cos(t) + cos(2*pi*t) + cos(5);
y = 1/3*lowpass(xa, (3*pi)/(6*pi), Fs);
plot(t,xa,'r-',t,y,'b-')
xlabel('t');ylabel('y');legend('输入','输出')

%% 绘制频谱图
n = 2^(floor(log2(numel(y))));
t = t(1:n);
x = xa(1:n);
L = numel(t);
x = fft(x,n);
P2 = abs(x/L);
k = floor(L/2);
f = (0:k-1)*Fs/L;
figure(20);clf 
plot(f(1:k),(P2(1:k)),'r-');hold on%输入信号
x = y(1:n);
y = fft(y,n);
P2 = abs(y/L);
k = floor(L/2);
f = (0:k-1)*Fs/L;
plot(f(1:k),(P2(1:k)),'b-');hold on%输出信号
legend('输入','输出')
xlabel('频率(Hz)');ylabel('幅值')