Matlab与相关工具箱解说.ppt

例1 解： 1、输入数据： x=[143 145 146 147 149 150 153 154 155 156 157 158 159 160 162 164]; X=[ones(16,1) x]; Y=[88 85 88 91 92 93 93 95 96 98 97 96 98 99 100 102]; 2、回归分析及检验： [b,bint,r,rint,stats]=regress(Y,X) b,bint,stats To MATLAB(liti11) 1、考察温度x对产量y的影响，测得下列10组数据： 求y关于x的线性回归方程，检验回归效果是否显著，并预测x=42℃时产量的估值及预测区间（置信度95%）. 2、某零件上有一段曲线，为了在程序控制机床上加工这一零件，需要求这段曲线的解析表达式，在曲线横坐标xi处测得纵坐标yi共11对数据如下： 求这段曲线的纵坐标y关于横坐标x的二次多项式回归方程. 例3-9 已知 ，当n=100时，求y的值。 程序如下： y=0; n=100; for i=1:n y=y+1/(2*i-1); end y 在实际MATLAB编程中，采用循环语句会降低其执行速度，所以前面的程序通常由下面的程序来代替： n=100; i=1:2:2*n-1; y=sum(1./i); y for语句更一般的格式为： for 循环变量=矩阵表达式 循环体语句 end 执行过程是依次将矩阵的各列元素赋给循环变量，然后执行循环体语句，直至各列元素处理完毕。 例3-10 写出下列程序的执行结果。 s=0; a=[12,13,14;15,16,17;18,19,20;21,22,23]; for k=a s=s+k; end disp(s); 2．while语句 while语句的一般格式为： while (条件) 循环体语句 end 其执行过程为：若条件成立，则执行循环体语句，执行后再判断条件是否成立，如果不成立则跳出循环。 例3-11 从键盘输入若干个数，当输入0时结束输入，求这些数的平均值和它们之和。 程序如下： sum=0; cnt=0; val=input(Enter a number (end in 0):); while (val~=0) sum=sum+val; cnt=cnt+1; val=input(Enter a number (end in 0):); end if (cnt 0) sum mean=sum/cnt end 3．break语句和continue语句 与循环结构相关的语句还有break语句和continue语句。它们一般与if语句配合使用。 break语句用于终止循环的执行。当在循环体内执行到该语句时，程序将跳出循环，继续执行循环语句的下一语句。 continue语句控制跳过循环体中的某些语句。当在循环体内执行到该语句时，程序将跳过循环体中所有剩下的语句，继续下一次循环。 例3-12 求[100，200]之间第一个能被21整除的整数。 程序如下： for n=100:200 if rem(n,21)~=0 continue end break end n 4．循环的嵌套 如果一个循环结构的循环体又包括一个循环结构，就称为循环的嵌套，或称为多重循环结构。 例3-13 若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6是完数。求[1,500]之间的全部完数。 for m=1:500 s=0; for k=1:m/2 if rem(m,k)==0 s=s+k; end end if m==s disp(m); end end 3.3 函数文件 3.3.1 函数文件的基本结构 函数文件由function语句引导，其基本结构为： function 输出形参表=函数名(输入形参表) 注释说明部分 函数体语句 其中以function开头的一行为引导行，表示该M文件是一个函数文件。函数名的命名规则与变量名相同。输入形参为函数的输入参数，输出形参为函数的输出参数。当输出形参多于一个时，则应该用方括号括起来。 例3-14 编写函数文件求半径为r的圆的面积和周长。 函数文件如下： function [s,p]=fcircle(r) %CIR