系統的一般分析方法.docVIP

  1. 1、本文档共6页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 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;

文档评论(0)

df9v4fzI + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档