Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 工程数学软件 > MATLAB论坛
MATLAB论坛 一切MATLAB相关问题在此讨论。
回复
 
主题工具 显示模式
旧 2007-12-20, 09:23   #1
czy515
初级会员
 
注册日期: 2007-12-20
帖子: 2
声望力: 0
czy515 正向着好的方向发展
默认 【求助】求教解非线性方程参数传递问题

我是matlab的初学者,碰到一个解非线性方程的问题解决不了,希望各位高手能够帮忙解决!
具体是:我要解一个非线性方程,方程中包含了一些参数,这些参数的值需要通过主程序提供,我编了一个子程序,然后在命令窗口编了一个主程序,子程序中参数Th,Ph的值由主程序传递:
主程序:
>> Th=298.15;Ph=1;
>> Vr=fzero('fz',10)
子程序:
function f=fz(Vr)
R=83.14472;
Tc=190.6;
Pc=46.41;
Vc=R.*Tc./Pc;
Tr=Th./Tc;
Pr=Ph.*10./Pc;
a1=8.72553928e-2;a2=-7.52599476e-01;a3=3.75419887e-1;a4=1.07291342e-2;
a5=5.49626360e-3;a6=-1.84772802e-2;a7=3.18993183e-4;a8=2.11079375e-4;
a9=2.01682801e-5;a10=-1.65606189e-5;a11=1.19614546e-4;a12=-1.08087289e-4;
a13=4.48262295e-2;a14=7.53970000e-1;a15=7.71670000e-2;
f=Pr.*Vr./Tr-(1+(a1+a2./Tr.^2+a3./Tr.^3)./Vr+(a4+a5./Tr.^2+a6./Tr.^3)./Vr.^2+(a7+a8./Tr.^2+a9./Tr.^3)./Vr.^4+(a10+a11./Tr.^2+a12./Tr.^3)./Vr.^5+(a13./(Tr.^3.*Vr.^2)).*(a14+a15./Vr.^2).*exp(-a15./Vr.^2));

在命令窗口运行上述程序,总是出现如下错误:??? Error using ==> fzero
FZERO cannot continue because user supplied function_handle ==> fz
failed with the error below.

Undefined function or variable 'Th'.

另外如将参数Th,Ph的值在子程序中赋值,则运行没有问题,但这有不符合我的整个程序的计算要求。希望各位高手能够帮忙!
czy515 当前离线   回复时引用此帖
旧 2007-12-22, 16:55   #2
czy515
初级会员
 
注册日期: 2007-12-20
帖子: 2
声望力: 0
czy515 正向着好的方向发展
默认 继续求助

各位Matlab高手,我将我上次求助的问题简化了一下,希望各位高手能够帮忙指点迷津!
具体是:我要解一个非线性方程,方程中包含了一些参数,这些参数的值需要通过主程序提供。

例如:我要求解方程:ax+b/x+c=0的解,其中参数a、b、c的值必须由主程序计算得出,我变了一个子程序,通过主程序调用来解方程(当然实际程序要复杂得多):



子程序:

function f=fz(x)
f=a.*x+b./x+c;



主程序在命令窗口输入:

a=1;b=2;c=3;

x=fzero('fz',-5)



运行后总是出现如下错误提示:

??? Error using ==> fzero
FZERO cannot continue because user supplied function_handle ==> fz
failed with the error below.

Undefined function or variable 'a'.



但如果我将参数a=1;b=2;c=3;移到子程序中,变为:

子程序:

function f=fz(x)

a=1;b=2;c=3;
f=a.*x+b./x+c;



主程序在命令窗口输入:

x=fzero('fz',-5)



运行结果正常:x = -2.0000



我现在存在的问题是参数a、b、c的值是由其它程序计算而来,不可能在子程序中输入,请问怎样将参数a、b、c的值传递到子程序?或者采用其它的解非线性方程的办法?
czy515 当前离线   回复时引用此帖
旧 2008-09-09, 08:40   #3
雪魄梅魂
初级会员
 
注册日期: 2008-09-05
年龄: 38
帖子: 2
声望力: 0
雪魄梅魂 正向着好的方向发展
默认 回复: 【求助】求教解非线性方程参数传递问题

要传递的参数数放到最后面
雪魄梅魂 当前离线   回复时引用此帖
回复


发帖规则
不可以发表新主题
不可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛启用 表情符号
论坛启用 [IMG] 代码
论坛禁用 HTML 代码



所有时间均为北京时间。现在的时间是 09:47


Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.