24小时热门版块排行榜    

查看: 918  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 独孤神宇 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料学求调剂 +6 Stella_Yao 2026-03-20 6/300 2026-03-25 00:37 by baoball
[考研] 289求调剂 +8 硕星赴 2026-03-23 8/400 2026-03-24 20:17 by peike
[考研] 085602 289分求调剂 +5 WWW西西弗斯 2026-03-24 5/250 2026-03-24 18:51 by jhhcooi
[考研] 材料专硕331求调剂 +4 鲜当牛 2026-03-24 4/200 2026-03-24 15:58 by JourneyLucky
[考研] 299求调剂 +7 某某某某位 2026-03-21 7/350 2026-03-24 15:24 by cuifj
[考研] 293求调剂 +6 加一一九 2026-03-24 6/300 2026-03-24 14:29 by JourneyLucky
[考研] 279分求调剂 一志愿211 +18 chaojifeixia 2026-03-19 20/1000 2026-03-24 10:34 by dolphin_ycj
[考研] 工科0856求调剂 +5 沐析汀汀 2026-03-21 5/250 2026-03-23 17:56 by 海瑟薇-
[考研] 一志愿070300浙大化学358分,求调剂! +4 酥酥鱼.. 2026-03-21 4/200 2026-03-23 08:12 by Iveryant
[考研] 311求调剂 +6 冬十三 2026-03-18 6/300 2026-03-22 20:18 by edmund7
[考研] 一志愿华中农业071010,总分320求调剂 +5 困困困困坤坤 2026-03-20 6/300 2026-03-22 17:41 by hxsm
[考研] 354求调剂 +7 Tyoumou 2026-03-18 10/500 2026-03-22 11:11 by 人来盛
[考研] 085600材料与化工306 +4 z1z2z3879 2026-03-21 4/200 2026-03-21 23:44 by ms629
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-18 3/150 2026-03-21 00:46 by JourneyLucky
[考研] 22408 344分 求调剂 一志愿 华电计算机技术 +4 solanXXX 2026-03-20 4/200 2026-03-20 23:49 by alg094825
[考研] 南京大学化学376求调剂 +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
信息提示
请填处理意见