24小时热门版块排行榜    

Znn3bq.jpeg
查看: 927  |  回复: 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)

1stOpt 拟合的曲线为什么和MATLAB里面绘出来的不一样?
1.jpg


1stOpt 拟合的曲线为什么和MATLAB里面绘出来的不一样?-1
2.jpg


1stOpt 拟合的曲线为什么和MATLAB里面绘出来的不一样?-2
3.jpg
回复此楼

» 猜你喜欢

数值计算
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
独孤神宇: 金币+3, 有帮助, 谢谢,已近发现问题所在了! 2015-08-26 17:05:16
独孤神宇: 金币+2 2015-11-15 19:07:17
在保证两种软件得到的拟合效果相同的前提下,输出的图形上肯定也是相似甚至完全相同的。否则,说明程序中有问题。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2015-08-26 16:35:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

1stopt是拟合效果图,如果针对公式进行1stopt专门画图取足够的点,也可以平滑。
一团和气
3楼2015-08-26 17:08:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
独孤神宇: 金币+5, 有帮助 2015-08-26 17:30:16
应该是一样的,1stOpt画的是折线(对应x的点数,14),Matlab画的点数多,显的光滑而已。
4楼2015-08-26 17:17:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 独孤神宇 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085801电气专硕272求调剂 +19 电气李 2026-04-13 21/1050 2026-04-15 13:37 by 黑科技矿业
[考研] 297工科调剂? +14 河南农业大学-能 2026-04-13 15/750 2026-04-15 13:25 by 黑科技矿业
[考研] 085404 22408 309分求调剂 +8 lzmk 2026-04-14 9/450 2026-04-15 13:23 by 黑科技矿业
[考研] 296求调剂 +10 汪!?! 2026-04-09 11/550 2026-04-15 09:17 by fenglj492
[考研] 考研求调剂 +6 ban班小七 2026-04-11 6/300 2026-04-14 14:06 by 哆啦A梦只是个梦
[考研] 一志愿哈工大 085600 277 12材科基求调剂 5+5 chenny174 2026-04-10 37/1850 2026-04-14 07:39 by Abskk
[考研] 一志愿211 0703化学 346分求调剂 +26 土豆er? 2026-04-09 29/1450 2026-04-13 15:15 by 独醉梦孤城
[找工作] 山东高校教师考核超级无底线,员工过不下去啦 +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[考研] 086003调剂求助 +21 苏弋万 2026-04-09 22/1100 2026-04-11 20:25 by dongdian1
[考研] 284求调剂 +11 archer.. 2026-04-09 12/600 2026-04-11 20:23 by 蓝云思雨
[考研] 求调剂 +11 翩翩一书生 2026-04-09 11/550 2026-04-11 19:57 by 逆水乘风
[考研] 求调剂 +6 archer.. 2026-04-09 8/400 2026-04-11 10:55 by zhq0425
[考研] 22408 352分求调剂0854类 +4 努力的夏末 2026-04-09 4/200 2026-04-11 09:57 by zhq0425
[考研] 0854调剂 +4 长弓傲 2026-04-09 4/200 2026-04-11 09:18 by 猪会飞
[考研] 337求调剂 +4 研s. 2026-04-10 4/200 2026-04-11 08:57 by zhq0425
[考研] 085506-求调剂-285分 +3 雷欧飞踢 2026-04-08 3/150 2026-04-11 08:37 by zhq0425
[考研] 265求调剂 +12 风说她早忘了 2026-04-10 13/650 2026-04-10 18:56 by chemisry
[考研] 一志愿矿大,材料工程专硕314分,0856可调都可以 +15 无懈可击的巨人 2026-04-09 15/750 2026-04-10 18:10 by hmn_wj
[考研] 求调剂 材料与工程 324分 专硕 +19 翩翩一书生 2026-04-10 21/1050 2026-04-10 11:41 by wp06
[考研] 调剂 +19 2261744733 2026-04-08 19/950 2026-04-09 19:11 by vgtyfty
信息提示
请填处理意见