24小时热门版块排行榜

返回列表

lhdcyx

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 44.5
帖子: 16
在线: 16.3小时
虫号: 3640354
注册: 2015-01-12
专业: 动力学与控制

[求助] 帮忙看个简单的matlab程序已有1人参与

一段求两点边值问题的matlab程序：
clc
clear all
f=@(t,x,v)[4*x(1)+v(1)*x(1)*x*(2);-2*x(2)+v(2)*x(1)*x(2)];
g=@(ya,yb,v)[ya(1)-2;ya(2)-1;yb(1)-4;yb(2)-2];
x1=[1;1];x2=[-1;1];
solinit=bvpinit(linspace(0,3,20),x1,x2);
sol=bvp5c(f,g,solinit);
sol.parameters
plot(sol.x,sol.y);
figure;
plot(sol.y(1,

,sol.y(2,

);
报错，出不来结果。不知道怎么回事？

回复此楼

» 猜你喜欢

表哥与省会女结婚，父母去帮带孩子被省会女气回家生重病了已经有12人回复
依托企业入选了国家启明计划青年人才。有无高校可以引进的。已经有14人回复
江汉大学解明教授课题组招博士研究生/博士后已经有3人回复
AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性已经有11人回复
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗已经有10人回复
依托企业入选了国家启明计划青年人才。有无高校可以引进的。已经有11人回复

1楼 2016-09-20 16:39:42

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)

应助: 591 (博士)
贵宾: 0.03
金币: 18247.2
红花: 97
帖子: 8891
在线: 1487.9小时
虫号: 2400059
注册: 2013-04-04
专业: 计算机应用技术

【答案】应助回帖

感谢参与，应助指数 +1

CODE:

>> f

f = 

    @(t,x,v)[4*x(1)+v(1)*x(1)*x*(2);-2*x(2)+v(2)*x(1)*x(2)]

>> g

g = 

    @(ya,yb,v)[ya(1)-2;ya(2)-1;yb(1)-4;yb(2)-2]

>> sol=bvp5c(f,g,solinit);

??? Error using ==> bvparguments at 121

Error in calling BVP5C(ODEFUN,BCFUN,SOLINIT):

  The derivative function ODEFUN should return a column

  vector of length 2.

Error in ==> bvp5c at 137

[neqn,nparam,nregions,atol,rtol,Nmax,xyVectorized,printstats]

= ...

>> whos 

  Name         Size            Bytes  Class              Attributes

  f            1x1                16  function_handle              

  g            1x1                16  function_handle              

  solinit      1x1              1146  struct                       

  x1           2x1                16  double                       

  x2           2x1                16  double

https://www.mathworks.com/matlabcentral/fileexchange/3819-tutorial-on-solving-bvps-with-bvp4c