- 1、本文档共10页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
1-1-1
clear
close all
dt=0.01;( dt=0.2;)
t=-2:dt:2;
x=sin(2*pi*t);
plot(t,x)
title('sinusoidal singnalx(t)')
xlabel('Time t(sec)')
1-2
t=-5:0.001:5;
x=exp(-0.5*t).*u(t);
y=exp(-0.5*(1.5*t+3)).*u(1.5*t+3);
subplot(211)
plot(t,x)
title('x(t)')
subplot(212)
plot(t,y)
title('y(t)=x(1,5*t+3)')
1-3
clear;
close all;
t0=0;t1=3.5;dt=0.001;
t=t0:dt:t1;
x1=u(t)-u(t-2);
x2=u(t-1)-u(t-3);
y=dt*conv(x1,x2);
subplot(221)
plot(t,x1)
grid on
title('Signal x1(t)')
axis([t0,t1,0,1.5])
subplot(222)
plot(t,x2)
grid on
title('signal x2(t)')
axis([t0,t1,0,1.5])
subplot(212)
t=t0:dt:2*t1;
plot(t,y)
grid on
title('The convolution of x1(t) and x2(t)')
axis([2*t0,2*t1,0,2.5])
xlabel('Time t(sec)')
+
2-1
clear
close all
num=input('Type in the right coeffcient vector of differential equation:');
den=input('Type in the left coefficient vector of differential equation:');
t=0:0.01:10;
x=input('Type in the expression of the input signal x(t):');
lsim(num,den,x,t)
运行过程:
Type in the right coeffcient vector of differential equation:8
Type in the left coefficient vector of differential equation:[1,3,2]
Type in the expression of the input signal x(t):(exp(-2*t)-exp(-3*t)).*u(t)
3-1
clear
close all
T=2;
dt=0.00001;
t=-2:dt:2;
x1=(u(t+1)-u(t-1-dt)).*abs(t);
x=0;
for m=-1:1
x=x+u(t-m*T)-u(t-1-m*T-dt);
end
w0=2*pi/T;
N=10;
L=2*N+1;
for k=-N:N;
ak(N+1+k)=(1/T)*x1*exp(-j*k*w0*t')*dt;
end
phi=angle(ak);
运行结果:
ak =
Columns 1 through 6
0.0000 - 0.0000i -0.0025 - 0.0000i 0.0000 + 0.0000i -0.0041 + 0.0000i
0.0000 + 0.0000i -0.0081 - 0.0000i
Columns 7 through 12
0.0000 - 0.0000i -0.0225 + 0.0000i 0.0000 - 0.0000i -0.2026 + 0.0000i
0.5000 -0.2026 - 0.0000i
Columns 13 through 18
0.0000 + 0.0000i -0.0225 - 0.0000i 0.0000 + 0.0000i -0.0081 + 0.0000i
0.0000 - 0.0000i -0.0041 - 0.0000i
Columns 19 through 21
0.0000 - 0.0000i -0.0025 + 0.0000i 0.0000 + 0.0000i
3-2
b=input('b=');
a=input('a=');
[H,w]=freqs(b,a);
Hm=abs(H);
phai=angle(H);
hr=real(H);
Hi=i
文档评论(0)