回归分析：非线性nlinfi

今天测试。这首题，真的很头疼，第一次看到这个题，就知道要用nlinefit函数做，但是我一个地方没搞清楚，

花了我40多分钟还没做也来。

最后终于是调用的函数出错了。主要是没有将一个列抽出来。一个二维的当一维做，能对吗？只能呵呵了！！！

在看题之前，先熟悉数组的操作：

创建数组

x=[a b c d e f ] 创建包含指定元素的行向量

x=first：last 创建从first开始，加1计数，到last结束的行向量

x=first:increment:last 创建从first开始，加increment计数，last结束的行向量

x=linspace(first，last，n） 创建从first开始，到last结束，有n个元素的行向量

x=logspace(first，last，n） 创建从开始，到结束，有n个元素的对数分隔行向量

数组元素的访问

访问一个元素： x(i)表示访问数组x的第i个元素.

访问一块元素： x(a:b:c)表示访问数组x的从第a个元素开始，以步长为b到c个元素（但不超过c），b可以为负数，b缺损时为1.

直接使用元素编址序号. x([a b c d]) 表示提取数组x的第a、b、c、d个元素构成一个新的数组

[x(a) x(b) x(c) x(d)].这点很重要的！！！！！！！！！

还有其它的操作，有兴趣的还可以找百度文库中的PPT:如下

在研究化学动力学反应过程中，建立了一个反应速度和反应物的数学模型,形式为

其中 为未知参数，是三种反应物（氢，n戊烷，异构戊烷）的含量，y是反应速度今测得一组数据如下表，试由此确定参数 。并给出置信区间。 的参考值（1，0.05，0.02，0.1，2）

序号反应速度y 氢 n戊烷 异构戊烷

1     8.55 470 300 10

2     3.79 285 80 10

3     4.82 470 300 120

4     0.02 470 80 120

5     2.75 470 80 10

6     14.39 100 190 10

7     2.54 100 80 65

8     4.35 470 190 65

9     13.00 100 300 54

10      8.5 100 300 120

11     0.05 100 80 120

12     11.32 285 300 10

13      3.13 285 190 120

用到的函数如下;

x=[470 300 10

285 80 10

470 300 120

470 80 120

470 80 10

100 190 10

100 80 65

470 190 65

100 300 54

100 300 120

100 80 120

185 190 10

285 190 120];

y=[8.55

3.79

4.82

0.02

2.75

14.39

2.54

4.35

13.00

8.5

0.05

11.32

3.13];

beta0=[1 0.05 0.02 0.1 2];

[beta,r,J]=nlinfit(x,y,'model',beta0);

[Y,DELTA]=nlpredci('model',x,beta,r,J);

Y+DELTA

Y-DELTA

 

 

 

二、 （6分）观察落叶松的树龄x和平均高度H有如下资料(见表14-4)

表14-4

xi	2	3	4	5	6	7	8	9	10	11
hi	5.6	8	10.4	12.8	15.3	17.8	19.9	21.4	22.4	23.2

若h对x的回归方程为抛物线型，试求出气方程中的未知参数，并作回归分析.

使用到的函数是polyfit和polyconf这两个函数。

 

这题读者可以试着求解！