1、 编写一段程序能够产生实指数信号 x(t)=e 。要求:t 的取值范围是从-5 到5,步长为 0.
01,使用函数 axis()控制图形的时间范围在0~3 秒之间。然后执行该程序,画出所得图形。
clearclear all variables
close allclose all figure windows
dt = 0.01specify the step of time variable
t = 5:dt:5specify the interval of time
x = exp(-2*tgenerate the signal
plot(t,x)
title('sinusoidal signal x(t)')
xlabel('time t ')
grid on
axis([0,3,0,3])%open a figure window and draw the plot of x(t)
2、编写一段程序,使之能够在同一个图形窗口中的两个子图中分别绘制信号 x(t)=0.5 和x(t)=cos(πt)[u(t)-u(t-2)]。要求:
t 的取值范围是从-5 到 5,步长为 0.01。
答:3、编写一段程序,由给定信号 x(t) =u(t),求信号 y(t) =x(1.5t+3),并在同一个窗口中绘制出 x(t) 和 y(t)的图形。
要求:t 的取值范围是从-5到5,步长为0.01。
信号 x(t) =e-0.5tu(t);求信号y(t) =x(1.5t+3),并绘制出x(t) 和y(t)的图形。
第一种:2、 4编写程序计算并绘制由如下微分方程表示的系统在输入信号为x(t) =e-e)u(t)时的零状态响应曲线。
5、编写一段程序,绘制下面的信号的波形图:x(t)=cos(t)- cos(3t)+ cos(5t)-…其中= 0.5π,要求绘出x(t)的波形图,其中x(t)为n=60项的和。
第一种:6、给定如下周期信号:周期为2
编写一段程序,以计算x(t)的傅里叶级数的系数,并求出n=5 时的11 个系数。
clear, close all
t = 2; dt = 0.001; t = 4:dt:4;
f1 = u(t+0.2) -u(t-0.2-dt);
f = 0;
for m = 1:1periodically extend f1(t) to form a periodic signal
f= f + u(t+0.2-m*t) -u(t-0.2-m*t-dt);
endw0 = 2*pi/t;
n = 5the number of the harmonic components
l = 2*n+1;
for n = n: nevaluate the fourier series coefficients fn
fn(n+1+n) =1/t)*f1*exp(-j*n*w0*t')*dt;
end结果:
在默认主页输入 fn 然后按回车键。
fn =columns 1 through 3
-0.0005 - 0.0000i 0.0464 + 0.0000i 0.1008 - 0.0000i
columns 4 through 6
0.1515 + 0.0000i 0.18750.2005
columns 7 through 9
0.18750.1515 - 0.0000i 0.1008 + 0.0000i
columns 10 through 11
0.0464 - 0.0000i -0.0005 + 0.0000i
7、给定时限信号:x(t)=cos()[u(t+1)-u(t-1)],编写matlab 程序,绘制出信号x(t)的幅度谱和相位谱。
先建立m文件,命名为。
function y = u(t)
y = t>=0); y = 1 for t > 0, else y = 0
在建立一个m文件,写以下程序。
t = 0.01; dw = 0.1; %时间和频率变化的步长。
t = 10:t:10;
w = 4*pi:dw:4*pi;
x=cos(pi*t/2).*u(t+1)-u(t-1));
x=x*exp(-j*t'*w)*t; %傅里叶变换。
x1=abs(x计算幅度谱。
phai=angle(x计算相位谱。
figure
plot(t,x)
figure
plot(w,x1)
title('the amplitude |x| of x(t)')
figure
plot(w,phai)
title('the phase phi(x) of x(t)')
8、给定时限信号:x ( t ) u (t + 2 ) u (t - 2 ),编写matlab 程序,绘制出信号x(t)的幅度谱和相位谱。
答:t = 0.01; dw = 0.1; %时间和频率变化的步长。
t = 10:t:10;
w = 4*pi:dw:4*pi;
x=u(t+2)-u(t-2);
x=x*exp(-j*t'*w)*t; %傅里叶变换。
x1=abs(x计算幅度谱。
phai=angle(x计算相位谱。
figure
plot(w,x1)
figure
plot(w,phai)
凡是题目中有u(t)的都要先建立一个名为的文件,再做题)是默认的文件名,不要更改。
9、某因果系统的系统函数为h(s)=绘出系统的极零点图。编写程序如下(绘制极零点图的子函数也要写出来)
答:先建立一个m文件。
function splane(num,den)
p = roots(dendetermine the poles
q = roots(numdetermine the zeros
p = p'; q = q';
x = max(abs([p qdetermine the range of real-axis
x = x+1;
y = xdetermine the range of imaginary-axis
plot([-x x],[0 0],'hold ondraw the real-axis
plot([0 0],[y y],'hold ondraw the imaginary-axis
plot(real(p),imag(p),'x');hold ondraw the poles
plot(real(q),imag(q),'o');hold ondraw the zeros
title('zero-pole plot');
xlabel('real part');ylabel('imaginal part')
axis([-x x -y y]);
保存上面那个程序后再建立一个m文件运行下面的。
clear, close all,num=[3 0] 可以采用键盘输入num=input(‘num=’)
den=[1 5 6den=input(‘den=’)
信号上机实验
信号上机实验三。一 实验目的。1 掌握信号时域运算的matlab实现方法。2 掌握信号时域变换的matlab实现方法。二 实验原理。信号的时域运算包括信号的相加 相乘,信号的时域变换包括信号的平移 反折 倒相及信号的尺度变换。三 实验内容。1 将实验原理中提到的例子在计算机上全部练习一遍 解 例3 ...
数字信号上机作业
数字信号处理上机报告。一 实验内容 1 利用傅立叶级数展开的方法,自由生成所需的x t 2 通过选择不同的采样间隔t 分别选t 或 1 2fc 从x t 获得相应的x n 作出x n 图形 3 对获得的不同x n 分别作傅立叶变换,分析其频率响应特性 给出幅频与相频特性曲线 4 利用巴特沃思 切比雪...
数字信号题和答案
填空。1线性时不变系统的表示方法 常系数线性差分方程 卷积 系统函数。2 线性时不变系统是因果系统的充要条件 h n 0,n 0 3离散时域求解法 迭代法 卷积计算法。4由u n 如何能得到u n 2 5 fft算法有俩种 按时间轴抽取法 dit 按频率抽取法 dif 6 级数收敛的充要条件是 满足...