内容发布更新时间 : 2024/12/23 8:32:47星期一 下面是文章的全部内容请认真阅读。
Q1-4已知一个序列为
x[n]???n,?0,0?n?4otherwise
编写MATLAB程序Q1_4,能够将x[n]以N = 8为周期进行周期延拓得到一个周期为N =8的周期序列y[n],并分别绘制x[n]和y[n]图形。
编写的程序Q1_4如下:
clear, close all; n1=-3:7;
x=n1.*(u(n1)-u(n1-5));
n2=-16:16; N=8; y=0; for k=-2:2;
y=y+(n2-k*N).*(u(n2-k*N)-u(n2-5-k*N)); end
subplot(211) stem(n1,x,'.')
title ('Original signal x[n]') xlabel ('Time n (sec)')
subplot(212) stem(n2,y,'.')
title ('Original signal y[n]') xlabel ('Time n (sec)')
信号x[n]的波形图
信号y[n]的波形图
Original signal x[n]43210-3-2-1023Time n (sec)Original signal y[n]1456743210-20-15-10-50Time n (sec)5101520
Q1-5 仿照范例程序Program1_7,编写程序Q1_5,计算并绘制由如下微分方程表示的系统
在输入信号为x(t) = (e-2t - e-3t)u(t)时的零状态响应和你手工计算得到的系统零状态响应曲线。
d2y(t)dy(t)?3?2y(t)?8x(t) 2dtdt手工计算得到的系统零状态响应的数学表达式是:y?(4e?t?8te?2t?4e?3t)u(t) 编写的程序Q1_5如下:
% Program1_7
% This program is used to compute the impulse response h(t) and the step response s(t) of a
% continuous-time LTI system clear, close all,
num = input('微分方程右边的系数:'); % 8 den = input('微分方程左边的系数:'); % [1 3 2]
t = 0:0.01:8;
x = input('表达式的输入信号x(t):'); %input (exp(-2*t)-exp(-3*t)).*u(t)
subplot(221);
impulse(num,den,8);
title ('单位冲激响应 x(t)') xlabel ('Time t (sec)')
subplot(222);
step(num,den,8); title ('单位阶跃响应 x(t)') xlabel ('Time t (sec)')
subplot(223);
plot(t,x); title ('零输入响应 x(t)') xlabel ('Time t (sec)')
subplot(224);
y=lsim(num,den,x,t); plot(t,y);
title ('零状态响应 x(t)') xlabel ('Time t (sec)')
手工计算得到的系统零状态响应曲线代码: clear, close all, dt=0.01; t=-5:dt:10;
y=(4*exp(-t)-8*t.*exp(-2*t)-4*exp(-3*t)).*u(t); plot(t,y);
axis([-2,10,0,0.3])
title ('微分方程解的图像 y(t)'); xlabel ('Time index t');
用MATLAB绘制的手工计算的系统响应
微分方程解的图像 y(t)0.250.20.150.10.050-2024Time index t6810
执行程序Q1_5得到的系统响应
单位冲激响应 x(t)21.543单位阶跃响应 x(t)Amplitude10.50Amplitude0246821002468Time t (sec) (sec)Time t (sec) (sec)零输入响应 x(t)0.20.150.10.0500246Time t (sec)80.40.30.20.100零状态响应 x(t)246Time t (sec)8
四、实验结论与体会
在这第一次实验中,由于对实验的不了解,导致做实验的时候不知道该怎么样去做,还是经过了一段时间的熟悉,才逐渐的明白要做什么,但是已经浪费了大量时间,使得原来的任务没有做完,只好推迟做事之前应当把准备工作做好、
实验二连续时间信号的频域分析
实验名称 信号与系统的频域分析 评分
实验日期 2015 年 12 月 22 日 指导教师 张鏊峰 姓名 易煌轩 专业班级 通信1401 学号 201403020130
一、实验目的
1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法; 2、观察截短傅里叶级数而产生的“Gibbs现象”,了解其特点以及产生的原因; 3、掌握连续时间傅里叶变换的分析方法及其物理意义;
4、掌握各种典型的连续时间非周期信号的频谱特征以及傅里叶变换的主要性质;
基本要求:掌握并深刻理傅里叶变换的物理意义,掌握信号的傅里叶变换的计算方法,掌握利用MATLAB编程完成相关的傅里叶变换的计算。
二、实验要求
学会利用MATLAB求解信号与系统的相关问题; 三、实验内容
Q2-1 编写程序Q2_1,绘制下面的信号的波形图:
?1n?11?sin()cosn(?0t) x(t)?cos?(t)?cos3?(t)?cos5?(t)???000n235n?1其中,?0 = 0.5π,要求将一个图形窗口分割成四个子图,分别绘制cos(?0t)、
cos(3?0t)、cos(5?0t) 和x(t) 的波形图,给图形加title,网格线和x坐标标签,并且程序能够接受从键盘输入的和式中的项数。 抄写程序Q2_1如下:
clear,
close all, dt = 0.001; t = -5:dt:5; w0 = 0.5*pi;
i = input('请输入和式中的项数:');
x = 0;
for r = 1:i;
x =
x+1/(r*2-1)*sin((2*r-1)*pi/2)*cos((2*r-1)*w0*t);
end;