24小时热门版块排行榜    

查看: 908  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 独孤神宇 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 298-一志愿中国农业大学-求调剂 +10 手机用户 2026-03-17 11/550 2026-03-23 16:30 by lingjue
[考研] 一志愿211 初试270分 求调剂 +3 谷雨上岸 2026-03-23 4/200 2026-03-23 15:53 by ACS Nano——
[考研] 0854 考研调剂 招生了!AI 方向 +3 pk3725069 2026-03-19 14/700 2026-03-23 14:44 by Jack_333
[考研] 333求调剂 +6 87639 2026-03-21 10/500 2026-03-23 10:41 by Iveryant
[考研] 070300化学求调剂 +5 苑豆豆 2026-03-20 5/250 2026-03-23 09:39 by guoliang1816
[考研] 070300,一志愿北航320求调剂 +3 Jerry0216 2026-03-22 5/250 2026-03-23 09:16 by 。。堂堂
[考研] 308求调剂 +3 墨墨漠 2026-03-21 3/150 2026-03-22 16:54 by i_cooler
[考研] 306求调剂 +5 来好运来来来 2026-03-22 5/250 2026-03-22 16:17 by BruceLiu320
[考研] 275求调剂 +6 shansx 2026-03-22 8/400 2026-03-22 15:27 by barlinike
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 302求调剂 +12 呼呼呼。。。。 2026-03-17 12/600 2026-03-21 17:29 by ColorlessPI
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 求调剂 +6 Mqqqqqq 2026-03-19 6/300 2026-03-21 08:04 by JourneyLucky
[考研] 化学求调剂 +4 临泽境llllll 2026-03-17 5/250 2026-03-21 02:23 by JourneyLucky
[考研] 265求调剂 +9 梁梁校校 2026-03-17 9/450 2026-03-21 02:17 by JourneyLucky
[考研] 296求调剂 +6 www_q 2026-03-18 10/500 2026-03-20 23:56 by JourneyLucky
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
信息提示
请填处理意见