matlab如何解这个差分方程
已知地表热通量G(T)，假设路面为均匀分布的均质体，根据热力学方程，t时刻单位质量的路面温度 T（t）和G(T)之间的关系为：
dT（t）/dt=dG(T)/C*dt
这里C 为常数,如果我已知t时刻的所有值,如何求得t+1时刻的路面温度呢?对时间进行差分吗?我想知道用matlab 或者其他的软件具体的做法,最好有代码!
正在做一个有关的论文，对这个很是迷惑,请大虾指点,不甚感激!!!
>> dsolve('DT=DG/C')
ans =
T(t)=1/C*G(t)+C1
用matlab怎么解差分方程?
syms z n u=(1\/5)^n; %设置u为输入信号 U=ztrans(u); %z变换 y=diff_eq([1 -5 6],[0],[1 0],U); %解析解。y=simplify(y)其运行结果，y(n)=3*2^n - 2*3^n 例2：试用matlab求解差分方程，y(t+2)+y(t+1)-6y(t)=5*2^t,y(0)=1，其中y(0)=1，y(1)=0...
怎么用matlab求解这道差分方程?
怎么用matlab求解这道差分方程？根据我的理解其基本思路是：1、根据差分方程式，写出迭代式，利用for循环语句和if判断语句来书写。即 for i=3:n if (k(i)-1)<0 uk1=0;else uk1=k(i)-1;end if (k(i)-2)<0 uk2=0;else uk2=k(i)-2;end y(i)=2*uk1-2*uk2+3*y(i-1)-...
怎么用matlab求差分
调用filter函数解差分方程。1）yn=filter(B,A,xn)是计算输入向量xn的零状态响应输出信号yn，yn与xn长度相等，其中B=[b0,b1，……bn], A=[a0,a1,……an]。其中a0=1。2）yn=filter(B,A,xn,xi)是计算全响应的函数。xi是等效初始条件的输入序列，xi能由初始条件确定。此时需要调用filtic函数。
在MATLAB中利用filter函数在给定输入和差分方程时求差分方程的解。给...
在MATLAB中，可以用函数y=filter(p,d,x)实现差分方程的仿真，也可以用函数 y=conv(x,h)计算卷积，用y=impz(p,d,N)求系统的冲激响应。实现差分方程 先从简单的说起：filter([1,2],1,[1,2,3,4,5])实现 y[k]=x[k]+2*x[k-1]y[1]=x[1]+2*0=1%(x[1]之前状态都用0)y[...
matlab如何解这个差分方程
这样可以吗？>> dsolve('DT=DG\/C')ans = T(t)=1\/C*G(t)+C1
matlab利用递归求解差分方程
首先，这个不是matlab利用递归求解差分方程，而是递推；差分方程其实就是递推关系式。然后这个循环：for i=N+1:N+length(n),y(i)= -a1*y(i-N:i-1)'+ b1*x(i-N:i-N+M)';end 其实是因为：y[n]+ a1*y[n-1]+ a2*y[n-2]...+ an*y[n-N]= b0*x[n]+ b1*x[n-1]+...
用MATLAB求解 已知系统的差分方程为y(n)=-a1y(n-1)-a2y(n-2)+bx(n...
这是一个二阶系统 a=[1,-0.8,0.64];% y(n)系数 b=[1,0,0];%x(n)系数 点序列 n=1:49;求单位脉冲响应序列 hn=impz(b,a,n);成图 stem(n,hn,'k','f');添加坐标及图形名称 ylabel('脉冲响应（h(n)');xlabel('序号（n）');title('单位脉冲响应');如果有问题问题，可以...
关于matlab filter这个函数到底怎么用?
一、在MATLAB中,可以用函数y=filter(p,d,x)实现差分方程的仿真,即y=filter(p,d,x)用来实现差分方程,d表示差分方程输出y的系数,p表示输入x的系数,而x表示输入序列。输出结果长度数等于x的长度。 举例:实现差分方程,先从简单的说起:filter([1,2],1,[1,2,3,4,5]),实现y[k]=x[k]+2*x[k-1]...
怎么用matlab解方程啊?有谁知道
MATLAB使用练习 l 掌握MATLAB的基本使用方法 实验1 数学建模初步 l 通过实例了解数学建模的一般步骤;l 在以后的数学实验中用数学建模方法解决经过简化的实际问题;l 自觉培养用数学方法解决实际问题的意识和能力。实验2 差分方程与数值微分 l 一阶常系数差分方程；l 高阶常系数差分方...
MATLAB 差分方程求解
这个问题要过细想，a1=20,a4=5,a2=a3=a5=0.1 s(t+1)=s(t)-a1*s(t)*i(t)-a2*a3*s(t)*i(t)\/(s(t)+e(t)+i(t));e(t+1)=e(t)+a1*s(t)*i(t)-e(t-a4)-a2*a3*(e(t)-i(t))\/(s(t)+e(t)+i(t));i(t+1)=i(t)+e(t-a4)-a2*i(t)^2\/(s(t)...
