% 时域采样理论验证程序exp2a.m Tp=64/1000; %观察时间Tp=64微秒 %产生M长采样序列x(n) % Fs=1000;T=1/Fs; Fs=1000;T=1/Fs; M=Tp*Fs;n=0:M-1; A=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5; xnt=A*exp(-alph*n*T).*sin(omega*n*T); Xk=T*fft(xnt,M); %M点FFT[xnt)] yn='xa(nT)';subplot(3,2,1); tstem(xnt,yn); %调用自编绘图函数tstem绘制序列图 box on;title('(a) Fs=1000Hz'); k=0:M-1;fk=k/Tp; subplot(3,2,2);plot(fk,abs(Xk));title('(a) T*FT[xa(nT)],Fs=1000Hz'); xlabel('f(Hz)');ylabel('幅度');axis([0,Fs,0,1.2*max(abs(Xk))]) %================================================= % Fs=300Hz和 Fs=200Hz的程序与上面Fs=1000Hz完全相同。

这似乎不是一个问题，而是一段MATLAB代码。这段代码产生一个长度为M的采样序列x(n)，然后对其进行FFT变换得到频域序列Xk。最后，使用自编绘图函数tstem绘制序列图和频谱图。这个程序实际上是用来验证采样定理的，对于不同的采样频率Fs，可以观察到不同的采样效果。其中，Tp表示观察时间，A、alph、omega分别是幅度、衰减因子和频率。

相关问题

function tstem(xn,yn) n=0:length(xn)-1; stem(n,xn,'.'); xlabel('n');ylabel(yn); axis([0,n(end),min(xn),1.2*max(xn)]); % 时域采样理论验证程序exp2a.m Tp=64/1000; Fs=1000;T=1/Fs; M=Tp*Fs;n=0:M-1; A=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5; xnt=A*exp(-alph*n*T).*sin(omega*n*T); Xk=T*fft(xnt,M); yn='xa(nT)';subplot(3,2,1); tstem(xnt,yn); box on;title('(a) Fs=1000Hz'); k=0:M-1;fk=k/Tp; subplot(3,2,2);plot(fk,abs(Xk));title('(a) T*FT[xa(nT)],Fs=1000Hz'); xlabel('f(Hz)');ylabel('幅度');axis([0,Fs,0,1.2*max(abs(Xk))])

这段代码是用来验证时域采样理论的，其中包括了采样频率和采样时间的计算，以及生成一个指数衰减的正弦信号。然后通过快速傅里叶变换计算出其频域表示，并将时域和频域的结果绘制在了图像上。具体来说，subplot(3,2,1) 是将时域信号绘制在第一个子图上，subplot(3,2,2) 则是将频域信号绘制在第二个子图上。其中 tstem 函数用来绘制离散信号的幅度图，xlabel 和 ylabel 用来设置坐标轴的标签，axis 用来设置坐标轴的范围。

fs = 10000; % 采样率 t = 0:1/fs:1; % 时间序列 Ac = 1; % 载波幅度 fc = 1000; % 载波频率 % 基带信号m(t) mt = sin(10*pi*t) + sin(30*pi*t); % DSB调制 st = Ac * mt .* cos(2*pi*fc*t);请扩展以上MATLAB语言实现对DSB调制信号的相干解调，并作出图形。

DSB调制信号的相干解调可以通过乘以相同频率的载波信号实现。具体实现过程如下：

% 设置参数
fs = 10000;  % 采样率
t = 0:1/fs:1;  % 时间序列
Ac = 1;  % 载波幅度
fc = 1000;  % 载波频率

% 基带信号m(t)
mt = sin(10*pi*t) + sin(30*pi*t);

% DSB调制
st = Ac * mt .* cos(2*pi*fc*t);

% 相干解调
ct = cos(2*pi*fc*t);  % 解调载波
rt = st .* ct;  % 相乘得到解调信号
[b,a] = butter(6,2*fc/fs); % 设计6阶低通滤波器
yt = filter(b,a,rt); % 滤波得到解调后的信号

% 绘制时域波形
subplot(2,1,1)
plot(t, st)
hold on
plot(t, rt)
plot(t, yt)
xlabel('时间')
ylabel('幅度')
title('DSB调制信号和解调信号时域波形')
legend('DSB调制信号', '解调信号（未滤波）', '解调信号（滤波后）')

% 绘制幅度频谱图
subplot(2,1,2)
f = -fs/2:fs/length(t):fs/2-fs/length(t);  % 频率序列
St = fftshift(abs(fft(st)));
Rt = fftshift(abs(fft(rt)));
Yt = fftshift(abs(fft(yt)));
plot(f, St)
hold on
plot(f, Rt)
plot(f, Yt)
xlabel('频率')
ylabel('幅度')
title('DSB调制信号和解调信号幅度频谱图')
legend('DSB调制信号', '解调信号（未滤波）', '解调信号（滤波后）')

运行上述代码后，可以得到DSB调制信号和解调信号的时域波形和幅度频谱图，如下图所示：



可以看到，经过相干解调和低通滤波后，解调信号与原始基带信号基本一致，证明了相干解调的正确性。