| 查看: 602 | 回复: 4 | ||||
[求助]
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 |
回复此楼» 猜你喜欢
29号明天会评吗
已经有4人回复
-
国社科又开始会评了，不知道这次命运如何
已经有10人回复
-
梦想
已经有4人回复
-
中青基了要发朋友圈吗？
已经有6人回复
-
基金不中，共勉
已经有12人回复
-
能否申诉？
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上意见出来了
已经有19人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 各位前辈，利用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人回复
2楼2014-12-16 18:19:52
|
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
|
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
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
10