信号系统实验报告
y=step(sys,t);
plot(t,y);
xlabel('Time(sec)')
ylabel('h(t)')
运行结果如有图所示:
3-4程序如下:
%program3_4
x=[0.85,0.53,0.21,0.67,0.84,0.12];
y=[0.68,0.37,0.83,0.52,0.71];
z=conv(x,y);
title('x2[k]')
A=1;w3=26*pi;phi=0;
t3=0:0.0001:1;
xt3=A*cos(w3*t3+phi);
subplot(2,3,3)
plot(t3,xt3)
title('x3(t)')
k3=0:20;
xk3=cos(13*pi/5*k3);
subplot(2,3,6)
stem(k3,xk3)
title('x3[k]')
运行结果如右上图:
2-8分别用square函数和
Swtooth画出题M2-8图所
示的离散周期序列。
程序如下:
%program2-81
k=0:1:40;
A=1;T=20;w0=2*pi/T;
xk=A*sawtooth(w0*k,0.5);
stem(k,xk)
ts=0;te=10;dt=0.001;
sys=tf([2 1],[1 3 2]);
t=ts:dt:te;
x=1*exp(-3*t)*1;
y=lsim(sys,x,t);
plot(t,y);
xlabel('Time(sec)')
ylabel('y(t)')
3-2程序如下:
%program3-2
ts=0;te=5;dt=0.01;
%example2_2
t=-2:0.001:2;
subplot(2,1,1)
plot(t,x2_2(t))
title('x(t)')
subplot(2,1,2)
A=1;w0=50;
t=-2:0.001:2;phi=0;
xt=[A*cos(w0*t+phi)].*[1*(t>=-1&t<0)+4*t.*(t>=0&t<=0.5)+(-4)*(t-1).*(t>0.5&t<=1)];
信号系统实验报告
———————————————————————————————— 作者:
———————————————————————————————— 日期:
ﻩ
信
号
实
验
报
告
信息科学与工程学院
电子信息工程
姓 名 :
学 号 :
实验一
一、实验内容:
1、连续信号的MATLAB表示:
MATLAB提供了大量的产生基本信号的函数。最常用的指数信号、正弦信号是MATLAB的内部函数,即不安装任何工具箱就可调用的函数。
title('x1[k]')
A=1;w2=14*pi;phi=0;
t2=0:0.0001:1;
xt2=A*cos(w2*t2+phi);
subplot(2,3,2)
plot(t2,xt2)
title('x2(t)')
k2=0:20;
xk2=cos(7*pi/5*k2);
subplot(2,3,5)
stem(k2,xk2)
n=length(z);
stem(-3:n-4,z);
运行结果如有图所示:
3-7某离散时间LTI系统
满足差分方程
Y[k]+0.42y[k-1]-0.19y
[k-2]=0.31x[k]+0.68x
[k-1],若x(t)=0.6ku(k)
%example2_3
t=-4:0.01:10;
subplot(3,1,1)
plot(t,x2_3(t))
title('x(t)')
axis([-4,10,-2,3])
t=-4:0.01:10;
subplot(3,1,2)
plot(t,x2_3(0.5*t))
title('x(0.5t)')
axis([-4,10,-2,3])
运行结果如右图所示:
%program2-82
k=0:1:40;
A=-1;T=10;w0=2*pi/T;
xk=A*square(w0*k,50);
stem(k,xk)
运行结果如右图:
实验二
3-1一个连续时间LTI
系统满足的微分方程为:
y″(t)+3y′(t)+2y(t)
=2x′(t)+x(t)
%program3_1
sys=tf([6],[1 3]);
t=ts:dt:te;
y=impulse(sys,t);
plot(t,y);
xlabel('Time(sec)')
ylabel('h(t)')
运行结果如有图所示:
%program3-2
ts=0;te=5;dt=0.01;
sys=tf([6],[1 3]);
t=ts:dt:te;
plot(t,xt)
title('x(t)cos(50Байду номын сангаас)')
运行结果如图所示:
2-3(1)编写表示题M2-3图所示的信号波形x(t)的MATLAB函数(2)试画出x(t)、x(0.5t)、和x(2-0.5t)的波形图。程序如下:
function yt=x2_3(t)
yt=1*t.*(t>=0&t<2)+2*(t>=2&t<3)+(-1)*(t>=3&t<5);
2、编程实现题
2-2已知信号x(t)如图m2-2图所示,分别用MATLAB表示信号x(t)、x(t)cos(50t),并画出波形,程序如下:
function yt=x2_2(t)
yt=1*(t>=-1&t<0)+4*t.*(t>=0&t<=0.5)+(-4)*(t-1).*(t>0.5&t<=1);
%program2-7
A=1;w1=6*pi;phi=0;
t1=0:0.0001:1;
xt1=A*cos(w1*t1+phi);
subplot(2,3,1)
plot(t1,xt1)
title('x1(t)')
k1=0:20;
xk1=cos(3*pi/5*k1);
subplot(2,3,4)
stem(k1,xk1)
t=-6:0.01:8;
subplot(3,1,3)
plot(t,x2_3(2-0.5*t))
title('x(2-0.5t)')
axis([-6,8,-2,3]
运行结果如有右图:
2-4 画出题M2-4图所示信号的奇分量和偶分量。程序如下:
2-7已知连续信号x1(t)=cos(6∏t)、x2(t)=cos(14∏t)、x3(t)=cos(26∏t),以抽样函数fs=10Hz对这三个信号进行抽样得离散序列x1[k]、x2[k]、x3[k].试在同一图上画出连续信号和其他对应的离散序列。程序如下: