数学建模Matlab实验程序

Matlab实验作业及答案

作业1：建立函数M文件

2.建立下面函数的M文件，并求f(x)值. (1)f(x)=log(x1+x22),x=(1,2) (2)f(x)=sin(x2)+exp(2x3),x=2 (1)

function f = fun1(x1,x2) f = log(x1 +x2^2); end (2)

function f = fun2( x) f = sin(x^2) + exp(2*x^3); end

3.试编写同时求sin(x),cos(x),exp(x),abs(x)的M函数文件.

function f = fun3(x)

f = [sin(x) cos(x) exp(x) abs(x)]; end

4.建立符号函数的M文件:当输入的变量为负数时，返回值-1；当输入的变量为正数时，返回

值1；而输入0时，返回值0.

function f = fun4(x) if x>0

f = 1; else if x == 0 f = 0; else f = -1; end end

x?1??e,x?0f??5.建立函数 的M文件。

2??x,x?0end

function f = fun5(x)

if x>0

f = exp(x-1); else

f = x^2; end end

6.通过帮助系统查询roots,poly,polyval,poly2str的用法，用这些命令解下面的问题：

已知一多项式的零点为{-1，1，2，3}，写出该多项式，并且计算多项式在点x=2.5处的值。

root = [-1 1 2 3]; p = poly(root); x = 2.5; a = polyval(p,x); eig(a)

计算多项式y=x3-3x+2的零点

P = [1 0 -3 2];a = company(p); eig(a)

7.查询sum,length的用法，建立一个求向量的平均值的M文件

a = [1 2 3 4]; b = sum(a);

c = length(a); d = b/c; eig(d)

8.查询input,disp的用法，建立M文件：输入x,y的值将其互换后输出

x = input('x='); y = input('y='); disp ([x,y]); t=x; x=y; y=t; disp ([x,y]);

作业题2：MATLAB矩阵的处理

1. 创建矩阵

A = [1 2 -1 3 5;1 -2 9 0 -6;-3 3 -4 7 1;9 8 0 7 6]; disp(A);

2.取A的1,2行与2,3列的交叉元素作子矩阵A1.

A([1,2],[2,3])

3.取A的1,3行，然后按行形成矩阵A2

A2 = A([1,3],:);

4.逆序提取A的1,2,3行，形成列矩阵A3.

A3 = [A(3,:) A(2,:) A(1,:)]’;(“’”为转置符号) 5.取A2的绝对值大于3的元素构成向量A4.

A4 = find(A>3);

6.求出A的最大值a及其所处的位置.

?1?1?A???3??92?238?19?4030775???6?1??6?a = max(max(A)); [row col v] = find(a); disp([row col v]); 7、设

?10?7?A??8??77565861097??5?9??10??32???23??b??33????31?用三种方法（克拉姆法则、矩阵的除法、逆矩阵）解方程组AX=b

A = [10 7 8 7;7 5 6 5;8 6 10 9;7 5 9 10]; b = [32;23;33;31];

逆矩阵法：

x = inv(A)*b; disp(x);

矩阵的除法

x =A\\b;

克拉姆法则

for n=1:4 B = A; B(:,n) = b;

x(n) = det(B)/det(A); end disp(x);

作业题3：Matlab语法控制结构的使用

(1) 用起泡法对10个数由小到大排序.即将相邻两个数比较,将小的调到前头.

a = [0 9 7 8 6 5 4 3 2 1];