帮忙用Matlab拟合下面的方程
帮忙用Matlab拟合下面的方程 P1 = linspace(-4,4,50); T1 = (1-P1+2*P1.^2).*exp(-.5*P1.^2)尽量用最小二乘拟合如果正确还有追加的100
最佳回答
直接用cftool就可以看到曲线拟合的界面了,输入数据后就可以选择拟合的方法了,直接可以得到系数。正着做:给一个r和xm,通过ode45求得t,x,然后再与你的数据对比。最终选择一个合适的r和xm。就要一直变参数,我觉得比较难。至少来说r和xm的选择范围太大了。而且ode45得到的t,x与你的实验t,x肯定不是同一个t下面的数据,也不好比较。不过基于你的方程比较简单,我们可以直接解出它的通解来:>> xx=dsolve('Dx/x=r-r/b*x') 得到:xx = b/(1+exp(-r*t)*C1*b) 自己想手算的话,见附录。附录:xm我用b代替 (dx/dt)/x/(r-r/b*x)=1 dx/[rx(1-1/b*x)]=dt [1/x+(1/b)/(1-x/b)]dx=r*dt 积分得 ln(x)-ln(1-x/b)=r*t+C -ln[(1-x/b)/x]=r*t+C ln[1/x-1/b]=-r*t-C 两边指数得 1/x-1/b=C1*exp(-r*t) 1/x=1/b+C1*exp(-r*t) x=b/[1+b*C1*exp(-r*t)] 即 x=b/(1+exp(-r*t)*C1*b)
最新回答共有2条回答
-
2026-03-30 17:07:40糟糕的巨人
回复直接用cftool就可以看到曲线拟合的界面了,输入数据后就可以选择拟合的方法了,直接可以得到系数。正着做:给一个r和xm,通过ode45求得t,x,然后再与你的数据对比。最终选择一个合适的r和xm。就要一直变参数,我觉得比较难。至少来说r和xm的选择范围太大了。而且ode45得到的t,x与你的实验t,x肯定不是同一个t下面的数据,也不好比较。不过基于你的方程比较简单,我们可以直接解出它的通解来:>> xx=dsolve('Dx/x=r-r/b*x') 得到:xx = b/(1+exp(-r*t)*C1*b) 自己想手算的话,见附录。附录:xm我用b代替 (dx/dt)/x/(r-r/b*x)=1 dx/[rx(1-1/b*x)]=dt [1/x+(1/b)/(1-x/b)]dx=r*dt 积分得 ln(x)-ln(1-x/b)=r*t+C -ln[(1-x/b)/x]=r*t+C ln[1/x-1/b]=-r*t-C 两边指数得 1/x-1/b=C1*exp(-r*t) 1/x=1/b+C1*exp(-r*t) x=b/[1+b*C1*exp(-r*t)] 即 x=b/(1+exp(-r*t)*C1*b)
热门文章
- 康达学院专转本五年制
- 高考一个考场分ab卷吗
- not only but also用法
- 某物体做自由落体运动,从释放开始计时,则物体在前2s内的平均速度为______m/s,物体下落2m时的速度大小为______m/s.
- 三角函数公式大全表格
- 地理中考必背知识点2022
- 2013-2014学年小学六年级科学上学期期末考试试卷及答案
- 人教版2014-2015学年小学五年级英语第二学期期中教学质量检测试卷及答案
- 【Linux驱动开发】设备树详解(二)设备树语法详解
- 别跟客户扯细节
- 在别的城市买房子能落户吗
- 卖房前要把装修贷还完吗
- 高中政治教学提高教学效果的方法探究
- “互联网+”背景下的初中英语课堂教学改革与创新策略研究
- 2022年终止合同范本
- 租房合同范本范文
- 如何挑选土豆
- 如何挑选土鸡
