实验二MATLAB中的数值计算与程序设计.doc

实验二 MATLAB中的数值计算与程序设计 一 实验目的： 完成MATLAB下数据分析与程序设计有关问题的练习。学习脚本文件的编写，数值输入，数组计算，程序调试等；练习利用MATLAB编程序解决基本数学问题。学习全局和局部变量的使用，流程控制结构，函数编写及其调试。 二 实验仪器和设备 装有MATLAB6.0以上计算机一台 三 实验原理 1. 常用运算命令 运算命令名 功能angle求复数的角min求最小值max求最大值sum求和roots求多项式的根poly由多项式的根求多项式的系数polyval求给定点多项式的值polyder多项式求导 例1 求a=S^3+2S^2+3S+4根. a=[1 2 3 4]; roots(a) ans = -1.6506 -0.1747 + 1.5469i -0.1747 - 1.5469i 或 roots([1 2 3 4]) ans = -1.6506 -0.1747 + 1.5469i -0.1747 - 1.5469i 例2 已知某多项式的根如下 -1.6506 -0.1747 + 1.5469i -0.1747 - 1.5469i 求该多项式. poly([-1.6506,-0.1747+1.5469i,-0.1747 - 1.5469i]) ans = 1.0000 2.0000 3.0001 4.0001 即所求多项式为 S^3+2S^2+3S+4 . 2. 无条件循环 当需要无条件重复执行某些命令时，可以使用for循环： for 循环变量t=表达式1 : 达式2 : 表达式3 语句体 end 说明：表达式1为循环初值，表达式2为步长，表达式3为循环终值；当表达式2省略时则默认步长为1；for语句允许嵌套。 例3： 例4：矩阵输入程序 生成3×4阶的Hiltber矩阵。 m=input(‘矩阵行数：m=’); for i=1 : 3 n= input(‘矩阵列数：n=’); for j=1 : 4 for i=1:m H（i，j）=1/（i+j-1）； for j=1:n end disp([‘输入第’,num2str(i),’行，第’,num2str(j),’ end 列元???’]) A(i, j) = input (‘ ’) end end 3. 条件循环 1) if-else-then语句 if-else-then语句的常使用三种形式为： (1) if 逻辑表达式 (3) if 逻辑表达式1 语句体 语句体1 end elseif 逻辑表达式2 语句体2 (2) if 逻辑表达式1 elseif 逻辑表达式3 语句体1 … else