24小时热门版块排行榜    

查看: 1031  |  回复: 15

月只蓝

主管区长 (职业作家)

引用回帖:
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

铁杆木虫 (职业作家)

【答案】应助回帖

深有同感,自己问问题都描述不清或模棱两可,让他人还得化时间去猜测。求助应该先学会如何简单明了地把自己的问题说清楚,大家不可能像你一样了解该问题及其背景。
16楼2015-08-20 16:00:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 夜雨tju 的主题更新
信息提示
请填处理意见