matlab程序设计实践

⑵ 例题

建立myf.m函数文件 建立dmyf.m导数的雅克比矩阵 [r,n]=mulDNewton([0,0]) 输出计算结果 输出迭代次数

2．编程解决以下科学计算和工程实际问题 ⑴

2-4在平炉炼钢过程中，由于矿石及炉气的氧化作用，铁水中的总含碳量在不断降低，一炉钢在冶炼初期（熔化期）中总的去碳量y与所加的两种矿石（天然矿石与烧结矿石）的量x1,x2及熔化时间x3有关（熔化时间愈长则去碳量愈多）。经实测，某号平炉的49炉数据如表2.6所示，试求出y与x1,x2,x3之间的线性回归方程。

解：

在MATLAB命令窗口中输入：

>> x1 =[2 7 5 12 1 3 2 6 7 0 3 0 8 6 0 3 7 16 6 0 9 4 0 9 2 9 12 6 12 0 5 4 0 6 4 10 4 5 9 6 5 5 8 2 7 4 10 3 4];

>> x2 =[18 9 14 3 20 12 17 5 8 23 16 18 4 14 21 14 12 0 16 15 0 6 17 0 16 6 5 13 7 24 12 15 20 16 17 4 14 13 8 13 8 11 6 13 8 10 5 17 15];

>> x3 =[50 40 46 43 64 40 64 39 37 55 60 49 50 51 51 51 56 48 46 52 40 32 47 44 39 39 61 41 47 61 37 48 45 42 48 48 36 36 51 54 100 44 63 55 50 45 40 64 72];

>> y=[4.3302 3.6485 4.4830 5.5468 5.4970 3.1125 5.1182 3.8759 4.6700 4.9536 5.0060 5.2701 5.8772 5.4849 4.5960 5.6645 6.0795 3.2194 5.8076 4.7306 4.6805 3.1272 2.6104 3.7174 3.8946 2.7066 5.6314 5.8152 5.1302 5.3910 4.4533 4.6569 4.5212 4.8650 5.3566 4.6096 2.3815 3.8746 4.5919 5.1588 5.4373 3.9960 4.3970 4.0622 2.2905 4.7115 4.5310 5.3637 6.0771]; >> y=y';

>> x=[x1',x2',x3',ones(49,1)]; >> regress(y,x) ans =

0.1504 0.1014 0.0370 0.7753 得到：

0.1504y=0.1014x1+0.0370x2+0.7753x3 流程图：

开始 输入y、x1、x2、x3数据 y=y' x=[x1',x2',x3',ones(49,1)]; regress(y,x) 输出y、x1、x2、x3的值，得到线性回归方程 结束