matlab微分方程求解中如何对方程进行编程求解啊

来源：蜘蛛抓取(WebSpider) 时间：2015-04-12 13:09 标签： matlab偏微分方程求解
•回答
•回答
•回答
•回答
•回答
终值xf%y0——初始值%h——步长%t——节点向量%y——数值解x1=0,ylabel(');,y1,'.01) %[t,y]=Euler(odefun,[0,xf)初值x0;3);y&#39,&#39,y]=Euler(plot(t,6],tspan,legend(&#39,x&#39.5,y)%tspan——(x0clear all.^1;),0;clcodefun=inline('y'y1=((x1;[t;k*-&#39,&#39,');x*y^(1&#47,1;r-'3)'x&#39.^2+2)/);精确解&#39,x1,&#39,xlabel(&#39,y0.1,h)%odefun——微分方程函数f(x:0,y;:6;解曲线'
猜你感兴趣
服务声明： 信息来源于互联网，不保证内容的可靠性、真实性及准确性，仅供参考，版权归原作者所有！Copyright &
.&&闽ICP备号&&&
Powered by视频: 【matlab数学建模代做】基于matlab编程解线性方程源码程序
分享给好友
您需要先安装&，才能下载视频哦
用优酷App或微信扫一扫，在手机上继续观看。
【matlab数学建模代做】基于matlab编程解线性方程源码程序
分享给站外好友
把视频贴到Blog或BBS
flash地址:
<input type="text" class="form_input form_input_s" id="link3" value=''>
<input id="link4" type="text" class="form_input form_input_s" value=''>
基于matlab编程解线性方程源码程序
公司主营：个人项目、学校项目、公司项目 matlab代做 matlab代写 matlab有偿编程【谷速软件】欢迎你
节目制作经营许可证京字670号
京公网安备号
药品服务许可证(京)-经营-如何在matlab下求解三角函数方程？
如何在matlab下求解三角函数方程？
09-12-24 & 发布
%%下面方法比较笨，但是可行。%%先看一下为什么你解不出来%%把你的k*1.57e-6看成一个整体x%%变成求解tan(x)./x-1/0.2159==0x=-2:0.0001:2;y=tan(x)./x-1/0.2159;plot(x,y,x,zeros(size(x)))title(&#39;[-2,2]&#39;)figurex=-1.5:0.;y=tan(x)./x-1/0.2159;plot(x,y,x,zeros(size(x)))title(&#39;[-1.5,1.5]&#39;)%%从上面两张图可以看出%%第一幅图中由于存在-1.57(-pi/2)时，tan变成无穷大%%这里精确解我也不知道怎么取到就是了。%%不过放大第二幅图(放大那个交点)，解应当在-1.420到-1.419之间figurex=-1.420:0..419;y=tan(x)./x-1/0.2159;plot(x,y,x,zeros(size(x)))for i=1:length(x)-1 if (y(i)*y(i+1)&0) &#160; endendvpa(x(i),10)得到结果为-1.对称性还可以得到1.你可以逐步求精，上面只是取到小数点后7位，你还可以再往下做，效率不是很高就是了。当然上面得到的是x你的k应当是k=x(i)/1.57e-6;vpa(k,10)tan(k*(1.57E-6))./(k*(1.57E-6))-1/0.2159另外-pi/2以外的区域，还有根，你可以类同地处理。
请登录后再发表评论!
%%下面方法比较笨，但是可行。%%先看一下为什么你解不出来%%把你的k*1.57e-6看成一个整体x%%变成求解tan(x)./x-1/0.2159==0x=-2:0.0001:2;y=tan(x)./x-1/0.2159;plot(x,y,x,zeros(size(x)))title(&#39;[-2,2]&#39;)figurex=-1.5:0.;y=tan(x)./x-1/0.2159;plot(x,y,x,zeros(size(x)))title(&#39;[-1.5,1.5]&#39;)%%从上面两张图可以看出%%第一幅图中由于存在-1.57(-pi/2)时，tan变成无穷大%%这里精确解我也不知道怎么取到就是了。%%不过放大第二幅图(放大那个交点)，解应当在-1.420到-1.419之间figurex=-1.420:0..419;y=tan(x)./x-1/0.2159;plot(x,y,x,zeros(size(x)))for i=1:length(x)-1 if (y(i)*y(i+1)&0) &#160; endendvpa(x(i),10)得到结果为-1.对称性还可以得到1.你可以逐步求精，上面只是取到小数点后7位，你还可以再往下做，效率不是很高就是了。当然上面得到的是x你的k应当是k=x(i)/1.57e-6;vpa(k,10)tan(k*(1.57E-6))./(k*(1.57E-6))-1/0.2159另外-pi/2以外的区域，还有根，你可以类同地处理。
请登录后再发表评论!用matlab求解下面方程，请编写程序，用m.文件
用matlab求解下面方程，请编写程序，用m.文件
这个问题我知道!functiony=zhongy(x)%保存为zhongy.m文件y=[2*x(1)-x(2)-exp(-x(1));-x(1)+2*x(2)-exp(-x(2))];------------------------------------------------------------[x,fval,exitflag]=fsolve('zhongy',[0.1;0.1])%x为所求的解%fval为解代入方程的验证结果%exitflag=1表示函数收敛于解；反之，=0
相关知识等待您来回答
编程领域专家

我要回帖

更多关于 matlab偏微分方程求解 的文章

 

随机推荐