目录。数字信号处理第一次上机实验报告 1
一、实验题目 1
二、实验目的 1
三、程序 1
四、运行结果 3
五、 结论 4
姓名:石星宇学号:02123010
1、考虑由连续时间信号,求采样后的离散时间信号及量化后的数字信号的matlab实现问题;
2、输出以下三种信号:
2 对采样,
3 对采样,
3、编写程序计算两序列线性卷积(不允许使用matlab自带的卷积函数)
1、通过实验熟悉matlab在数字信号处理中的应用。
2、熟悉数字信号处理中常用的函数以及语句。
3、加深对采样定理的理解。
4、加深对线性卷积计算过程的理解。
1、连续时间信号,求采样后的离散时间信号及量化后的数字信号的matlab实现。
设连续时间信号为。
采样时间间隔为,量化步长为。
程序如下:dt=0.001;tf=6;t=0:dt:tf;
xa=sqrt(t)+cos(t);
t=0.5;n=0:tf/t;
x=sqrt(n*t)+cos(n*t);
deltax=0.5;
xq=round(x/deltax)*deltax;
subplot(121)
plot(t,xa,':hold on ,grid on
plot(n*t,x,'o')
stem(n*t,xq,'*grid on
legend('连续时间信号xa','离散时间信号x','数字信号xq')
subplot(122)
stairs(n*t,xq),grid on
legend('数字信号采样保持')
set(gcf,'color','w')
2、输出给定的三种信号。
t=linspace(0,10,5000);%连续自变量。
tf=10;
xa=sin(2.5*pi*t);%原始连续时间信号xa
fs1=2.5;fs2=4;%两个采样率。
n1=0:tf*fs1; n2=0:tf*fs2;%离散自变量向量。
x1=sin(2.5*pi*n1*(1/fs1));采样率为2.5hz的离散时间信号x1
x2=sin(2.5*pi*n2*(1/fs2));采样率为4hz的离散时间信号x2
figure(1)
plot(t,xa,':hold on ,grid on
stem(n1/fs1,x1,'ro')
stem(n2/fs2,x2,'m*')
legend('连续时间信号xa','采样信号x1(fs=2.5hz)',采样信号x2(fs=4hz)')
title('连续信号及其采样信号')
xlabel('t'),ylabel('xa(t)、x1(n)、x2(n)')
3、编写程序计算两序列线性卷积。
function w=my_conv(x,y)
len=length(x)+length(y)-1;
fn=zeros(1,len);
gn=zeros(1,len);
fn(length(y):end)=x;
for m=1:length(y)
g1(m)=y(length(y)-m+1);
end gn(1:length(y))=g1;
for n=1:len
gn1=circshift(gn,[0,n-1]);
w(n)=sum(fn.*gn1);
end1、连续时间信号、离散时间信号及量化后的数字信号如右图所示。
2、三种给定的信号输出如下图所示。
3、计算两序列的线性卷积。
设输入序列如下:
用编写的程序计算结果如下,用matlab自带的卷积函数计算,并与之对比,结果是正确的。
计算机无法严格地表示连续量和模拟量,就自变量来说,只能把取的很密。将自变量时间轴按采样周期进行分割,得到离散时间自变量。通常从零开始,取顺序增加的整数,直到终点,matlab语言写成。
数字信号处理上机作业
目录。数字信号处理第二次上机实验报告 3 实验目的 3 实验一 3一 实验题目 3 二 程序 3 三 运行结果 4 1 序列长度为n 16时 4 2 序列长度为n 32时 5 四 结论 7 实验二 7一 实验题目 7 二 程序 7 1 计算及画图函数 7 2 序列输入脚本 8 三 运行结果 9 1 ...
数字信号上机作业
数字信号处理上机报告。一 实验内容 1 利用傅立叶级数展开的方法,自由生成所需的x t 2 通过选择不同的采样间隔t 分别选t 或 1 2fc 从x t 获得相应的x n 作出x n 图形 3 对获得的不同x n 分别作傅立叶变换,分析其频率响应特性 给出幅频与相频特性曲线 4 利用巴特沃思 切比雪...
数字信号处理作业
题一 令,其中,u n 是白噪声。1 使用matlab中的有关文件产生均值为0,功率为0.1的均匀分布白u n 求其自相关函数并画出其波形。解 1 产生均值为0,功率为0.1的均匀分布的白u n 的程序如下 答 1 matlab的rand函数可用来产生均值微0.5,方差 即功率 为1 12,在 0,...