24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

» 猜你喜欢

垃圾破二本职称评审标准已经有19人回复
职称评审没过，求安慰已经有53人回复
毕业后当辅导员了，天天各种学生超烦已经有5人回复
26申博自荐已经有3人回复
A期刊撤稿已经有4人回复

1楼 2011-05-27 21:11:31

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

695561947

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 331
帖子: 14
在线: 12.5小时
虫号: 2139944
注册: 2012-11-21
专业: 信息处理方法与技术

引用回帖:

411328楼: Originally posted by lds588 at 2011-05-28 00:27:57
clear,clc
p=;
c=;
cc=c';
pp=p';
F=@(q,x)(x(:,1)-x(:,2)/q(1)).*(10.^(q(2)*(x(:,1)-x(:,2)/q(1)))+q(3))-x(:,2);
x=;

q0=;
warning off
q=nlinfit(x,zeros(size(x,1),1),F,q0);
disp('H A B分别为 ...

你好，我想拟合下微分方程的系数
已知函数形式dx/dt=x*(1+x^m)(1+x^n)
已知大量（x,t）数据
用Matlab拟合系数m n，因为急用来不及研究，能不能帮忙写下，谢谢

赞一下

回复此楼

16楼2013-01-09 17:33:21

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

查看全部 17 个回答

lds588

铁杆木虫 (著名写手)

愚钝

应助: 15 (小学生)
金币: 10949.6
散金: 2018
红花: 5
帖子: 1444
在线: 206.3小时
虫号: 890720
注册: 2009-11-01
专业: 凝聚态物性I:结构、力学和

【答案】应助回帖

★ ★ ★ ★ ★
臭水沟(金币+5): 给力！给5个金币！理由：我们就需要像你这样的热心的虫子！ 2011-05-28 09:07:21

CODE:

clear,clc

p=[0  0.001  0.002  0.003  0.004  0.005  0.006  0.008  0.011  0.014  0.019  0.023  0.03  0.037  0.048  0.06  0.08  0.1  0.136  0.17  0.217  0.279  0.388  0.444];

c=[0  0.307598804  0.538975997  0.777628765  0.959350561  1.069744162  1.190592446  1.340828541  1.512368753  1.643259589  1.769609554  1.903576465  2.027165658  2.151392907  2.294367083  2.425558327 2.598588163  2.723808495  2.912380082  3.049909679  3.206076964 3.386393547  3.644346885  3.758352826];

cc=c';

pp=p';

F=@(q,x)(x(:,1)-x(:,2)/q(1)).*(10.^(q(2)*(x(:,1)-x(:,2)/q(1)))+q(3))-x(:,2);   

x=[cc pp];

q0=[1    1     1];

warning off

q=nlinfit(x,zeros(size(x,1),1),F,q0);

disp('H A B分别为')

disp(num2str(q));

plot(x(:,1),x(:,2),'r^');hold on; 

ezplot(@(x,y)F(q,[x,y]),[0,3.8,-0.1,0.5]);