君，已阅读到文档的结尾了呢~~
扫扫二维码，随身浏览文档
手机或平板扫扫即可继续访问
非线性方程几种数值解法的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秒自动关闭窗口查看: 8324|回复: 15|关注: 0
【求助】怎样用matlab解这样一个非线性方程
<h1 style="color:# 麦片财富积分
新手, 积分 5, 距离下一级还需 45 积分
各位热心的朋友，请大家帮帮忙，不胜感激！
方程的具体表达式：sqrt(0.067*(0.127-y)/(0.059*y))-tan(sqrt(0.059*9.1*10^(-31)*y*1.6*10^(-19)*(x*10^(-9))^2/2/(6.63*10^(-34)/(2*3.14))^2))=0
为了大家看起来方便，上面的常数用符号来表示，可将方程写成这种形式：
& && &sqrt(A*(B-y)/y)-tan(sqrt(C*y*x^2))=0
其中A、B、C是已知常数，x，y为未知数。
需要解决的问题是给定x的一个值（例如x=8），怎样求出y在0——0.13之间的实数解？希望大家能给出程序！
我用solve命令计算过，结果是复数，而且不对。
[ 本帖最后由 zhls1209 于
20:54 编辑 ]
<h1 style="color:#0 麦片财富积分
关注者： 9
eq1=sqrt(0.067*(0.127-y)/(0.059*y))-tan(sqrt(0.059*9.1*10^(-31)*y*1.6*10^(-19)*(x*10^(-9))^2/2/(6.63*10^(-34)/(2*3.14))^2));
y=vpa(solve(eq1),4)
1.022+.5378*i
<h1 style="color:# 麦片财富积分
回复 2# ww711 的帖子
不好意思，我没说清楚，我只要实数解，不要复数的，而且是在0——0.13之间的解，最好是不用初值的，如果实在没办法，带初值的也可以。我代进去验算了一下好像式子不为0啊。
[ 本帖最后由 zhls1209 于
20:55 编辑 ]
<h1 style="color:#4 麦片财富积分
关注者： 6
x=8时，有许多根吧？
& &4.779e-002& && & 0.6102& && & 0.6107& && & 0.6107& && & 0.6106& && & 0.6106& && & 0.6106& && & 0.6106& && & 0.6106& && & 0.6106& && & 0.6106& && & 0.6105& && & 0.6105& && & 0.6105& && & 0.6105& && & 0.6105& && & 0.6105& && & 0.6104& && & 0.6104& && & 0.6104& && & 0.6104& && & 0.6104& && & 0.6104
还可能更多。
<h1 style="color:#4 麦片财富积分
关注者： 6
4.94e-002& && & 0.8963& && & 0.8886& && & 0.8886& && & 0.8886& && & 0.8886& && & 0.8887& && & 0.8887& && & 0.8887& && & 0.8887& && & 0.8887& && & 0.8887& && & 0.8888& && & 0.8887& && & 0.8888& && & 0.8888& && & 0.8888& && & 0.8888& && & 0.8888& && & 0.8888& && & 0.8889& && & 0.8889& && & 0.8889& && & 0.8889& && & 0.8889& && & 0.8889& && &&&0.889& && &&&0.889& && &&&0.889& && &&&0.889& && &&&0.889& && &&&0.889& && & 0.8891& && &&&0.889& && & 0.8891& && & 0.8891& && & 0.8891& && & 0.8891& && & 0.8891& && & 0.8891& && & 0.8892& && & 0.8892& && & 0.8892& && & 0.8892& && & 0.8892& && & 0.8892& && & 0.8893& && & 0.8893& && & 0.8893& && & 0.8893& && & 0.8893& && & 0.8893& && & 0.8894& && & 0.8893& && & 0.8894& && & 0.8894& && & 0.8894& && & 0.8894& && & 0.8894& && & 0.8894& && & 0.8895& && & 0.8895& && & 0.8895& && & 0.8895& && & 0.8895& && & 0.8895& && & 0.8896& && & 0.8896& && & 0.8896& && & 0.8896& && & 0.8896& && & 0.8896& && & 0.8897& && & 0.8896& && & 0.8897& && & 0.8897& && & 0.8897& && & 0.8897& && & 0.8897& && & 0.8897& && & 0.8898& && & 0.8898& && & 0.8898& && & 0.8898& && & 0.8898& && & 0.8898& && & 0.8899& && & 0.8899& && & 0.8899& && & 0.8899& && & 0.8899& && & 0.8899& && && &0.89& && & 0.8899& && && &0.89& && && &0.89& && && &0.89& && && &0.89& && && &0.89& && && &0.89& && & 0.8901& && & 0.8901& && & 0.8901& && & 0.8901& && & 0.8901& && & 0.8901& && & 0.8902& && & 0.8902& && & 0.8902& && & 0.8902& && & 0.8902& && & 0.8902& && & 0.8903& && & 0.8902& && & 0.8903& && & 0.8903& && & 0.8903& && & 0.8903& && & 0.8903& && & 0.8903& && & 0.8904& && & 0.8904& && & 0.8904& && & 0.8904& && & 0.8904& && & 0.8904& && & 0.8905& && & 0.8905& && & 0.8905& && & 0.8905& && & 0.8905& && & 0.8905& && & 0.8906& && & 0.8905& && & 0.8906& && & 0.8906& && & 0.8906& && & 0.8906& && & 0.8906& && & 0.8906& && & 0.8907& && & 0.8907& && & 0.8907& && & 0.8907& && & 0.8907& && & 0.8907& && & 0.8908& && & 0.8908& && & 0.8908& && & 0.8908& && & 0.8908& && & 0.8908& && & 0.8909& && & 0.8908& && & 0.8909& && & 0.8909& && & 0.8909& && & 0.8909& && & 0.8909& && & 0.8909& && &&&0.891& && &&&0.891& && &&&0.891& && &&&0.891& && &&&0.891& && &&&0.891& && & 0.8911& && & 0.8911& && & 0.8911& && & 0.8911& && & 0.8911& && & 0.8911& && & 0.8912& && & 0.8911& && & 0.8912& && & 0.8912& && & 0.8912& && & 0.8912& && & 0.8912& && & 0.8912& && & 0.8913& && & 0.8913& && & 0.8913& && & 0.8913& && & 0.8913& && & 0.8913& && & 0.8914& && & 0.8914& && & 0.8914& && & 0.8914& && & 0.8914& && & 0.8914& && & 0.8915& && & 0.8914& && & 0.8915& && & 0.8915& && & 0.8915& && & 0.8915& && & 0.8915& && & 0.8915& && & 0.8916& && & 0.8916& && & 0.8916
<h1 style="color:#4 麦片财富积分
关注者： 6
不好意思，好多重复的。
<h1 style="color:# 麦片财富积分
给一个x值都能得到0-0.13的y值？
感觉问题有点怪？
<h1 style="color:# 麦片财富积分
回复 6# forcal 的帖子
谢谢你帮我，麻烦你能不能把程序发上来，我不知道你是怎么算的啊！谢谢啊！
<h1 style="color:#4 麦片财富积分
关注者： 6
验证结果有些问题，我再看一下。
<h1 style="color:#4 麦片财富积分
关注者： 6
似乎只有第一个是对的。应该先绘图再求解，又犯了一次错误。
站长推荐 /2
快速搭建新能源汽车整车模型及其性能优化
MATLAB中文论坛是全球最大的 MATLAB & Simulink 中文社区。用户免费注册会员后，即可下载代码，讨论问题，请教资深用户及结识书籍作者。立即注册加入我们吧！
MATLAB官方社交平台
MATLAB中文论坛微社区后使用快捷导航没有帐号？
中国汽车工程师之家--聚集了汽车行业80%专业人士&
论坛口号：知无不言，言无不尽！QQ:&
本站手机访问：直接在浏览器中输入本站域名即可&
本站广告咨询
微信:cartech8_cq
电话：023-
含变参数的单变量非线性方程求解-MATLAB
该用户从未签到
发表于 29-5-:56
& & & && &
& &&& & &&
& & && & & &
请教各位朋友，我有个非线性问题纠结了很久，请各位帮忙，问题是这样的：我想要求解一个含变参数的单变量非线性方程-MATLAB，我的调用函数是
global a cmod D GI s1 sa1 s2 sa2 fa
cmod=0.1214;D=600;GI=66.2454;
s1=int(sqrt((1-x/a)^2+(1.081-1.149*a/D)*(x/a-(x/a)^2)),x,240,a);s2=int((1-x/a)^2+(1.081-1.149*a/D)*(x/a-(x/a)^2),x,240,a);
sa1=simple(s1);sa2=simple(s2);
a=fzero('fz',335)
M文件中的定义函数为
function fa=fz(a)
&&syms a cmod D GI
&&syms s1 sa1 s2 sa2 fa
&&cmod=0.1214;GI=66.2454;
&&fa=1000/(a-240)*(2.79*cmod*sa1-19.4811*cmod^2*sa2)-GI;
单变量为a，可是求解总是提示错误如下:
??? Function 'isfinite' is not defined for values of class 'sym'.
Error in ==& fzero at 309
& & elseif ~isfinite(fx) || ~isreal(fx)
经查看，309行显示如下：
elseif ~isfinite(fx) || ~isreal(fx)
& && &&&error('MATLAB:fzero:ValueAtInitGuessComplexOrNotFinite',...
& && && && &'Function value at starting guess must be finite and real.');
请问是怎么回事，该怎样改动，才能求得a的值？
& & & & & & & & & &&
扫描订阅汽车工程师之家 & & & & & & & & & 汽车行业通信录-设计制造
& & & &微信公共账号 & & & & & & & & & & & & & & & & & & &零部件采购销售
上汽维修群
卡车维修群
该用户从未签到
发表于 29-5-:23
晕，俺也纠结，飘过~~
发动机QQ群：
电动车QQ群：
底盘QQ群：
工作时间：9:30-17:00
客服热线：023-
官方微信扫一扫
站务中心QQ:
GMT+8, 16-9-
Powered by&#xe621; 上传我的文档
&#xe602; 下载
&#xe60c; 收藏
该文档贡献者很忙，什么也没留下。
&#xe602; 下载此文档
用matlab求解非线性方程组的几种方法之程序
下载积分：1200
内容提示：用matlab求解非线性方程组的几种方法之程序
文档格式：PDF|
浏览次数：396|
上传日期： 10:54:07|
文档星级：&#xe60b;&#xe60b;&#xe60b;&#xe60b;&#xe60b;&#xe60b;
全文阅读已结束，如果下载本文需要使用
&#xe71b; 1200 积分
&#xe602;下载此文档
该用户还上传了这些文档
用matlab求解非线性方程组的几种方法之程序
官方公共微信