- 1、本文档共6页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
系統的一般分析方法
信号的表示与卷积实验报告
实验名称 系统的一般分析方法 实验时间 2013年 06月 18日 姓名 班级 成绩
一、实验目的
1.掌握差分方程的迭代求解,通过将微分化为差分的方法,把微分方程变成差。
2.掌握符号法求解微分方程解析解的方法。
3.掌握连续与离散系统冲激、阶跃响应求解的常用函数。
实验内容
1.微分和差分方程的求解
2.冲激响应(单位样值响应)和阶跃响应
3.系统分析
三、相关背景知识介绍
自然科学与社会科学中,很多动态变化的系统可以采用微分或差分方程建模。一些饶有兴味的建模案例将在第三部分进行研究。本节研究数学形式的微分和差分方程的计算机求解方法。
差分方程适合迭代求解,而求解微分方程的一般思路是将微分方程转化为差分方程。MATLAB提供的求解微分方程的常用指令有用于符号方程求解的dsolve和用于书追结算的ode指令集。
四、
1.解方程
(1),y(-1)=1,y(-2)=2, f(k)=ε(k),
求系统的h(k)、g(k)、y(k)。
输入: n=3:13;
F(n)=1;
F(1)=0;F(2)=0;
for k=3:13
Y(1)=2;
Y(2)=1;
for k=3:13
Y(k)=F(k)+F(k-1)-2*Y(k-1)-2*Y(k-2);
end
m=n-3;
stem(m,Y(m+1));xlabel(k);ylabel(y(k));grid on;
F(n)=0;
F(3)=1;
for k=3:13
Y(k)=F(k)+F(k-1)-2*Y(k-1)-2*Y(k-2);
end
stem(m,Y(m+1));xlabel(k);ylabel(h(k));grid on;
(2) (a)求解微分方程,初始条件为:y(0+)=y(0+)=1,。理论计算全响应,用ezplot绘结果图。范围为区间[0,4]。
(b)用数值法求解需先化为差分方程。化出对应的差分方程。步长取T=0.1。
(c)确定此差分方程的初始条件。
输入:
Ts=0.1;
n=1:4/Ts+3;
F(n)=1;F(1)=0;F(2)=0;
Y(1)=0.9;Y(2)=2;
for k=3:length(n)
Y(k)=(13*F(k)-10*F(k-1)-100*Y(k-2)+240*Y(k-1))/143;
end
m=(n-3)/Ts;
subplot(2,1,1);plot(m,Y);
xlabel(t);ylabel(y(t));title(全响应(数值计算));grid on;...
subplot(2,1,2);ezplot(-exp(-3*t)+exp(-t)+1,[0,4]);...
xlabel(t);ylabel(y(t));title(全响应);grid on;
电压源激励的RLC串联回路微分方程为:。求当R=1Ω,C=1F,L=1H时二阶电路的冲激响应和阶跃响应并绘波形图。
输入:
Ts=0.1;
n=1:4/Ts+3;
F(n)=0;
F(3)=1;
Y(1)=0;Y(2)=0;
for k=3:length(n)
Y(k)=[Ts^2*F(k)+(Ts+2)*Y(k-1)-Y(k-2)]/(Ts^2+Ts+1);
end
m=(n-3)/Ts;
plot(m,Y); %冲击响应
F(n)=1;
F(1)=0;F(2)=0;
Y(1)=0;Y(2)=0;
for k=3:length(n)
Y(k)=[Ts^2*F(k)+(Ts+2)*Y(k-1)-Y(k-2)]/(Ts^2+Ts+1);
end
m=(n-3)/Ts;
plot(m,Y); %阶跃响应
求解微分方程y+3y=3f(t),初始条件y(0)=1.5,y(0)=0,f(t)=ε(t)。
输入:
Ts=0.1;
n=1:4/Ts+3;
F(n)=1;F(1)=0;F(2)=0;
您可能关注的文档
最近下载
- 2024风力发电装置技术规范.docx VIP
- 歐姆龍AOI_OMRON自動光學檢測機編程培訓手冊(完美精華版).pdf VIP
- 人教精通版6年级上下册重点单词和句型默写.pdf VIP
- 社区获得性肺炎诊疗指南PPT课件.pptx VIP
- 钢筋笼吊装验算书.docx VIP
- 小学语文教学中人工智能辅助作文批改与评价策略研究教学研究课题报告.docx
- 美国仪表协会标准 ISA TR 88.95.01-2008 一起使用ISA-88和ISA-95.pdf VIP
- DB32T 4797-2024 旅游警务服务规范.docx VIP
- 《固定式钢梯及平台安全要求 第2部分:钢斜梯》GB 4053.2-2009.docx VIP
- 影视传播概论4.ppt VIP
文档评论(0)