应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » matlab求解非线性方程组
162/2返回列表上一页12
查看: 1268  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

月只蓝

主管区长 (职业作家)
引用回帖:
7楼: Originally posted by 夜雨tju at 2015-08-16 20:11:03
我是化工类的，不太擅长这个，大神可不可以指点细一些，我自己去查资料也可以，大神觉得什么函数可以，我之前也看了智能算法的书。

8楼的dingd是这类问题的大神，求助于他，能让你事半功倍。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
11楼2015-08-19 10:25:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

夜雨tju

金虫 (小有名气)
引用回帖:
9楼: Originally posted by dingd at 2015-08-16 22:38:40
抱歉，应该是“已知常数”

你好，我刚刚才看到你的回复，一直常数很多是公式算出来的，这里面我用了两个循环
赞一下 回复此楼
12楼2015-08-19 21:16:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

夜雨tju

金虫 (小有名气)
引用回帖:
9楼: Originally posted by dingd at 2015-08-16 22:38:40
抱歉，应该是“已知常数”

global K CA0 CB0 CC0 CD0 k1 U L V xa1 xb1 xc1 xd1 ya1 HETP1 HETP2 HETP3 HETP4 mi M1 M2 M3 M4 K1
L=1.168;
V=0.7642;
xa1=0.07039;
xb1=0.07198;
xc1=0.85393;
xd1=0.0037;
ya1=0.02176;
yb1=0.02531;
yc1=0.92063;
yd1=0.03229;
T=370.72;
HETP1=0.4;
HETP2=0.2;
HETP3=0.42;
HETP4=0.2;
mi=803956.8;            
mig=922.45;            
M1=60;
M2=88;
M3=130;
M4=18;
Ts=40;               
lambda=286.76;        
wi1=0.5;                  
wi2=1;
Le=0.02772;
CA0=465.25;
CB0=475.76;
CC0=5644.11;
CD0=24.46;
k1=0.00000206;
K=3.621;
K1=0.85;
K2=0.9;
K3=0.95;
g=9.8;
for m=1:0.5:2;
    for R=0.01:0.01:0.4;
                S=pi*R^2;  
                a=acos((1+2*m)/(3+2*m));
                S1=a*R^2-((1+2*m)*(R^2)*sin(a))/(3+2*m);
                S5=a*R^2-((1+2*m)*(R^2)*sin(a))/(3+2*m);
                b=acos(1/(3+2*m));
                S2=(b-a)*R^2+((1+2*m)*R^2*sin(a)-(R^2)*sin(b))/(3+2*m);
                S4=(b-a)*R^2+((1+2*m)*R^2*sin(a)-(R^2)*sin(b))/(3+2*m);
                S3=(1-2*b/pi)*pi*R^2+(2*R^2/(3+2*m))*sin(b);
                U=(xa1*M1+xb1*M2+xc1*M3+xd1*M4)*L/(mi*pi*R^2);
                UG=(ya1*M1+yb1*M2+yc1*M3+yd1*M4)*V/(mig*pi*R^2);
                U1=(ya1*M1+yb1*M2+yc1*M3+yd1*M4)*V/(mig*pi*R^2);
                U3=(ya1*M1+yb1*M2+yc1*M3+yd1*M4)*V/(mig*pi*R^2);
                U5=(ya1*M1+yb1*M2+yc1*M3+yd1*M4)*V/(mig*pi*R^2);
                A2=S2+S4;
                A4=S1+S3+S5;
                U2=S*U1/A2;
                U4=S*U1/A4;
                L1=3600*U;
                F=U2*sqrt(mig/1000);
                Lc1=Ts*U/2;
                dpz=(7.38*10^(4.46*L1*10^(-3)))*(UG*sqrt(mig/1000))^(1.72+L1*3.8*10^(-3));
                hL1=(dpz*Lc1*1000+lambda*(Lc1+Le)*U*mi-Lc1*mi)/(mi-dpz*1000);
                h1=Lc1+hL1;
                Lc2=Lc1;
                h3=h1;
                h2=HETP2*(0.5-h1/HETP1);
                h4=h2;
                lin=L1/S;
赞一下 回复此楼
13楼2015-08-19 21:17:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็
引用回帖:
13楼: Originally posted by 夜雨tju at 2015-08-19 21:17:31
global K CA0 CB0 CC0 CD0 k1 U L V xa1 xb1 xc1 xd1 ya1 HETP1 HETP2 HETP3 HETP4 mi M1 M2 M3 M4 K1
L=1.168;
V=0.7642;
xa1=0.07039;
xb1=0.07198;
xc1=0.85393;
xd1=0.0037;
ya1=0.02176;
yb1=0.02531 ...

代码没贴全吧？
两个for 循环都不封闭？？
赞一下 回复此楼
ส็็็็็็็็็็็็็็็็็็็็
14楼2015-08-20 13:43:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็
引用回帖:
13楼: Originally posted by 夜雨tju at 2015-08-19 21:17:31
global K CA0 CB0 CC0 CD0 k1 U L V xa1 xb1 xc1 xd1 ya1 HETP1 HETP2 HETP3 HETP4 mi M1 M2 M3 M4 K1
L=1.168;
V=0.7642;
xa1=0.07039;
xb1=0.07198;
xc1=0.85393;
xd1=0.0037;
ya1=0.02176;
yb1=0.02531 ...

看了你的代码，让人更加困惑了，你既然给常数赋值都是标量，为什么放在嵌套循环里面做无用功？

如果这段代码正确的写法是后面缺少两个end

那么你的代码相当于把 两个 end 去掉之后

for m=1:0.5:2;
    for R=0.01:0.01:0.4;

语句换成
m=2; R=0.4;

效果相同，而且少做了循环嵌套中的很多无意义的计算。

建议你不要用代码，而是用word文档和公式的形式，老老实实详细表达自己的问题，否则别人很难理解你的问题是什么
赞一下 回复此楼
ส็็็็็็็็็็็็็็็็็็็็
15楼2015-08-20 14:13:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

深有同感，自己问问题都描述不清或模棱两可，让他人还得化时间去猜测。求助应该先学会如何简单明了地把自己的问题说清楚，大家不可能像你一样了解该问题及其背景。
赞一下(1人) 回复此楼
16楼2015-08-20 16:00:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 夜雨tju 的主题更新
162/2返回列表上一页12
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 a2tycdlnq1 2026-05-23 5/250 2026-05-24 16:21 by hhx1yx9evi
[基金申请] 河北省自然科学基金 +6 Peterchao 2026-05-18 9/450 2026-05-24 16:02 by 130067131
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 pmo95bazuy 2026-05-23 8/400 2026-05-24 15:56 by 1uy1ht2y9r
[基金申请] 西安交大新媒学院副院长用撤稿论文结题 +3 bjvtcliu 2026-05-24 5/250 2026-05-24 10:16 by kudofaye
[教师之家] 论文撤稿了 +3 bjvtcliu 2026-05-24 5/250 2026-05-24 10:06 by Equinoxhua
[教师之家] 某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！ +4 zju2000 2026-05-21 5/250 2026-05-24 09:35 by songwz
[考博] 26/27申博自荐 10+4 ZXW0202 2026-05-22 9/450 2026-05-24 08:47 by bjvtcliu
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 hvkbtfonbv 2026-05-23 3/150 2026-05-24 08:01 by 9ps9vgkqva
[基金申请] 揭秘青基评审内幕：几个A才能顺利中标 +3 国自然国社科中 2026-05-23 4/200 2026-05-23 15:37 by 2000zf36392
[考博] 博士申请 +3 焦晓明 2026-05-21 3/150 2026-05-23 11:26 by mlc840311
[文学芳草园] 献血感触 +7 呀呀好傻 2026-05-19 13/650 2026-05-21 20:15 by 呀呀好傻
[基金申请] 面上本子正文33页，违规吗？会被低分嘛？ +14 1234567wang 2026-05-17 16/800 2026-05-21 17:58 by 脆脆的饼干
[基金申请] 国自然评分 +4 无名者登山 2026-05-20 5/250 2026-05-21 16:35 by swuq
[基金申请] 国自然上会要求 +7 无名者登山 2026-05-18 11/550 2026-05-21 15:50 by draco1987
[基金申请] 提交了我也来说说感想 +9 fummck 2026-05-20 10/500 2026-05-21 14:17 by draco1987
[基金申请] 评审有感 +15 popular289 2026-05-18 26/1300 2026-05-21 10:35 by 西葫芦炒鸡蛋
[有机交流] 反应很差，大量原料没有反应 5+3 Mr.Zot 2026-05-19 8/400 2026-05-20 22:19 by Equinoxhua
[考博] 如果工作了想读博，可以边工作边读全日制嘛？ 30+3 铁达火车 2026-05-18 5/250 2026-05-20 09:33 by tfang
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
[硕博家园] 我在等一个没有答案的答案 +3 Love_MH 2026-05-17 3/150 2026-05-18 02:22 by 竹林孤影
信息提示
关闭
请填处理意见
关闭