MATLAB课时训练5

1.设计目的 

利用 MATLAB 实现信号的幅度调制， 绘制出原信号x(t) = sin(2πft)  以及已调信号y(t) = x(t) . cos(wt) 的时域波形图、频谱图以及功率谱。 将f 、w 作为参数变化，观察三种波形的变化情况。

2.设计题目涉及的理论知识

利用matlab绘制基本信号，绘制不同类型的信号特征图，并在改变参数的基础上，理解参数的意义。

（1）绘制时域波形图

（1）绘制时域波形图
源代码
 
 
subplot(221)
t=-1:0.01:5;                              %设置x的取值范围
f = 5;
y = sin(2*pi*t*f);                           
plot(t,y);
grid on;                                    %plot函数是绘制二维图形的函数
title('原信号');                              %title是指标题
xlabel('时间(s)','FontSize',10);                %xlabel是指横坐标
ylabel('幅度(V)','FontSize',10);                %ylabel是指纵坐标
axis([-1 5 -0.5 1.5])                       %axis是设置坐标轴的最大值，最小值
 
 
 
 
subplot(222)
t=-1:0.01:5;                              %设置x的取值范围
f  =  50;
y = sin(2*pi*t*f2);                           
plot(t,y);
grid on;                                    %plot函数是绘制二维图形的函数
title('原信号（2）');                              %title是指标题
xlabel('时间(s)','FontSize',10);                %xlabel是指横坐标
ylabel('幅度(V)','FontSize',10);                %ylabel是指纵坐标
axis([-1 5 -0.5 1.5])                       %axis是设置坐标轴的最大值，最小值
 
subplot(223)
t=-1:0.01:5;                              %设置x的取值范围
f =  5;
w = 5;
xt = sin(2*pi*t*f2);
y = xt.*cos(w*t);
plot(t,y);
grid on;                                %plot函数是绘制二维图形的函数
title('调频信号');                              %title是指标题
xlabel('时间(s)','FontSize',10);                %xlabel是指横坐标
ylabel('幅度(V)','FontSize',10);                %ylabel是指纵坐标
axis([-1 5 -0.5 1.5])                       %axis是设置坐标轴的最大值，最小值
 
subplot(224)
t=-1:0.01:5;                              %设置x的取值范围
f =  5;
w = 50;
xt = sin(2*pi*t*f2);
y = xt.*cos(w*t);
plot(t,y);
grid on;                                %plot函数是绘制二维图形的函数
title('调频信号（2）');                              %title是指标题
xlabel('时间(s)','FontSize',10);                %xlabel是指横坐标
ylabel('幅度(V)','FontSize',10);                %ylabel是指纵坐标
axis([-1 5 -0.5 1.5])                       %axis是设置坐标轴的最大值，最小值
截图

（2）频谱图；

（2）频谱图；
源代码
subplot(221)
t=-1:0.01:5;                              %设置x的取值范围
f0 = 5;
y0 = sin(2*pi*t*f0);                           
plot(t,y0);
grid on;                                    %plot函数是绘制二维图形的函数
title('原信号');                              %title是指标题
xlabel('时间(s)','FontSize',10);                %xlabel是指横坐标
ylabel('幅度(V)','FontSize',10);                %ylabel是指纵坐标
axis([-1 5 -0.5 1.5])                       %axis是设置坐标轴的最大值，最小值
 
subplot(222)
t=-1:0.01:5;                              %设置x的取值范围
f1  =  50;
y1 = sin(2*pi*t*f1);    
y0 = abs(fft(y1));
plot(t,y0);
grid on;                                    %plot函数是绘制二维图形的函数
title('频谱图');                              %title是指标题
 
subplot(223)
t=-1:0.01:5;                              %设置x的取值范围
f3 =  5;
w1 = 50;
xt = sin(2*pi*t*f3);
y1 = xt.*cos(w1*t);
plot(t,y1);
grid on;                                %plot函数是绘制二维图形的函数
title('调频信号（2）');                              %title是指标题
xlabel('时间(s)','FontSize',10);                %xlabel是指横坐标
ylabel('幅度(V)','FontSize',10);                %ylabel是指纵坐标
axis([-1 5 -0.5 1.5])                       %axis是设置坐标轴的最大值，最小值
 
 
subplot(224)
t=-1:0.01:5;                              %设置x的取值范围
f2 =  50;
w0 = 5;
xt = sin(2*pi*t*f2);
y2 = xt.*cos(w0*t);
y0 = abs(fft(y2));
plot(t,y0);
grid on;                                %plot函数是绘制二维图形的函数
title('频谱图');                              %title是指标题
 

（3）功率图；

（3）功率图；
源代码
subplot(221)
t=-1:0.01:5;                              %设置x的取值范围
f0 = 5;
y0 = sin(2*pi*t*f0);                           
plot(t,y0);
grid on;                                    %plot函数是绘制二维图形的函数
title('原信号');                              %title是指标题
xlabel('时间(s)','FontSize',10);                %xlabel是指横坐标
ylabel('幅度(V)','FontSize',10);                %ylabel是指纵坐标
axis([-1 5 -0.5 1.5])                       %axis是设置坐标轴的最大值，最小值
 
subplot(222)
t=-1:0.01:5;                              %设置x的取值范围
f1  =  50;
y1 = sin(2*pi*t*f1);    
fs=10;                      %采样频率
NFFT=50;                    %即FFT数据点的个数
window=50;                  %window 是计算功率谱每个窗口的信号长度，
noverlap=30;                %noverlap 是每个窗口之间重叠的长度
[pxx,f] = pwelch(y1,window,noverlap,NFFT,fs);%调用pwelch函数
plot(f,10*log10(pxx))
xlabel('Frequency (Hz)')
ylabel('PSD (dB/Hz)')
grid on;                                    %plot函数是绘制二维图形的函数
title('功率图');                              %title是指标题
 
subplot(223)
t=-1:0.01:5;                              %设置x的取值范围
f3 =  5;
w1 = 50;
xt = sin(2*pi*t*f3);
y1 = xt.*cos(w1*t);
plot(t,y1);
grid on;                                %plot函数是绘制二维图形的函数
title('调频信号（2）');                              %title是指标题
xlabel('时间(s)','FontSize',10);                %xlabel是指横坐标
ylabel('幅度(V)','FontSize',10);                %ylabel是指纵坐标
axis([-1 5 -0.5 1.5])                       %axis是设置坐标轴的最大值，最小值
 
 
subplot(224)
t=-1:0.01:5;                              %设置x的取值范围
f2 =  50;
w0 = 5;
xt = sin(2*pi*t*f2);
y2 = xt.*cos(w0*t);
y0 = abs(fft(y2));
fs=10;                             %采样频率
NFFT=50;                           %即FFT数据点的个数
window=50;                         %window 是计算功率谱每个窗口的信号长度，
noverlap=30;                         %noverlap 是每个窗口之间重叠的长度
[pxx,f] = pwelch(y0,window,noverlap,NFFT,fs);  %调用pwelch函数
plot(f,10*log10(pxx))
xlabel('Frequency (Hz)')
ylabel('PSD (dB/Hz)')
grid on;                                %plot函数是绘制二维图形的函数
title('调频信号的功率图');                              %title是指标题

三  结论（结果分析）

在改变参数的情况下，对于信号会有明显的变化。由我所作来看，对于f而言，f的变化对于图像的影响似乎不大，对于w而言，w越大，图像越密集。

四  心得体会

ATLAB的内部函数库提供了相当丰富的函数，这些函数可以解决许多基本问题，如矩阵的输入。在其他语言中（如C语言），要输入一个矩阵，先要编写一个矩阵的子函数，而MATLAB语言则提供了一个人机交互的数学系统环境，该系统的基本数据结构是矩阵，在生成矩阵对象的时候，不要求做明确的维数说明。与利用C语言或FORTRAN语言编写数值计算的程序设计相比，利用MATLAB可以节省大量的编程时间。这给用户节省很多的时间，使用户能够把自己的精力放在创造方面，而把繁琐的问题交给内部函数来解决。