题目
帮忙用Matlab拟合下面的方程
P1 = linspace(-4,4,50);
T1 = (1-P1+2*P1.^2).*exp(-.5*P1.^2)
尽量用最小二乘拟合
如果正确还有追加的100
P1 = linspace(-4,4,50);
T1 = (1-P1+2*P1.^2).*exp(-.5*P1.^2)
尽量用最小二乘拟合
如果正确还有追加的100
提问时间:2021-11-26
答案
直接用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)
正着做:给一个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)
举一反三
已知函数f(x)=x,g(x)=alnx,a∈R.若曲线y=f(x)与曲线y=g(x)相交,且在交点处有相同的切线,求a的值和该切线方程.
我想写一篇关于奥巴马的演讲的文章,写哪一篇好呢?为什么好
奥巴马演讲不用看稿子.为什么中国领导演讲要看?
想找英语初三上学期的首字母填空练习……
英语翻译
最新试题
- 1英语翻译
- 2比较根号3+根号2与根号2乘根号5的大小,并说明理由
- 3如图所示,一很长的、不可伸长的柔软轻绳跨过光滑定滑轮(定滑轮质量不计),绳两端各系一小球a和b.a球质量为m,静置于地面;b球质量为3m,用手托住,高度为h,此时轻绳刚好拉紧.从静
- 4若|a+3|-|b+1|=0则a-b+1/2
- 5英语交通拥挤怎么写
- 6英语翻译
- 7杜甫生活的时代特点是什么?这在他的诗中有何表现?
- 8一根原长为l0=0.1m的轻弹簧,一端拴住质量为m=0.5kg的小球,以另一端为圆心在光滑的水平面上做匀速圆周运动,如图所示,角速度为ω=10rad/s,弹簧的劲度系数k=100N/m,求:小球做匀速
- 9化简20的结果是( ) A.52 B.25 C.210 D.45
- 10酱油为什么不像水一样会结冰?
热门考点
- 1are there any pandas in the mountains?和are there any fish in the rivers? 中文是什么
- 2父子俩早晨从家跑步去公园晨练,儿子每小时9千米,父亲每小时15千米.儿子出发半小时父亲才出发,结果二人
- 3notice什么意思
- 4甲乙两辆车的速度比为5:7,两车同时从甲乙两地出发,相向而行,两车的相遇地点距离中点80km甲乙两地距?
- 5请在括号里填上一个3个字母的单词:I like bananas ( ) ice cream.
- 6已知发动机的排放量为3升,燃气对活塞的平均压强为8×10的五次方帕斯卡,燃气对活塞
- 7英语:The colour of the carrots are red.有哪里错了?
- 8一道简单的判断题
- 9同一样工作,小李用了5个小时,小张用了4个小时.请问小李的效率比小张的慢了几分之几?
- 10将一张足够大的报纸对折51次后其厚度会超过地球到太阳的距离,请问这个结果是怎么算出来的?