扫扫二维码，随身浏览文档
手机或平板扫扫即可继续访问
《MATLAB程序设计教程(第二版)》第6章__MATLAB解方程与最优化问题求解
举报该文档为侵权文档。
举报该文档含有违规或不良信息。
反馈该文档无法正常浏览。
举报该文档为重复文档。
推荐理由：
将文档分享至：
分享完整地址
文档地址：
粘贴到BBS或博客
flash地址：
支持嵌入FLASH地址的网站使用
html代码：
&embed src='/DocinViewer-4.swf' width='100%' height='600' type=application/x-shockwave-flash ALLOWFULLSCREEN='true' ALLOWSCRIPTACCESS='always'&&/embed&
450px*300px480px*400px650px*490px
支持嵌入HTML代码的网站使用
您的内容已经提交成功
您所提交的内容需要审核后才能发布，请您等待！
3秒自动关闭窗口温馨提示！由于新浪微博认证机制调整，您的新浪微博帐号绑定已过期，请重新绑定！&&|&&
LOFTER精选
阅读(1395)|
用微信&&“扫一扫”
将文章分享到朋友圈。
用易信&&“扫一扫”
将文章分享到朋友圈。
历史上的今天
loftPermalink:'',
id:'fks_',
blogTitle:'matlab中采用ode45解方程的问题',
blogAbstract:'&&&&&& 帮师兄调程序来的，就一个ode45中参数传递的问题，捣鼓了好一会儿。
&&&&&& 现将ode45参数之间的对应关系示例如下：
&&&&&& function dy=rigid(t,y,vv)
&&&&&& %%%%%%%%%%格式为：function 输出=函数名字（输入）
&&&&&& A=1;
&&&&&& afa=1;
&&&&&& gama=0.5;
&&&&&& n=0.5;
&&&&&& dy=A*vv-(afa+gama)*vv*y^n;
&&&&&& 在main程序中调用情况如下：
&&&&&& options=。。。。。
blogTag:'',
blogUrl:'blog/static/',
isPublished:1,
istop:false,
modifyTime:2,
publishTime:7,
permalink:'blog/static/',
commentCount:0,
mainCommentCount:0,
recommendCount:1,
bsrk:-100,
publisherId:0,
recomBlogHome:false,
currentRecomBlog:false,
attachmentsFileIds:[],
groupInfo:{},
friendstatus:'none',
followstatus:'unFollow',
pubSucc:'',
visitorProvince:'',
visitorCity:'',
visitorNewUser:false,
postAddInfo:{},
mset:'000',
remindgoodnightblog:false,
isBlackVisitor:false,
isShowYodaoAd:false,
hostIntro:'',
hmcon:'0',
selfRecomBlogCount:'0',
lofter_single:''
{list a as x}
{if x.moveFrom=='wap'}
{elseif x.moveFrom=='iphone'}
{elseif x.moveFrom=='android'}
{elseif x.moveFrom=='mobile'}
${a.selfIntro|escape}{if great260}${suplement}{/if}
{list a as x}
推荐过这篇日志的人：
{list a as x}
{if !!b&&b.length>0}
他们还推荐了：
{list b as y}
转载记录：
{list d as x}
{list a as x}
{list a as x}
{list a as x}
{list a as x}
{if x_index>4}{break}{/if}
${fn2(x.publishTime,'yyyy-MM-dd HH:mm:ss')}
{list a as x}
{if !!(blogDetail.preBlogPermalink)}
{if !!(blogDetail.nextBlogPermalink)}
{list a as x}
{if defined('newslist')&&newslist.length>0}
{list newslist as x}
{if x_index>7}{break}{/if}
{list a as x}
{var first_option =}
{list x.voteDetailList as voteToOption}
{if voteToOption==1}
{if first_option==false},{/if}&&“${b[voteToOption_index]}”&&
{if (x.role!="-1") },“我是${c[x.role]}”&&{/if}
&&&&&&&&${fn1(x.voteTime)}
{if x.userName==''}{/if}
网易公司版权所有&&
{list x.l as y}
{if defined('wl')}
{list wl as x}{/list}matlab用龙格库塔法求解微分方程组题目是这样的：将微分方程组重写为一阶以后是这样的初始矢量不知道在matlab中应该怎么写,分数不多,但是是我有的全部了 TAT_百度作业帮
matlab用龙格库塔法求解微分方程组题目是这样的：将微分方程组重写为一阶以后是这样的初始矢量不知道在matlab中应该怎么写,分数不多,但是是我有的全部了 TAT
解方程很简单,调用ode45即可关键是求运动时间,假设z方向坐标变为0即为落至地面运动结束,可求得运动时间的近似值及落地点坐标代码如下：clear allclcc=0.07;u0=25;a=pi/6;f=@(t,x)([x(2);-c*sqrt(x(2)^2+(x(4)-7-.35*x(5))^2+x(6)^2)*x(2);x(4);-c*sqrt(x(2)^2+(x(4)-7-.35*x(5))^2+x(6)^2)*(x(4)-7-.35*x(5));x(6);-9.81-c*sqrt(x(2)^2+(x(4)-7-.35*x(5))^2+x(6)^2)*x(6)]);[t,Y]=ode45(f,0:.01:2,[0 u0*cos(a) 0 0 1.4 u0*sin(a)]);n=max(find(Y(:,5)>0));t_end=t(n)s=[Y(n,1),Y(n,3),0]plot3(Y(1:n,1),Y(1:n,3),Y(1:n,5))其中t_end为运动时间s为落地点坐标看了半天,怎么都觉得方程不太对再检查一下方程吧,尤其是第二个方程等号右边到底有没有负号最后一个方程等号右边也有问题,请仔细检查matlab solve global bglobal nglobal yt1=0.1;t2=0.001;t3=0.005;t4=0.01;t5=0.05;t6=0.1;t7=0.002;t8=0.006;t9=0.07;t10=0.1;shizi1='b/n*((t1-y)^(b-1)+(t2-y)^(b-1)+(t3-y)^(b-1)+(t4-y)^(b-1)+(t5-y)^(b-1)+(t6-y)^(b-1)+(t7-y)^(b-1)+(t8-y)^(b-1)+(t9-y)^(b-1)+(t10-y)^(b-1))-(b-1)*((t1-y)^(-_百度作业帮
matlab solve global bglobal nglobal yt1=0.1;t2=0.001;t3=0.005;t4=0.01;t5=0.05;t6=0.1;t7=0.002;t8=0.006;t9=0.07;t10=0.1;shizi1='b/n*((t1-y)^(b-1)+(t2-y)^(b-1)+(t3-y)^(b-1)+(t4-y)^(b-1)+(t5-y)^(b-1)+(t6-y)^(b-1)+(t7-y)^(b-1)+(t8-y)^(b-1)+(t9-y)^(b-1)+(t10-y)^(b-1))-(b-1)*((t1-y)^(-1)+(t2-y)^(-1)+(t3-y)^(-1)+(t4-y)^(-1)+(t5-y)^(-1)+(t6-y)^(-1)+(t7-y)^(-1)+(t8-y)^(-1)+(t9-y)^(-1)+(t10-y)^(-1))=0';shizi2='-j/n+1/n^2*((t1-y)^b+(t2-y)^b+(t3-y)^b+(t4-y)^b+(t5-y)^b+(t6-y)^b+(t7-y)^b+(t8-y)^b+(t9-y)^b+(t10-y)^b)=0';shizi3='n/b+(log(t1-y)+log(t2-y)+log(t3-y)+log(t4-y)+log(t5-y)+log(t6-y)+log(t7-y)+log(t8-y)+log(t9-y)+log(t10-y))-1/n*((t1-y)^b*log(t1-y)+(t2-y)^b*log(t2-y)+(t3-y)^b*log(t3-y)+(t4-y)^b*log(t4-y)+(t5-y)^b*log(t5-y)+(t6-y)^b*log(t6-y)+(t7-y)^b*log(t7-y)+(t8-y)^b*log(t8-y)+(t9-y)^b*log(t9-y)+(t10-y)^b*log(t10-y))=0';[b,y,n]=solve(shizi1,shizi2,shizi3,b,y,n);哪儿错了?错误提示：Error using mupadengine/feval (line 157)MuPAD error:Error:']' expected.[line 1,col 35]Error in solve (line 160)sol = eng.feval('symobj::solvefull',eqns,vars);Error in disanti (line 17)[a,d,f]=solve(shizi1,shizi2,shizi3,b,y,
将前3行全局变量定义改为符号变量定义：syms b n y.希望你的shizi表达式有意义.我的机子求解你的问题一直busy
没有报错啊，结果是：ans =
-((COH0^2*COHadd^2*Kb^4 + 2*COH0^2*COHadd*Kb^3 + 4*COH0^2*Cblg*Kb^3*N1 + 4*COH0^2*Cblg*Kb^3*N2 + 4*COH0^2*Cblg*Kb^3*N3 + 4*COH0^2*Cblg*Kb^3*N4 + 4*COH0^2*Cblg*Kb^3*N5 +...matlab解方程解得问题：如何知道得到的解对应哪个未知量_百度作业帮
matlab解方程解得问题：如何知道得到的解对应哪个未知量
不要用[x,y]=这个和方程中的变量x,y没有任何联系,不是传址调用的标准用法是a=solve('x+y=2','x-y=1')然后你会看到a = x:[1x1 sym]y:[1x1 sym]然后输入a.x就是你的xa.y是你的y这样清楚了吧不然你直接定义向量你不知道哪个是哪个的一般是按照字符顺序x,y,z来排的所以一般你知道第一个出来的是x的值但是你赋给y没有任何问题,就是为什么你改[x,y],[y,x]第一个值都是3/2的原因
你help solve就知道了它的用法了，