[求助]怎么用ＭＡＴＬＡＢ解Logistic模型 - 编程论坛

#2

wangweicupl2007-09-22 06:34

　　　　如果各位高手有空能不能告诉我怎么解：t=[0 1 2 3 4 5];x=[111 117 122 127 134 139];求递增常数Ｋ和最大理论值Ｍ

#3

wangweicupl2007-09-22 06:36

　　我知道应该很简单的，各位就帮帮忙告诉我怎么解吧．

#4

ydgsl2007-09-22 09:05

这样的数据应该是无法解的。你若是用y=M/(1+a*exp(-K*t))进行拟合，其t取值不当。t=0时，y值是起始生长值。y值犹如人的体重增长过程，随着时间t的增加而增长。当t=0时，该婴儿已达111（斤），这是完全不可能的。另外，数据点也太少，不能很好地反映整个过程。

#5

wangweicupl2007-09-22 09:22

回复：（ydgsl）这样的数据应该是无法解的。你若是用y...

哦？

果然是个高手啊！
那这组数据能解吗？
t=[1 2 3 4 5 6 ];
x=[12.767 12.845 12.9227 12.9988 13.0756];
还是用的Ｌogistic来解的．

#6

wangweicupl2007-09-22 09:27

再加一句：求递增常数Ｋ和最大理论值Ｍ
本人不盛感激啊．

#7

ydgsl2007-09-22 09:41

对于前一组数据，结果如下：
SSy = 550
[M a K] = [345.4143660, 2.106613013, .7035725550e-1]
RSS = 1.094542689
MSe = 0.3648
R^2 = .998010
对于后一组数据：
SSy = 0.0594
[M a K] = [18.68430621, .4725336711, .1934850297e-1]
RSS = .2883640565e-6
MSe = 1.4418e-007
R^2 = .999995

#8

wangweicupl2007-09-22 09:51

回复：（ydgsl）对于前一组数据，结果如下：SSy = 550...

不好意思啊，我刚才少给了一个数啊！

应该是x=[12.6743 12.7627 12.8453 12.8453 12.9227 12.9988 13.0756]
这是我的一次作业，你能给过程给我吗？
高手？

#9

wangweicupl2007-09-22 09:52

老师让我们一个月交，可是我最近还要花时间看四级的．所以这里不太懂写方程．

#10

wangweicupl2007-09-22 10:01

t=[1 2 3 4 5 6 ];
x=[12.6743 12.7627 12.8453 12.8453 12.9227 12.9988 13.0756]
求递增常数Ｋ和最大理论值Ｍ
知道你可能忙，我就整理了一下
辛苦你了啊！

#11

wangweicupl2007-09-22 13:29

回复：（wangweicupl）再加一句：求递增常数Ｋ和最大理...

t=[1 2 3 4 5 6 ];
x=[12.6743 12.7627 12.8453 12.8453 12.9227 12.9988 13.0756]
求递增常数Ｋ和最大理论值Ｍ
知道你可能忙，我就整理了一下
辛苦你了啊！
帮帮忙吧，这个很重要啊！

#12

ydgsl2007-09-22 14:53

[M a K] = [14.68771968, .1669652412, .5039047747e-1]
RSS = .1041610104e-4
MSe = 3.4720e-006
R^2 =.999907
上面的数据肯定还有错误，我删掉了一个12.8453。
又及，只取生长过程中间一段的数据来拟合是不应该的，不恰当的，其结果也基本无意义。

#13

wangweicupl2007-09-22 17:03

，非常感谢啊　．

#14

ydgsl2007-09-22 22:36

[M a K] =[14.68771952, .1669652294, .5039048117e-1]
RSS = .1041610104e-4
MSe = 3.4720e-006
R2 = 0.9999
（数据仍然有错，我删掉了一个12.8453）
由于上面所说的原因，拟合结果可能没有多大的实际意义。请参考使用。

#15

wangweicupl2007-09-23 01:45

,现在是半夜了啊！
呵呵，非常谢谢你啊！
我又去看一下原题：t=[1 2 3 4 5];
y=[12.7627 12.8453 12.9227 12.9988 13.0756].
你给的答案和参考答案很接近，你真的很有水平啊．
　　　你有空能加我吗？ＱＱ：５７４９４９０２２
　　　　还有可是给运算过程：用ＭＡＴＬＡＢ求的方程你能ＣＯＰＹ给我吗？

#16

wangweicupl2007-09-23 14:33

周末快过了啊！

#17

ydgsl2007-09-23 15:51

用lsqcurvefit 或 nlinfit 命令，以我给的参数值作为初值即可。

#18

wangweicupl2007-09-23 19:15

谢谢　　
有空能加我Ｑ吗？
以后还望多指教啊

#19

wangweicupl2007-09-23 19:16