| 查看: 1055 | 回复: 3 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
独孤神宇版主 (知名作家)
|
[求助]
1stOpt 拟合的曲线为什么和MATLAB里面绘出来的不一样？ 已有2人参与
|
||
|
1stOpt 拟合得到的参数，在MATLAB里面绘制出来图形为什么不一样？ SharedModel; Variable x,y(2);//y1,y2; Parameter a=[-2,2],b=[-1,1],c=[-2,2],d=[-1,1]; Constant epsilon=0.3921,h0=6.3,R=125,Sf=0.04558,VR=1308.33; ConstStr h1=h0*(1-epsilon),Delta_h=h0*epsilon, L=sqrt(R*Delta_h-(Delta_h^2)/4), z=Delta_h*x^2/(2*L^2)+h1/2,V1=VR*(1+Sf),Vx=h1*V1/(2*z),V0=h1*V1/h0, delta_f=15+9.2*((h0-2*z)/h0)^0.11, p1=(exp(2*b*(a*Vx-V0)*x/z^2))*(4/3*b*delta_f*Delta_h/((L^2)*z^3)*(V0-a*Vx)*x^3+delta_f*Delta_h*x^2/(z*L^2)+delta_f*(1-2*Delta_h/h0)); Function y1=p1; y2=delta_f*Delta_h*x/(L^2)+2*d*(c*Vx-V0)*y1/z; Data; //x,y1,y2 1.25 66.16352201 32.45283019 2.5 99.87421384 34.96855346 3.75 86.28930818 41.00628931 5 80.75471698 25.91194969 6.25 76.72955975 21.88679245 7.5 67.67295597 18.36477987 8.75 60.62893082 1.257861635 10 56.60377358 -15.8490566 11.25 54.59119497 -28.42767296 12.5 53.58490566 -28.93081761 13.75 54.08805031 -33.96226415 15 48.55345912 -20.37735849 16.25 44.02515723 -10.81761006 17.5 45.03144654 -5.283018868 —————————————————————————————————— MATLAB clear; h0=6.32;epsilon=0.304;R=125;VR=1046.67;Sf=0.03456;h1=h0*(1-epsilon); Delta_h=h0-h1;L=sqrt(R*Delta_h-(Delta_h^2)/4); x=0:0.01:L;b=0.0012; z=Delta_h*x.^2/(2*L^2)+h1/2; Vx=h1*VR*(1+Sf)./(2*z);V1=(1+Sf)*VR;V0=h1*V1/h0; delta_f=10+9*(((h0-2*z)/h0).^0.06); a=0.85;b=0.004; p2=(exp(2*b*(a*Vx-V0).*x./z.^2)).*(4/3*b*delta_f*Delta_h./(L^2*z.^3).*(V0-a*Vx).*x.^3+delta_f*Delta_h.*x.^2./(L^2*z)+delta_f*(1-2*Delta_h/h0)); plot(x,p2) c=0.68;d=0.014; t2=delta_f*Delta_h.*x/(L^2)+2*d*(c*Vx-V0).*p2./z; hold on plot(x,t2)  1.jpg  2.jpg  3.jpg |
回复此楼» 猜你喜欢
93BebMhtakh前后11位开头都是大写
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有86人回复
-
感觉是下周放榜了
已经有11人回复
-
filecode=后面第一个是大写字母
已经有10人回复
-
时间戳变了，能看出什么问题？
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
lipenggg
专家顾问 (正式写手)
-
专家经验: +65 - 应助: 122 (高中生)
- 金币: 4841.3
- 红花: 43
- 帖子: 722
- 在线: 435.7小时
- 虫号: 3929201
- 注册: 2015-06-17
- 性别: GG
- 专业: 流体力学
- 管辖: 计算模拟
3楼2015-08-26 17:08:45
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68166.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.3小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2015-08-26 16:35:54
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2015-08-26 17:17:54
