24小时热门版块排行榜

返回列表

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

wangj25

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 775.7
散金: 21
红花: 1
帖子: 552
在线: 152.4小时
虫号: 348682
注册: 2007-04-18
专业: 污染控制化学

[求助] 50金币求一MATLAB的拟合程序

这是那个微分方程，有一组数据：
dS/dt=-a*S/(b+S+S的平方/c)
注意前面有一个负号。

T S
0 48.1
8 46.97
16 44.63
24 39.98
32 38.32
40 37.69
48 37.11
56 35.98
64 35.66
72 33.09
80 32.11
88 31.96
96 31.89
120 31.34
144 31.86
168 32.51
192 31.66

现在要将曲线拟合出来，求出a、b、 c 的值。并求得相关系数R。
求matlab实现程序。请写下具体程序

[ Last edited by wangj25 on 2011-8-24 at 19:42 ]

回复此楼

» 猜你喜欢

交叉科学部支持青年基金，对三无青椒是个机会吗？已经有5人回复
国家基金申请书模板内插入图片不可调整大小? 已经有6人回复
招博士已经有6人回复
限项规定已经有8人回复
国家级人才课题组招收2026年入学博士已经有5人回复
Fe3O4@SiO2合成已经有6人回复
青年基金C终止已经有4人回复
青椒八年已不青，大家都被折磨成啥样了？已经有7人回复
26申博求博导推荐-遥感图像处理方向已经有4人回复
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向已经有3人回复

1楼 2011-08-24 19:30:29

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

信彼南山

木虫 (著名写手)

应助: 33 (小学生)
金币: 4142.9
散金: 1221
红花: 16
帖子: 1178
在线: 233.5小时
虫号: 1133529
注册: 2010-10-27
专业: 导航、制导与传感技术

【答案】应助回帖

★
wangj25(金币+40): 已经按照你的方法弄了，你的思路是完全正确的，是我的数据有点问题，谢谢了 2011-08-25 22:02:08
xiegangmai(金币+1): 辛苦了！ 2011-08-26 07:09:33

CODE:

% wangj25 

clear;clc;

%% 初始数据

z=[0 48.1 

8 46.97 

16 44.63 

24 39.98 

32 38.32 

40 37.69 

48 37.11 

56 35.98 

64 35.66 

72 33.09 

80 32.11 

88 31.96 

96 31.89 

120 31.34 

144 31.86 

168 32.51 

192 31.66 ];

%% main

T=z(:,1);

S=z(:,2);

figure(1);plot(T,S);

h=1;

T0=0;TT=192;

T1=T0:h:TT;

S1=spline(T,S,T1);

z1=[T1;S1];

figure(2);plot(T1,S1);

dsdt=[];

for i=1:TT-1

    dsdt(i)=(S1(i+2)-S1(i))/(T1(i+2)-T1(i));

end

t=T1(2:TT);

figure(3);plot(t,dsdt);