版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

18704487833

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 57
帖子: 26
在线: 8.2小时
虫号: 2674109
注册: 2013-09-23

[求助] 1stopt程序，求修改

Title "Type your title here";
LoopConstant r=[0:0.02:1];
Constant kr=1.6,lamda0=3,u=0.3,theta0=0.1,theta1c=0.12,theta1s=0.15;//Parameters ;
Variable t=[0.0:3.14],b=0,b'=1;
ODEOptions=[SS=0.01,A=0,P=5];
Plot b[t], b';
ChartType =3;
ODEFunction b''=((r*cos(b)+u*sin(t))^2+(b'*r+(lamda0*(1+kr*r*cos(t)))*cos(b)+u*sin(beta)*cos(t))^2)*sin(2*(theta0+theta1c*cos(t)+theta1s*sin(t))-2*(atan((b'*r+(lamda0*(1+kr*r*cos(t)))*cos(b)+u*sin(beta)*cos(t))/((r*cos(b)+u*sin(t))))))/(13.04651163*0.04*(3*r^2+0.02^2/4)/3.3/0.5273)-sin(b)*cos(b); //Function ;
这个程序哪里有错误？求大神指点一二

捕获.PNG

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

各位前辈，利用1stOpt对非线性方程拟合每次运算结果都不一样怎么办呢？已经有9人回复
看了下1stOpt的介绍，感觉略显浮夸已经有14人回复
求大哥大姐帮跑一下1stopt~新手啊，版本低啊已经有6人回复
1stopt自变量超多时是否有简单语法？AND 自变量多数据量大的话1stopt都会算的特慢么已经有7人回复
求大神跑下1stOpt，一个sharemodel，谢谢！已经有7人回复
求高手给做一道题，用MATLAB或1stOpt，请把结果直接给出来，我的软件就不行，谢谢已经有5人回复
求助1stopt含约束非线性拟合已经有5人回复
1stopt解方程组已经有12人回复
求回归方程参数的最优化算法matlab代码已经有17人回复
1stopt求最小值已经有37人回复
一个1stOpt的初级问题：复数方程组计算问题（有代码）已经有5人回复
rmse，sse计算求助已经有3人回复
1stopt中使用sharedmodel拟合曲线出现问题，求指示！！已经有6人回复
求大哥大姐帮跑一下1stopt~新手啊，版本低啊已经有9人回复
1stOPT1.5版每次仿真结果不一样，求助！急！已经有3人回复
运算过程有复数出现，求高版本1stOpt 已经有6人回复
大家在使用1stOpt进行数据拟合的时候有没有遇到过这种情况？已经有8人回复
学1stopt要多久？已经有12人回复
求回归方程参数估计的最优化算法matlab代码已经有9人回复
求高人指点用matlab求解非线性方程组，解决了追加100金币；已经有11人回复

1楼 2014-12-14 08:09:04

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

雕刻你的拾光

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 143.4
帖子: 34
在线: 15小时
虫号: 3596576
注册: 2014-12-16
专业: 病原真菌、真菌感染与宿主

//Parameters a, b, u,;
//Variables t, x,;
//QuickReg ;
//Function  x=a*b*exp(u*t)/(b-a+a*exp(u*t));
//RowData ;
      t=0,4,8,12,16,20,24,28,32,36;
      x=0.05,0.16,0.6515,1.5426,2.4047,3.1066,3.6693,3.8165,3.8874,3.9618;
我想求出a b u 的最佳估算值，求帮忙指点一下

赞一下

回复此楼

2楼2014-12-16 18:19:52

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

雕刻你的拾光

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 143.4
帖子: 34
在线: 15小时
虫号: 3596576
注册: 2014-12-16
专业: 病原真菌、真菌感染与宿主

Parameters X0, Xm, um,a,b,P0,Yx,Yp,m,S0;
Variables t, X, P, S;
Function X=X0*Xm*exp(um*t)/(Xm-X0+X0*exp(um*t))
         P=P0+a*(X0*Xm*exp(um*t)/(Xm-X0+X0*exp(um*t))-X0)+(b*Xm/um)*ln((Xm-X0+X0*exp(um*t))/Xm)
         S=S0-(1/Yx+a/Yp)*(X0*Xm*exp(um*t)/(Xm-X0+X0*exp(um*t)-X0)-Xm/um*(b/Yp+m)*ln((Xm-X0+X0*exp(um*t))/Xm);
Data;
t X       P       S
0 0.05       0       79
4 0.16       0       78
8 0.6515    0       76
12  1.5426    1.5    74
16  2.4047    4.9    69
20  3.1066    8.7    60
24  3.6693    12.4    50
28  3.8165    16.1    43
32  3.8874    17.2    40
36  3.9618    17.9    38
就是这个，方程应该没啥问题，求专家帮我看看

赞一下

回复此楼

3楼2014-12-16 21:37:23

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

雕刻你的拾光

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 143.4
帖子: 34
在线: 15小时
虫号: 3596576
注册: 2014-12-16
专业: 病原真菌、真菌感染与宿主

Parameters X0, Xm, um,a,b,P0,Yx,Yp,m,S0;
Variables t, X, P, S;
Function X=X0*Xm*exp(um*t)/(Xm-X0+X0*exp(um*t))
         P=P0+a*(X0*Xm*exp(um*t)/(Xm-X0+X0*exp(um*t))-X0)+(b*Xm/um)*ln((Xm-X0+X0*exp(um*t))/Xm)
         S=S0-(1/Yx+a/Yp)*(X0*Xm*exp(um*t)/(Xm-X0+X0*exp(um*t)-X0)-Xm/um*(b/Yp+m)*ln((Xm-X0+X0*exp(um*t))/Xm);
Data;
t X       P       S
0 0.05       0       79
4 0.16       0       78
8 0.6515    0       76
12  1.5426    1.5    74
16  2.4047    4.9    69
20  3.1066    8.7    60
24  3.6693    12.4    50
28  3.8165    16.1    43
32  3.8874    17.2    40
36  3.9618    17.9    38
三个函数表达式分开算没问题，在一起算就说有误，什么情况啊

赞一下

回复此楼

4楼2014-12-16 21:40:52

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

lipenggg

专家顾问 (正式写手)

专家经验: +65
应助: 122 (高中生)
金币: 4841.3
红花: 43
帖子: 722
在线: 435.7小时
虫号: 3929201
注册: 2015-06-17
性别: GG
专业: 流体力学
管辖: 计算模拟

Parameters X0, Xm, um,a,b,P0,Yx,Yp,m,S0;
Variables t, X, P, S;
Function X=X0*Xm*exp(um*t)/(Xm-X0+X0*exp(um*t));
         P=P0+a*(X0*Xm*exp(um*t)/(Xm-X0+X0*exp(um*t))-X0)+(b*Xm/um)*ln((Xm-X0+X0*exp(um*t))/Xm);
         S=S0-(1/Yx+a/Yp)*(X0*Xm*exp(um*t)/(Xm-X0+X0*exp(um*t)-X0)-Xm/um*(b/Yp+m)*ln((Xm-X0+X0*exp(um*t))/Xm));
Data;
t X       P       S
0 0.05       0       79
4 0.16       0       78
8 0.6515    0       76
12  1.5426    1.5    74
16  2.4047    4.9    69
20  3.1066    8.7    60
24  3.6693    12.4    50
28  3.8165    16.1    43
32  3.8874    17.2    40
36  3.9618    17.9    38

赞一下

回复此楼

一团和气

5楼2015-06-19 16:45:34

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

相关版块跳转我要订阅楼主 18704487833 的主题更新

返回列表