题目
matlab求解泛函微分方程数值解.
我想用matlab求解一个泛函微分方程的数值解,编了个程序总是运行不出来,
lags=[1];
sol=dde23('ddefun',lags,'history',[0,100]);
plot(sol.x,sol.y);
title('ddefun');xlabel('t');ylabel('y');
legend('y_1','y_2','y_3','y_4',2);
tx=[4.3,4.7];y=deval(sol,tx);hold on;plot(tx,y,'x');
function dydt=ddefun(t,y,Z)
ylag1=Z(:,1);
dydt=[(-0.5)*y(1)-2*y(2)+ylag1(3)
2*y(1)-0.5*y(2)+ylag1(4)
(-0.5)*y(3)-0.5*y(4)+ylag1(1)
0.5*y(3)-0.5*y(4)+ylag1(2)];
function S=history(t)
S=ones(3,1);
原方程是这样的
dy1/dt=-0.5y1-2y2+y3(t-1)
dy2/dt=2y1-0.5y2+y4(t-1)
dy3/dt=-0.5y3-0.5y4+y1(t-1)
dy4/dt=0.5y3-0.5y4+y2(t-1)
还是不对,运行后出现了这样的提示:
In an assignment A(:,matrix) = B,the number of rows in A and B
must be the same.
Error in ==> C:\MATLAB6.5\toolbox\matlab\funfun\dde23.m (lagvals)
On line 656 ==> Z(:,j) = temp(:);
Error in ==> C:\MATLAB6.5\toolbox\matlab\funfun\dde23.m
On line 270 ==> Z = lagvals(t,lags,history,t,y,[],varargin{:});
我想用matlab求解一个泛函微分方程的数值解,编了个程序总是运行不出来,
lags=[1];
sol=dde23('ddefun',lags,'history',[0,100]);
plot(sol.x,sol.y);
title('ddefun');xlabel('t');ylabel('y');
legend('y_1','y_2','y_3','y_4',2);
tx=[4.3,4.7];y=deval(sol,tx);hold on;plot(tx,y,'x');
function dydt=ddefun(t,y,Z)
ylag1=Z(:,1);
dydt=[(-0.5)*y(1)-2*y(2)+ylag1(3)
2*y(1)-0.5*y(2)+ylag1(4)
(-0.5)*y(3)-0.5*y(4)+ylag1(1)
0.5*y(3)-0.5*y(4)+ylag1(2)];
function S=history(t)
S=ones(3,1);
原方程是这样的
dy1/dt=-0.5y1-2y2+y3(t-1)
dy2/dt=2y1-0.5y2+y4(t-1)
dy3/dt=-0.5y3-0.5y4+y1(t-1)
dy4/dt=0.5y3-0.5y4+y2(t-1)
还是不对,运行后出现了这样的提示:
In an assignment A(:,matrix) = B,the number of rows in A and B
must be the same.
Error in ==> C:\MATLAB6.5\toolbox\matlab\funfun\dde23.m (lagvals)
On line 656 ==> Z(:,j) = temp(:);
Error in ==> C:\MATLAB6.5\toolbox\matlab\funfun\dde23.m
On line 270 ==> Z = lagvals(t,lags,history,t,y,[],varargin{:});
提问时间:2021-01-02
答案
改动:
function S=history(t)
S=ones(3,1);
改为 S=ones(4,1);
OVER!
PS.多翻翻help!
function S=history(t)
S=ones(3,1);
改为 S=ones(4,1);
OVER!
PS.多翻翻help!
举一反三
已知函数f(x)=x,g(x)=alnx,a∈R.若曲线y=f(x)与曲线y=g(x)相交,且在交点处有相同的切线,求a的值和该切线方程.
我想写一篇关于奥巴马的演讲的文章,写哪一篇好呢?为什么好
奥巴马演讲不用看稿子.为什么中国领导演讲要看?
想找英语初三上学期的首字母填空练习……
英语翻译
最新试题
- 1如图:计算下面各个图形的表面积和体积
- 2一个正方形被分成了3个长方形,每个长方形的周长都是80厘米,这个正方形的周长
- 3He has got much milk.这句话该怎么翻译,又是什么时态?
- 4夫子何命焉为的“焉为”是什么意思
- 5微风吹拂着孩子们的脸改成拟人句
- 61分钟的跳绳比赛,龙龙跳了120下亮亮跳的是龙龙跳的8分之5丽丽跳的是亮亮跳的5分之3谁得冠军
- 7已知a,b、c三数在数轴上的位置如图所示,化简|a|a+|b|b+|c|c.
- 8现有1克,2克,4克和8克的砝码各一个,在天平上可以称出多少种不同重量的物体?要一一列举
- 9力透纸背,入木三分,龙飞凤舞这些成语最早都是形容什么艺术的
- 10计算面积怎么算
热门考点
- 1有两组数,第一组平均数为12.6,第二组平均数为10.8,而两组总的平均数是11.5,两组的个数之和是几(写明过程)
- 240比()多百分之十
- 3在夏季,我们看得到北极星吗
- 4咖啡属于酸性还是碱性
- 5行星跟恒星的形成有怎样的区别?为什么它们的组成成分差异那么大?
- 6一mol硝酸中含有多少个氢原子,多少个氧原子?
- 7一个长方体长、宽、高各缩小到原来的二分之一,表面积缩小到原来的【 】,体积缩小到原来的【 】.
- 8一种节能型冰箱,商店按成本价加价20%作为标价,一顾客与销售员讨价还价后,以标价的
- 9在平面直角坐标系中,已知点A(3,0),P是圆x^2-y^2=1上一个动点,且角AOP的平分线交PA
- 10最后一课小弗朗士是从哪些方面得出整个教室有一种不平常的严肃的气愤的