MATLAB数学实验第二版答案

数学实验答案

Chapter 1

Page20,ex1

(5) 等于[exp(1),exp(2);exp(3),exp(4)]

(7) 3=1*3, 8=2*4

(8) a为各列最小值，b为最小值所在的行号

(10) 1>=4,false, 2>=3,false, 3>=2, ture, 4>=1,ture

(11) 答案表明：编址第2元素满足不等式(30>=20)和编址第4元素满足不等式(40>=10)

(12) 答案表明：编址第2行第1列元素满足不等式(30>=20)和编址第2行第2列元素满足不等式(40>=10)

Page20, ex2

(1)a, b, c的值尽管都是1，但数据类型分别为数值，字符，逻辑，注意a与c相等，但他们不等于b

(2)double(fun)输出的分别是字符a,b,s,(,x,)的ASCII码

Page20,ex3

>> r=2;p=;n=12;

>> T=log(r)/n/log(1+*p)

Page20,ex4

>> x=-2::2;f=x.^4-2.^x;

>> [fmin,min_index]=min(f)

最小值最小值点编址

>> x(min_index)

ans =

最小值点

>> [f1,x1_index]=min(abs(f)) 求近似根--绝对值最小的点

f1 =

x1_index =

24

>> x(x1_index)

ans =

>> x(x1_index)=[];f=x.^4-2.^x; 删去绝对值最小的点以求函数绝对值次小的点

>> [f2,x2_index]=min(abs(f)) 求另一近似根--函数绝对值次小的点

f2 =

x2_index =

65

>> x(x2_index)

ans =

>> z=magic(10)

z =

92 99 1 8 15 67 74 51 58 40

98 80 7 14 16 73 55 57 64 41

4 81 88 20 22 54 56 63 70 47

85 87 19 21 3 60 62 69 71 28

86 93 25 2 9 61 68 75 52 34

17 24 76 83 90 42 49 26 33 65

23 5 82 89 91 48 30 32 39 66

79 6 13 95 97 29 31 38 45 72

10 12 94 96 78 35 37 44 46 53

11 18 100 77 84 36 43 50 27 59

>> sum(z)

>> sum(diag(z))

>> z(:,2)/sqrt(3)

>> z(8,:)=z(8,:)+z(3,:)

Chapter 2

Page 45 ex1

先在编辑器窗口写下列M函数，保存为function [xbar,s]=ex2_1(x)

n=length(x);

xbar=sum(x)/n;

s=sqrt((sum(x.^2)-n*xbar^2)/(n-1));

例如

>>x=[81 70 65 51 76 66 90 87 61 77]; >>[xbar,s]=ex2_1(x)

Page 45 ex2

s=log(1);n=0;

while s<=100

n=n+1;

s=s+log(1+n);

end

m=n

Page 40 ex3

clear;

F(1)=1;F(2)=1;k=2;x=0;

e=1e-8; a=(1+sqrt(5))/2;

while abs(x-a)>e

k=k+1;F(k)=F(k-1)+F(k-2); x=F(k)/F(k-1); end

a,x,k

计算至k=21可满足精度

clear;tic;s=0;

for i=1:1000000

s=s+sqrt(3)/2^i;

end

s,toc

tic;s=0;i=1;

while i<=1000000

s=s+sqrt(3)/2^i;i=i+1;

end

s,toc

tic;s=0;

i=1:1000000;

s=sqrt(3)*sum(1./2.^i);

s,toc

Page 45 ex5

t=0:24;

c=[15 14 14 14 14 15 16 18 20 22 23 25 28 ...

31 32 31 29 27 25 24 22 20 18 17 16];

plot(t,c)

Page 45 ex6

(1)

x=-2::2;y=x.^2.*sin(x.^2-x-2);plot(x,y)

y=inline('x^2*sin(x^2-x-2)');fplot(y,[-2 2])

(2)参数方法

t=linspace(0,2*pi,100);

x=2*cos(t);y=3*sin(t); plot(x,y)

(3)

x=-3::3;y=x;

[x,y]=meshgrid(x,y);

z=x.^2+y.^2;

surf(x,y,z)

(4)

x=-3::3;y=-3::13;

[x,y]=meshgrid(x,y);

z=x.^4+3*x.^2+y.^2-2*x-2*y-2*x.^2.*y+6;

surf(x,y,z)

(5)

t=0::2*pi;

x=sin(t);y=cos(t);z=cos(2*t);

plot3(x,y,z)

(6)

theta=linspace(0,2*pi,50);fai=linspace(0,pi/2,20); [theta,fai]=meshgrid(theta,fai);

x=2*sin(fai).*cos(theta);

y=2*sin(fai).*sin(theta);z=2*cos(fai);

surf(x,y,z)

(7)

x=linspace(0,pi,100);

y1=sin(x);y2=sin(x).*sin(10*x);y3=-sin(x);

plot(x,y1,x,y2,x,y3)

page45, ex7

x=::;

y=*(x>+x.*(x<=.*(x>=*(x<;

plot(x,y)

page45,ex9

clear;close;

x=-2::2;y=x;

[x,y]=meshgrid(x,y);

a=;b=;

p=a*exp*y.^*x.^*x).*(x+y>1);

p=p+b*exp(-y.^2-6*x.^2).*(x+y>-1).*(x+y<=1);

p=p+a*exp*y.^*x.^2+*x).*(x+y<=-1);

mesh(x,y,p)

page45, ex10

lookfor lyapunov

help lyap

>> A=[1 2 3;4 5 6;7 8 0];C=[2 -5 -22;-5 -24 -56;-22 -56 -16]; >> X=lyap(A,C)

X =

Chapter 3

Page65 Ex1

>> a=[1,2,3];b=[2,4,3];a./b,a.\b,a/b,a\b

ans =

ans =

2 2 1

ans =

一元方程组x[2,4,3]=[1,2,3]的近似解

ans =

0 0 0

0 0 0

矩阵方程[1,2,3][x11,x12,x13;x21,x22,x23;x31,x32,x33]=[2,4,3]的特解

Page65 Ex 2

(1)

>> A=[4 1 -1;3 2 -6;1 -5 3];b=[9;-2;1];

>> rank(A), rank([A,b]) [A,b]为增广矩阵

ans =

3

ans =

3 可见方程组唯一解

>> x=A\b

x =

(2)

>> A=[4 -3 3;3 2 -6;1 -5 3];b=[-1;-2;1];

>> rank(A), rank([A,b])

ans =

3

ans =

3 可见方程组唯一解

>> x=A\b

x =

(3)

>> A=[4 1;3 2;1 -5];b=[1;1;1];

>> rank(A), rank([A,b])

ans =

2

ans =

3 可见方程组无解

>> x=A\b

x =

最小二乘近似解

(4)

>> a=[2,1,-1,1;1,2,1,-1;1,1,2,1];b=[1 2 3]';%注意b的写法>> rank(a),rank([a,b])

ans =

3

ans =

3 rank(a)==rank([a,b])<4说明有无穷多解

>> a\b

ans =

1

1

0 一个特解

Page65 Ex3

>> a=[2,1,-1,1;1,2,1,-1;1,1,2,1];b=[1,2,3]'; >> x=null(a),x0=a\b

x =

x0 =

1

1

通解kx+x0

Page65 Ex 4

>> x0=[ ]';a=[ ; ];

>> x1=a*x, x2=a^2*x, x10=a^10*x

>> x=x0;for i=1:1000,x=a*x;end,x

x =

>> x0=[ ]';

>> x=x0;for i=1:1000,x=a*x;end,x

x =

>> [v,e]=eig(a)

v =

e =

>> v(:,1)./x

ans =

成比例，说明x是最大特征值对应的特征向量Page65 Ex5

用到公式

>> B=[6,2,1;,1,;3,,];x=[25 5 20]';

>> C=B/diag(x)

C =

>> A=eye(3,3)-C

A =

>> D=[17 17 17]';x=A\D

x =

Page65 Ex 6

(1)

>> a=[4 1 -1;3 2 -6;1 -5 3];det(a),inv(a),[v,d]=eig(a) ans =

-94

ans =

v =

d =

0 0

0 0

0 0

(2)

>> a=[1 1 -1;0 2 -1;-1 2 0];det(a),inv(a),[v,d]=eig(a) ans =

1

ans =

v =

+ -

- +

d =

0 0

0 + 0

0 0 -

(3)

>> A=[5 7 6 5;7 10 8 7;6 8 10 9;5 7 9 10]

A =

5 7

6 5

7 10 8 7

6 8 10 9

5 7 9 10

>> det(A),inv(A), [v,d]=eig(A)

ans =

1

ans =

v =

d =

0 0 0

0 0 0

0 0 0

0 0 0

(4)(以n=5为例)

方法一（三个for）

n=5;

for i=1:n, a(i,i)=5;end

for i=1:(n-1),a(i,i+1)=6;end

for i=1:(n-1),a(i+1,i)=1;end

a

方法二（一个for）

n=5;a=zeros(n,n);

a(1,1:2)=[5 6];

for i=2:(n-1),a(i,[i-1,i,i+1])=[1 5 6];end a(n,[n-1 n])=[1 5];

a

方法三（不用for）

n=5;a=diag(5*ones(n,1));

b=diag(6*ones(n-1,1));

c=diag(ones(n-1,1));

a=a+[zeros(n-1,1),b;zeros(1,n)]+[zeros(1,n);c,zeros(n-1,1)] 下列计算

>> det(a)

ans =

665

>> inv(a)

ans =

>> [v,d]=eig(a)

v =

d =

0 0 0 0

0 0 0 0

0 0 0 0

0 0 0 0

0 0 0 0

Page65 Ex 7

(1)

>> a=[4 1 -1;3 2 -6;1 -5 3];[v,d]=eig(a)

v =

d =

0 0

0 0

0 0

>> det(v)

ans =

%v行列式正常, 特征向量线性相关，可对角化

>> inv(v)*a*v 验算

ans =

>> [v2,d2]=jordan(a) 也可用jordan

v2 =

特征向量不同

d2 =

0 0

0 - 0

0 0 +

>> v2\a*v2

ans =

>> v(:,1)./v2(:,2) 对应相同特征值的特征向量成比例ans =

(2)

>> a=[1 1 -1;0 2 -1;-1 2 0];[v,d]=eig(a)

v =

+ -

- +

d =

0 0

0 + 0

0 0 -

>> det(v)

ans =

v的行列式接近0, 特征向量线性相关，不可对角化>> [v,d]=jordan(a)

v =

1 0 1

1 0 0

1 -1 0

d =

1 1 0

0 1 1

0 0 1 jordan标准形不是对角的，所以不可对角化

>> A=[5 7 6 5;7 10 8 7;6 8 10 9;5 7 9 10]

A =

5 7

6 5

7 10 8 7

6 8 10 9

5 7 9 10

>> [v,d]=eig(A)

v =

d =

0 0 0

0 0 0

0 0 0

0 0 0

>> inv(v)*A*v

ans =

本题用jordan不行, 原因未知

(4)

参考6(4)和7(1)

Page65 Exercise 8

只有(3)对称, 且特征值全部大于零, 所以是正定矩阵. Page65 Exercise 9

(1)

>> a=[4 -3 1 3;2 -1 3 5;1 -1 -1 -1;3 -2 3 4;7 -6 -7 0] >> rank(a)

ans =

3

>> rank(a(1:3,:))

ans =

2

>> rank(a([1 2 4],:)) 1,2,4行为最大无关组

ans =

3

>> b=a([1 2 4],:)';c=a([3 5],:)';

>> b\c 线性表示的系数

ans =

Page65 Exercise 10

>> a=[1 -2 2;-2 -2 4;2 4 -2]

>> [v,d]=eig(a)

v =

d =

0 0

0 0

0 0

>> v'*v

ans =

0 v确实是正交矩阵

Page65 Exercise 11

设经过6个电阻的电流分别为i1, ..., i6. 列方程组如下

20-2i1=a; 5-3i2=c; a-3i3=c; a-4i4=b; c-5i5=b; b-3i6=0;

i1=i3+i4;i5=i2+i3;i6=i4+i5;

计算如下

>> A=[1 0 0 2 0 0 0 0 0;0 0 1 0 3 0 0 0 0;1 0 -1 0 0 -3 0 0 0; 1 -1 0 0 0 0 -4 0 0;

0 -1 1 0 0 0 0 -5 0;0 1 0 0 0 0 0 0 -3; 0 0 0 1 0 -1 -1 0 0;0 0 0 0 -1 -1 0 1 0;

0 0 0 0 0 0 -1 -1 1];

>>b=[20 5 0 0 0 0 0 0 0]'; A\b

ans =

Page65 Exercise 12

>> A=[1 2 3;4 5 6;7 8 0];

>> left=sum(eig(A)), right=sum(trace(A))

left =

6

>> left=prod(eig(A)), right=det(A) 原题有错, (-1)^n应删去left =

right =

27

>> fA=(A-p(1)*eye(3,3))*(A-p(2)*eye(3,3))*(A-p(3)*eye(3,3)) fA =

*

>> norm(fA) f(A)范数接近0

ans =

Chapter 4

Page84 Exercise 1

(1)

roots([1 1 1])

(2)

roots([3 0 -4 0 2 -1])

(3)

p=zeros(1,24);

p([1 17 18 22])=[5 -6 8 -5];

roots(p)

(4)

p1=[2 3];

p2=conv(p1, p1);

p3=conv(p1, p2);

p3(end)=p3(end)-4; %原p3最后一个分量-4

roots(p3)

Page84 Exercise 2

fun=inline('x*log(sqrt(x^2-1)+x)-sqrt(x^2-1)*x');

fzero(fun,2)

Page84 Exercise 3

fun=inline('x^4-2^x');

fplot(fun,[-2 2]);grid on;

fzero(fun,-1),fzero(fun,1),fminbnd(fun,,

Page84 Exercise 4

fun=inline('x*sin(1/x)','x');

fplot(fun, [ ]);

x=zeros(1,10);for i=1:10, x(i)=fzero(fun,*;end;

Page84 Exercise 5

fun=inline('[9*x(1)^2+36*x(2)^2+4*x(3)^2-36;x(1)^2-2*x(2)^2-20*x(3);16*x(1)-x(1 )^3-2*x(2)^2-16*x(3)^2]','x');

[a,b,c]=fsolve(fun,[0 0 0])

Page84 Exercise 6

fun=@(x)[x(1)*sin(x(1))*cos(x(2)),x(2)*cos(x(1))+*sin(x(2))];

[a,b,c]=fsolve(fun,[ ])

Page84 Exercise 7

clear; close; t=0:pi/100:2*pi;

x1=2+sqrt(5)*cos(t); y1=3-2*x1+sqrt(5)*sin(t);

x2=3+sqrt(2)*cos(t); y2=6*sin(t);

plot(x1,y1,x2,y2); grid on; 作图发现4个解的大致位置，然后分别求解

y1=fsolve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[,2])

y2=fsolve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[,-2]) y3=fsolve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[,-5]) y4=fsolve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[4,-4]) Page84 Exercise 8

(1)

clear;

fun=inline('x.^2.*sin(x.^2-x-2)');

fplot(fun,[-2 2]);grid on; 作图观察

x(1)=-2;

x(3)=fminbnd(fun,-1,;

x(5)=fminbnd(fun,1,2);

fun2=inline('-x.^2.*sin(x.^2-x-2)');

x(2)=fminbnd(fun2,-2,-1);

x(4)=fminbnd(fun2,,;

x(6)=2

feval(fun,x)

答案: 以上x(1)(3)(5)是局部极小，x(2)(4)(6)是局部极大，从最后一句知道x(1)全局最小， x(2)最大。

(2)

clear;

fun=inline('3*x.^5-20*x.^3+10');

fplot(fun,[-3 3]);grid on; 作图观察

x(1)=-3;

x(3)=fminsearch(fun,;

fun2=inline('-(3*x.^5-20*x.^3+10)');

x(2)=fminsearch(fun2,;

x(4)=3;

feval(fun,x)

(3)

fun=inline('abs(x^3-x^2-x-2)');

fplot(fun,[0 3]);grid on; 作图观察

fminbnd(fun,,

fun2=inline('-abs(x^3-x^2-x-2)');

fminbnd(fun2,,

Page84 Exercise 9

close;

x=-2::1;y=-7::1;

[x,y]=meshgrid(x,y);

z=y.^3/9+3*x.^2.*y+9*x.^2+y.^2+x.*y+9;

mesh(x,y,z);grid on; 作图观察

fun=inline('x(2)^3/9+3*x(1)^2*x(2)+9*x(1)^2+x(2)^2+x(1)*x(2)+9');

x=fminsearch(fun,[0 0]) 求极小值

fun2=inline('-(x(2)^3/9+3*x(1)^2*x(2)+9*x(1)^2+x(2)^2+x(1)*x(2)+9)'); x=fminsearch(fun2,[0 -5]) 求极大值

Page84 Exercise 10

clear;t=0:24;

c=[15 14 14 14 14 15 16 18 20 22 23 25 28 ...

31 32 31 29 27 25 24 22 20 18 17 16];

p2=polyfit(t,c,2)

p3=polyfit(t,c,3)

fun=inline('a(1)*exp(a(2)*(t-14).^2)','a','t');

a=lsqcurvefit(fun,[0 0],t,c) 初值可以试探

f=feval(fun, a,t)

norm(f-c) 拟合效果

plot(t,c,t,f) 作图检验

fun2=inline('b(1)*sin(pi/12*t+b(2))+20','b','t'); 原题修改f(x)+20 b=lsqcurvefit(fun2,[0 0],t,c)

figure

f2=feval(fun2, b,t)

norm(f2-c) 拟合效果

plot(t,c,t,f2) 作图检验

Page84 Exercise 11

fun=inline('(1-x)*sqrt+x)*x*sqrt(1+x)*sqrt(5)');

x=fzero(fun, 0, 1)

Page84 Exercise 12

r=12/100;N=20*12;

x=7500*180 房屋总价格

y=x* 首付款额

x0=x-y 贷款总额

a=(1+r)^N*r*x0/((1+r)^N-1) 月付还款额

r1=12/100;x1=10*10000; 公积金贷款

a1=(1+r1)^N*r1*x1/((1+r1)^N-1)

x2=x0-x1 商业贷款

a2=(1+r)^N*r*x2/((1+r)^N-1)

a=a1+a2

Page84 Exercise 13

列方程th*R^2+(pi-2*th)*r^2-R*r*sin(th)=pi*r^2/2

化简得sin(2*th)-2*th*cos(2*th)=pi/2

以下Matlab计算

clear;fun= inline('sin(2*th)-2*th*cos(2*th)-pi/2','th')

th=fsolve(fun,pi/4)

R=20*cos(th)

Page84 Exercise 14

先在Editor窗口写M函数保存

function x=secant(fname,x0,x1,e)

while abs(x0-x1)>e,

x=x1-(x1-x0)*feval(fname,x1)/(feval(fname,x1)-feval(fname,x0)); x0=x1;x1=x;

end

再在指令窗口

fun=inline('x*log(sqrt(x^2-1)+x)-sqrt(x^2-1)*x');

secant(fun,1,2,1e-8)

Page84 Exercise 15

作系数为a,初值为xo,从第m步到第n步迭代过程的M函数：

function f=ex4_15fun(a,x0,m,n)

x(1)=x0; y(1)=a*x(1)+1;x(2)=y(1);

if m<2, plot([x(1),x(1),x(2)],[0,y(1),y(1)]);hold on; end

for i=2:n

y(i)=a*x(i)+1; x(i+1)=y(i);

if i>m, plot([x(i),x(i),x(i+1)],[y(i-1),y(i),y(i)]); end

end

hold off;

M脚本文件

subplot(2,2,1);ex4_15fun,1,1,20);

subplot(2,2,2);ex4_15fun,1,1,20);

subplot(2,2,3);ex4_15fun,1,1,20);

subplot(2,2,4);ex4_15fun,1,1,20);

Page84 Exercise 16

设夹角t, 问题转化为 min f=5/sin(t)+10/cos(t)

取初始值pi/4, 计算如下

fun=@(t)5/sin(t)+10/cos(t);

[t,f]=fminsearch(fun, pi/4)

t =

f =

Page84 Exercise 17

提示：x(k+2)=f(x(k))=a^2*x(k)*(1-x(k))*(1-a*x(k)*(1-x(k)))

计算平衡点x

|f'(x)|<1则稳定

Page84 Exercise 18

先写M文件

function f=ex4_18(a,x0,n)

x=zeros(1,n);y=x;

x(1)=x0;

y(1)=a*x(1)+1;

x(2)=y(1);

plot([x(1),x(1),x(2)],[0,y(1),y(1)],'r'); hold on;

for i=2:n

y(i)=a*x(i)+1;

x(i+1)=y(i);

plot([x(i),x(i),x(i+1)],[y(i-1),y(i),y(i)]) end

hold off;

再执行指令

>> ex4_18,1,20)

>> ex4_18,1,20)

>> ex4_18,1,20)

>> ex4_18,1,20)

Page84 Exercise 19

clear; close; x(1)=0; y(1)=0;

for k=1:3000

x(k+1)=1+y(k)*x(k)^2; y(k+1)=*x(k);

end

plot(x(1000:1500),y(1000:1500),'+g');hold on plot(x(1501:2000),y(1501:2000),'.b');

plot(x(2001:2500),y(2001:2500),'*y');

plot(x(2501:3001),y(2501:3001),'.r'); Chapter 5

Page101 Exercise 1

x=[0 4 10 12 15 22 28 34 40];

y=[0 1 3 6 8 9 5 3 0];

trapz(x,y)

Page101 Exercise 2

x=[0 4 10 12 15 22 28 34 40];

y=[0 1 3 6 8 9 5 3 0];

diff(y)./diff(x)

Page101 Exercise 3

xa=-1::1;ya=0::2;

[x,y]=meshgrid(xa,ya);

z=x.*exp(-x.^2 -y.^3);

[px,py] = gradient(z,xa,ya);

px

Page101 Exercise 4

t=0::;

x=log(cos(t));

y=cos(t)-t.*sin(t);

dydx=gradient(y,x)

[x_1,id]=min(abs(x-(-1)));%找最接近x=-1的点

dydx(id)

Page101 Exercise 5

(1)

Fun=inline(‘1/(sqrt(2*pi)).*exp(-x.^2./2)’);

Quadl(fun,0,1)

(2)

fun=inline('exp(2*x).*cos(x).^3');

quadl(fun,0,2*pi)

或用trapz

x=linspace(0,2*pi,100);

y=exp(2*x).*cos(x).^3;

trapz(x,y)

(3)

fun=@(x)x.*log(x.^4).*asin(1./x.^2);

quadl(fun,1,3)

或用trapz

x=1::3;

y=feval(fun,x);

trapz(x,y)

(4)

fun=@(x)sin(x)./x;

quadl(fun,1e-10,1) %注意由于下限为0，被积函数没有意义，用很小的1e-10代替 (5)

%参考Exercise 5(4)

(6)

fun=inline('sqrt(1+r.^2.*sin(th))','r','th');

dblquad(fun,0,1,0,2*pi)

(7)

首先建立84页函数dblquad2

clear;

fun=@(x,y)1+x+y.^2;

clo=@(x)-sqrt(2*x-x.^2);

dup=@(x)sqrt(2*x-x.^2);

dblquad2(fun,0,2,clo,dhi,100)

Page101 Exercise 6

t=linspace(0,2*pi,100);

x=2*cos(t);y=3*sin(t);

dx=gradient(x,t);dy=gradient(y,t);

f=sqrt(dx.^2+dy.^2);

trapz(t,f)

Page101 Exercise 7

xa=-1::1;ya=0::2;

[x,y]=meshgrid(xa,ya);

z=x.*exp(x.^2+y.^2);

[zx,zy]=gradient(z,xa,ya);

f=sqrt(1+zx.^2+zy.^2);

s=0;

for i=2:length(xa)

for j=2:length(ya)

s=s+(xa(i)-xa(i-1))*(ya(j)-ya(j-1))*(f(i,j)+f(i-1,j)+f(i,j-1)+f(i-1,j-1))/4; end

end

s

Page101 Exercise 8

funl=inline('-(-x).^.*cos(x)');

funr=inline('x.^.*cos(x)');

quadl(funl,-1,0)+quadl(funr,0,1)

Page101 Exercise 9 (以I32为例)

fun=@(x)abs(sin(x));

h=;x=0:h:32*pi;y=feval(fun,x);t1=trapz(x,y)

h=pi;x=0:h:32*pi;y=feval(fun,x);t2=trapz(x,y)%步长与周期一致，结果失真

q1=quad(fun,0,32*pi)

q2=quadl(fun,0,32*pi)

Page101 Exercise 10

(2)

先在程序编辑器，写下列函数，保存为ex5_10_2f

function d=ex5_10_2f(fname,a,h0,e)

h=h0;d=(feval(fname,a+h)-2*feval(fname,a)+feval(fname,a-h))/(h*h);

d0=d+2*e;

while abs(d-d0)>e

d0=d;h0=h;h=h0/2;

d=(feval(fname,a+h)-2*feval(fname,a)+feval(fname,a-h))/(h*h);

end

再在指令窗口执行

fun=inline('x.^2*sin(x.^2-x-2)','x');

d=ex5_10_2f(fun,,,1e-3)

Page101 Exercise 11

提示：f上升时，f'>0;f下降时，f'<0; f极值， f'=0.

Page101 Exercise 12

在程序编辑器，写下列函数，保存为ex5_12f

function I=ex5_12(fname,a,b,n)

x=linspace(a,b,n+1);

y=feval(fname,x);

I=(b-a)/n/3*(y(1)+y(n+1)+2*sum(y(3:2:n))+4*sum(y(2:2:n)));

再在指令窗口执行

ex5_12(inline('1/sqrt(2*pi)*exp(-x.^2/2)'),0,1,50)

Page101 Exercise 13

fun=inline('5400*v./*v.^2+2000)','v');

quadl(fun,15,30)

Page101 Exercise 14

重心不超过凳边沿。1/2, 2/3, 3/4, ...,n/(n+1)

Page101 Exercise15

利润函数fun=inline('(p-c0+k*log(M*exp(-a*p)))*M*exp(-a*p)','p');

求p使fun最大

Page101 Exercise 16

clear; x=-3/4::3/4;

y=(3/4+x)*2.*sqrt(1-16/9.*x.^2)*;

P=trapz(x,y) %单位：千牛

Page101 Exercise 17

clear; close;

fplot('17-t^(2/3)-5-2*t^(2/3)',[0,20]); grid;

t=fzero('17-x^(2/3)-5-2*x^(2/3)',7)

t=0::8; y=17-t.^(2/3)-5-2*t.^(2/3);

trapz(t,y)-20 %单位：百万元

Page101 Exercise 18

曲面面积计算

Chapter 6

Page121 Exercise 1

(1)

fun=inline('x+y','x','y');

[t,y]=ode45(fun,[0 1 2 3],1) %注意由于初值为y(0)=1,[0 1 2 3]中0不可缺 (3)

令y(1)=y,y(2)=y',化为方程组

y(1)'=y(2),y(2)'=*y(2)^2-2*y(1)+sin(t)

运行下列指令

clear;close;

fun=@(t,y)[y(2);*y(2)^2-2*y(1)+sin(t)];

[t,y]=ode45(fun,[0 5],[0;1]);

plot(t,y(:,1))

(5)

令y(1)=y,y(2)=y',化为方程组

y(1)'=y(2),y(2)'=-mu*(y(1)^2-1)*y(2)-y(1)

Matlab数学实验报告一

数学软件课程设计 题目非线性方程求解 班级数学081 姓名曹曼伦

实验目的：用二分法与Newton迭代法求解非线性方程的根； 用Matlab函数solve、fzero、fsolve求解非线性方程（组）的解。 编程实现二分法及Newton迭代法； 学会使用Matlab函数solve、fzero、fsolve求解非线性方程（组）的解。 通过实例分别用二分法及迭代法解非线性方程组并观察收敛速度。 实验内容： 比较求exp(x)+10*x-2的根的计算量。(要求误差不超过十的五次方） (1)在区间（0，1)内用二分法； (2)用迭代法x=(2-exp(x))/10，取初值x=0 。 试验程序 (1)二分法: format long syms x s=exp(x)+10*x-2 a=0; b=1; A=subs(s,a) B=subs(s,b) f=A*B %若f<0，则为由根区间 n=0; stop=1.0e-5; while f<0&abs(a-b)>=stop&n<=100; Xk=(a+b)/2; %二分 M= subs(s, Xk); if M* A<0 symbol=1 %若M= subs(s, Xk)为正，则与a二分 b= Xk else symbol=0 % 若M= subs(s, Xk)为负，则与b二分 a= Xk end n=n+1 end Xk n (2)牛顿迭代法； format long

syms x s= (2-exp(x))/10; %迭代公式 f=diff(s); x=0; %迭代初值 a=subs(f,x); %判断收敛性(a是否小于1) s=(2-exp(x))/10; stop=1.0e-5; %迭代的精度 n=0; while a<1&abs(s-x)>=stop&n<=100; x=s %迭代 s=(2-exp(x))/10; n=n+1 end 实验结果： （1）二分法： symbol =1 b =0.50000000000000 n =1 symbol =1 b =0.25000000000000 n =2 symbol =1 b =0.12500000000000 n =3 symbol =0 a =0.06250000000000 n =4 symbol =1 b =0.09375000000000 n =5 symbol =0 a =0.07812500000000 n =6 symbol =1 b =0.09054565429688 n =15 symbol =1 b =0.09053039550781 n =16 symbol =0 a =0.09052276611328 n =17 Xk =0.09052276611328 n =17 （2）迭代法 由x =0.10000000000000 n =1 x =0.08948290819244 n =2 x =0.09063913585958 n =3 x =0.09051261667437 n =4 x =0.09052646805264 n =5 试验结果可见用二分法需要算17次，而用迭代法求得同样精度的解仅用5次，但由于迭代法一般只具有局部收敛性，因此通常不用二分法来求得非线性方程的精确解，而只用它求得根的一个近似解，再用收敛速度较快的迭代法求得其精确解。

MATLAB数学实验第二版答案(胡良剑)

数学实验答案 Chapter 1 Page20,ex1 (5) 等于[exp(1),exp(2);exp(3),exp(4)] (7) 3=1*3, 8=2*4 (8) a为各列最小值，b为最小值所在的行号 (10) 1>=4,false, 2>=3,false, 3>=2, ture, 4>=1,ture (11) 答案表明：编址第2元素满足不等式(30>=20)和编址第4元素满足不等式(40>=10) (12) 答案表明：编址第2行第1列元素满足不等式(30>=20)和编址第2行第2列元素满足不等式(40>=10) Page20, ex2 (1)a, b, c的值尽管都是1，但数据类型分别为数值，字符，逻辑，注意a与c相等，但他们不等于b (2)double(fun)输出的分别是字符a,b,s,(,x,)的ASCII码 Page20,ex3 >> r=2;p=0.5;n=12; >> T=log(r)/n/log(1+0.01*p) Page20,ex4 >> x=-2:0.05:2;f=x.^4-2.^x; >> [fmin,min_index]=min(f) 最小值最小值点编址 >> x(min_index) ans = 0.6500 最小值点 >> [f1,x1_index]=min(abs(f)) 求近似根--绝对值最小的点 f1 = 0.0328 x1_index = 24 >> x(x1_index) ans = -0.8500 >> x(x1_index)=[];f=x.^4-2.^x; 删去绝对值最小的点以求函数绝对值次小的点 >> [f2,x2_index]=min(abs(f)) 求另一近似根--函数绝对值次小的点 f2 = 0.0630 x2_index = 65 >> x(x2_index) ans = 1.2500

数学软件MATLAB实验作业

数学软件与数学实验作业 一.《数学软件》练习题（任选12题，其中19－24题至少选2题）： 3.对下列各式进行因式分解. (1). syms x y >> factor(x^5-x^3) (2). syms x y >> factor(x^4-y^4) (3). syms x >> factor(16-x^4) (4). syms x >> factor(x^3-6*x^2+11*x-6) (5). syms x y >> factor((x+y)^2-10*(x+y)+25) (6). syms x y >> factor(x^2/4+x*y+y^2) (7). syms x y a b >> factor(3*a*x+4*b*y+4*a*y+3*b*x) (8). syms x >> factor(x^4+4*x^3-19*x^2-46*x+120) 5.解下列方程或方程组. (1).solve('(y-3)^2-(y+3)^3=9*y*(1-2*y)') (2). solve('3*x^2+5*(2*x+1)') (3). solve('a*b*x^2+(a^4+b^4)*x+a^3*b^3','x') (4). solve('x^2-(2*m+1)*x+m^2+m','x') (5). [x,y]=solve('4*x^2-9*y^2=15','2*x-3*y=15') 6.计算极限. (1). syms x f=(exp(x)-exp(-x))/sin(x); limit(f,x,0) (2) syms x >> f=(x/(x-1)-1/log(x)); >> limit(f,x,1) (3). syms x >> f=(1-cos(x))/x^2; >> limit(f,x,0)

MATLAB基础教程 薛山第二版 课后习题答案

《MATLAB及应用》实验指导书《MATLAB及应用》实验指导书 班级：T1243-7 姓名：柏元强 学号：20120430724 总评成绩： 汽车工程学院 电测与汽车数字应用中心

目录 实验04051001 MATLAB语言基础 (1) 实验04051002 MATLAB科学计算及绘图 (18) 实验04051003 MATLAB综合实例编程 (31)

实验04051001 MATLAB语言基础 1实验目的 1)熟悉MATLAB的运行环境 2)掌握MATLAB的矩阵和数组的运算 3)掌握MATLAB符号表达式的创建 4)熟悉符号方程的求解 2实验内容 第二章 1.创建double的变量，并进行计算。 (1)a=87，b=190，计算 a+b、a-b、a*b。 clear,clc a=double(87); b=double(190); a+b,a-b,a*b (2)创建 uint8 类型的变量，数值与(1)中相同，进行相同的计算。 clear,clc a=uint8(87); b=uint8(190); a+b,a-b,a*b 2．计算：

(1) () sin 60 (2) e3 (3) 3cos 4??π ??? clear,clc a=sind(60) b=exp(3) c=cos(3*pi/4) 3．设2u =，3v =，计算： (1) 4 log uv v (2) () 2 2 e u v v u +- (3) clear,clc u=2;v=3; a=(4*u*v)/log(v) b=((exp(u)+v)^2)/(v^2-u) c=(sqrt(u-3*v))/(u*v) 4．计算如下表达式： (1) ()() 3542i i -+ (2) () sin 28i - clear,clc (3-5*i)*(4+2*i) sin(2-8*i)

matlab数学实验练习题

Matlab 数学实验 实验一 插值与拟合 实验内容： 预备知识：编制计算拉格朗日插值的M 文件。 1. 选择一些函数，在n 个节点上（n 不要太大，如5 ~ 11）用拉格朗日、分段线性、三次样条三种插值方法，计算m 个插值点的函数值（m 要适中，如50~100）。通过数值和图形输出，将三种插值结果与精确值进行比较。适当增加n ，再做比较，由此作初步分析。下列函数任选一种。 （1）、 ;20,sin π≤≤=x x y （2）、;11,)1(2/12≤≤--=x x y （3）、;22,c o s 10 ≤≤-=x x y （4）、22),exp(2≤≤--=x x y 2.用电压V=10伏的电池给电容器充电，电容器上t 时刻的电压为 ) (0)()(t e V V V t v ---=，其中0V 是电容器的初始电压，τ是充电常数。试由下面 一组t ，V 数据确定0V 和τ。 实验二 常微分方程数值解试验 实验目的： 1. 用MATLAB 软件求解微分方程，掌握Euler 方法和龙格-库塔方法； 2. 掌握用微分方程模型解决简化的实际问题。 实验内容：

实验三地图问题 1.下图是一个国家的地图，为了计算出它的国土面积，首先对地图作如下测量： 以由西向东方向为x轴，由南到北方向为y轴，选择方便的原点，并将从最西边界点到最东边界点在x轴上的区间适当地划分为若干段，在每个分点的y方向测出南边界点和北边界点的y坐标y1和y2，这样就得到了表中的数据（单位mm）。 根据地图的比例我们知道18mm相当于40km，试由测量数据计算该国土 的近似面积，并与它的精确值41288km2比较。

MATLAB实验报告

实验一 MATLAB 环境的熟悉与基本运算 一、实验目的及要求 1．熟悉MATLAB 的开发环境； 2．掌握MATLAB 的一些常用命令； 3．掌握矩阵、变量、表达式的输入方法及各种基本运算。 二、实验内容 1.熟悉MATLAB 的开发环境： ① MATLAB 的各种窗口： 命令窗口、命令历史窗口、工作空间窗口、当前路径窗口。 ②路径的设置： 建立自己的文件夹，加入到MATLAB 路径中，并保存。 设置当前路径，以方便文件管理。 2.学习使用clc 、clear ，了解其功能和作用。 3.矩阵运算： 已知:A=[1 2;3 4]; B=[5 5;7 8]; 求:A*B 、A.*B ，并比较结果。 4.使用冒号选出指定元素: 已知:A=[1 2 3;4 5 6;7 8 9]; 求:A 中第3列前2个元素；A 中所有列第2，3行的元素； 5.在MATLAB 的命令窗口计算: 1) )2sin(π 2) 5.4)4.05589(÷?+ 6.关系及逻辑运算 1)已知:a=[5:1:15]; b=[1 2 8 8 7 10 12 11 13 14 15]，求: y=a==b ，并分析结果 2)已知:X=[0 1;1 0]; Y=[0 0;1 0]，求: x&y+x>y ，并分析结果 7.文件操作 1)将0到1000的所有整数，写入到D 盘下的文件 2)读入D 盘下的文件，并赋给变量num

8.符号运算 1)对表达式f=x 3 -1 进行因式分解 2)对表达式f=（2x 2*(x+3)-10)*t ，分别将自变量x 和t 的同类项合并 3）求 3(1)x dz z +? 三、实验报告要求 完成实验内容的3、4、5、6、7、8，写出相应的程序、结果

MATLAB数学实验100例题解

一元函数微分学 实验1 一元函数的图形（基础实验） 实验目的 通过图形加深对函数及其性质的认识与理解, 掌握运用函数的图形来观察和分析 函数的有关特性与变化趋势的方法,建立数形结合的思想; 掌握用Matlab 作平面曲线图性的方法与技巧. 初等函数的图形 2 作出函数x y tan =和x y cot =的图形观察其周期性和变化趋势. 解：程序代码： >> x=linspace(0,2*pi,600); t=sin(x)./(cos(x)+eps); plot(x,t);title('tan(x)');axis ([0,2*pi,-50,50]); 图象： 程序代码： >> x=linspace(0,2*pi,100); ct=cos(x)./(sin(x)+eps); plot(x,ct);title('cot(x)');axis ([0,2*pi,-50,50]); 图象： cot(x) 4在区间]1,1[-画出函数x y 1 sin =的图形. 解：程序代码： >> x=linspace(-1,1,10000); y=sin(1./x); plot(x,y); axis([-1,1,-2,2]) 图象：

二维参数方程作图 6画出参数方程???==t t t y t t t x 3cos sin )(5cos cos )(的图形： 解：程序代码： >> t=linspace(0,2*pi,100); plot(cos(t).*cos(5*t),sin(t).*cos(3*t)); 图象： 极坐标方程作图 8 作出极坐标方程为10/t e r =的对数螺线的图形. 解：程序代码： >> t=0:0.01:2*pi; r=exp(t/10); polar(log(t+eps),log(r+eps)); 图象： 90270 分段函数作图 10 作出符号函数x y sgn =的图形. 解：

Matlab程序设计与应用第二版刘卫国课后实验答案

. 实验一： T1： %%第一小题 z1=2*sin(85*pi/180)/(1+exp(2)) %%第二小题 x=[2,1+2i;-0.45,5]; z2=1/2*log(x+sqrt(1+x.^2)); z2 %%第三小题 a=-3.0:0.1:3.0; z3=1/2*(exp(0.3*a)-exp(-0.3*a)).*sin(a+ 0.3)+log((0.3+a)/2) %%第四题 t=0:0.5:2.5 z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t .^2-1)+(t>=2&t<3).*(t.^2-2*t+1) T2： A=[12,34,-4;34,7,87;3,65,7] B=[1,3,-1;2,0,3;3,-2,7] disp ('A+6*B='); disp(A+6*B); disp('A-B+I=');disp(A-B+eye(3)); disp('A*B='); disp(A*B); disp('A.*B='); disp(A.*B); disp('A^3='); disp(A^3); disp('A.^3='); disp(A.^3); disp('A/B='); disp(A/B); disp('B\A='); disp(B\A); disp('[A,B]='); disp([A,B]); disp('[A([1,3],:);B^2]='); disp([A([1,3],:);B^2]); T3: z=1:25; A=reshape(z,5,5)'; B=[3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11]; C=A*B

浅析Matlab数学实验报告

数学实验报告 姓名: 班级: 学号: 第一次实验任务 过程: a=1+3i; b=2-i; 结果: a+b =3.0000 + 2.0000i a-b =-1.0000 + 4.0000i a*b = 5.0000 + 5.0000i a/b = -0.2000 + 1.4000i 过程: x=-4.5*pi/180; y=7.6*pi/180; 结果: sin(abs(x)+y)/sqrt(cos(abs(x+y))) =0.2098 心得:对于matlab 中的角度计算应转为弧度。 (1)过程: x=0:0.01:2*pi; y1=sin(x); y2=cos(x); y3=exp(x); y4=log(x); plot(x,y1,x,y2,x,y3,x,y4) plot(x,y1,x,y2,x,y3,x,y4) 结果: (2)过程:>> subplot(2,2,1) >> plot(x,y1) >> subplot(2,2,2) >> plot(x,y2) ./,,,,2,311b a b a b a b a i b i a ?-+-=+=计算、设有两个复数 6,7,5.4)

cos()sin(2=-=++y x y x y x ,其中、计算的图形。 下分别绘制)同一页面四个坐标系)同一坐标系下(、在( x y e y x y x y x ln ,,cos ,sin 213==== >> subplot(2,2,3) >> plot(x,y3) >> subplot(2.2.4) >> subplot(2,2,4) >> plot(x,y4) 结果: 心得:在matlab中,用subplot能够实现在同一页面输出多个坐标系的图像,应注意将它与hold on进行区别,后者为在同一坐标系中划出多条曲线。 5、随机生成一个3x3矩阵A及3x2矩阵B,计算(1)AB,(2)对B中每个元素平方后得到的矩阵C,(3)sinB,(4)A的行列式,(5)判断A是否可逆,若可逆,计算A的逆矩阵,(6)解矩阵方程AX=B,(7)矩阵A中第二行元素加1,其余元素不变,得到矩阵D,计算D。 过程:A=fix(rand(3,3).*10) ; B=fix(rand(3,3).*10);

matlab与数学实验大作业

《数学实验与MATLAB》 ——综合实验报告 实验名称：不同温度下PDLC薄膜的通透性 与驱动电压的具体关系式的研究学院：计算机与通信工程学院 专业班级： 姓名： 学号： 同组同学： 2014年 6月10日

一、问题引入 聚合物分散液晶(PDLC)是将低分子液晶与预聚物Kuer UV65胶相混合，在一定条件下经聚合反应，形成微米级的液晶微滴均匀地分散在高分子网络中，再利用液晶分子的介电各向异性获得具有电光响应特性的材料，它主要工作在散射态和透明态之间并具有一定的灰度。聚合物分散液晶膜是将液晶和聚合物结合得到的一种综合性能优异的膜材料。该膜材料能够通过驱动电压来控制其通透性，可以用来制作PDLC型液晶显示器等，具有较大的应用范围。已知PDLC薄膜在相同光强度及驱动电压下，不用的温度对应于不同的通透性，不同温度下的阀值电压也不相同。为了尽量得到不同通透性的PDLC薄膜，有必要进行温度对PDLC薄膜的特性的影响的研究。现有不同温度下PDLC 薄膜透过率与驱动电压的一系列数据，试得出不同温度下PDLC薄膜通透性与驱动电压的具体关系式，使得可以迅速得出在不同温度下一定通透性对应的驱动电压。 二、问题分析 想要得到不同温度下PDLC薄膜通透性与驱动电压的具体关系式可以运用MATLAB多项式农合找出最佳函数式，而运用MATLAB多项式插值可以得出在不同温度下一定通透性所对应的驱动电压。 三、实验数据 选择10、20、30摄氏度三个不同温度，其他条件一致。

（1）、10摄氏度 实验程序： x=2:2:40; y=[5.2,5.4,5.8,6.4,7.2,8.2,9.4,10.8,12.2,14.0,16.6,22.0, 30.4,39.8,51.3,55.0,57.5,58.8,59.6,60.2]; p3=polyfit(x,y,3); p5=polyfit(x,y,5); p7=polyfit(x,y,7); disp('三次拟合函数'),f3=poly2str(p3,'x') disp('五次拟合函数'),f5=poly2str(p5,'x') disp('七次拟合函数'),f7=poly2str(p7,'x') x1=0:1:40; y3=polyval(p3,x1); y5=polyval(p5,x1); y7=polyval(p7,x1); plot(x,y,'rp',x1,y3,'--',x1,y5,'k-.',x1,y7); legend('拟合点','三次拟合','五次拟合','七次拟合') 实验结果：

南邮MATLAB数学实验答案(全)

第一次练习 教学要求：熟练掌握Matlab 软件的基本命令和操作，会作二维、三维几何图形，能够用Matlab 软件解决微积分、线性代数与解析几何中的计算问题。 补充命令 vpa(x,n) 显示x 的n 位有效数字，教材102页 fplot(‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形 在下面的题目中m 为你的学号的后3位（1-9班）或4位（10班以上） 1.1 计算30sin lim x mx mx x →-与3 sin lim x mx mx x →∞- syms x limit((902*x-sin(902*x))/x^3) ans = 366935404/3 limit((902*x-sin(902*x))/x^3,inf) ans = 0 1.2 cos 1000 x mx y e =，求''y syms x diff(exp(x)*cos(902*x/1000),2) ans = (46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算 22 11 00 x y e dxdy +?? dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1) ans = 2.1394 1.4 计算4 2 2 4x dx m x +? syms x int(x^4/(902^2+4*x^2)) ans = (91733851*atan(x/451))/4 - (203401*x)/4 + x^3/12 1.5 (10)cos ,x y e mx y =求 syms x diff(exp(x)*cos(902*x),10) ans = -356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x) 1.6 0x =的泰勒展式(最高次幂为4).

MATLAB程序设计及应用(第二版)课后实验答案

Matlab 课后实验题答案 实验一 MATLAB 运算基础 1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。 (1) 0 12 2sin851z e =+ (2) 221 ln(1)2z x x = ++，其中2120.45 5i x +??=??-?? (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+= ++=-- (4) 22 42011 122123t t z t t t t t ?≤=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3) .*(t.^2-2*t+1) 2. 已知：

1234413134787,2033657327A B --???? ????==???? ????-???? 求下列表达式的值： (1) A+6*B 和A-B+I （其中I 为单位矩阵） (2) A*B 和A.*B (3) A^3和A.^3 (4) A/B 及B\A (5) [A,B]和[A([1,3],:);B^2] 解： M 文件： A=[12 34 -4;34 7 87;3 65 7];B=[1 3 -1;2 0 3;3 -2 7]; A+6.*B A-B+eye(3) A*B A.*B A^3 A.^3 A/B B\A [A,B] [A([1,3],:);B^2] 3. 设有矩阵A 和B 1234 53 166789101769,11 121314150 23416171819209 7021222324254 1311A B ???? ????-??? ?????==-??? ? ???????????? (1) 求它们的乘积C 。 (2) 将矩阵C 的右下角3×2子矩阵赋给D 。 (3) 查看MATLAB 工作空间的使用情况。 解：. 运算结果： E=(reshape(1:1:25,5,5))';F=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11]; C= E*F H=C(3:5,2:3) C = 93 150 77

MATLAB数学实验报告

Matlab 数学实验报告

一、实验目的 通过以下四组实验,熟悉MATLAB的编程技巧,学会运用MATLAB的一些主要功能、命令，通过建立数学模型解决理论或实际问题。了解诸如分岔、混沌等概念、学会建立Malthu模型和Logistic 模型、懂得最小二乘法、线性规划等基本思想。 二、实验内容 2.1实验题目一 2.1.1实验问题 Feigenbaum曾对超越函数y=λsin(πｘ)（λ为非负实数）进行了分岔与混沌的研究，试进行迭代格式x k+1=λsin(πx k),做出相应的Feigenbaum图 2.1.2程序设计 clear;clf; axis([0,4,0,4]); hold on for r=0:0.3:3.9 x=[0.1]; for i=2:150 x(i)=r*sin(3.14*x(i-1)); end pause(0.5) for i=101:150

plot(r,x(i),'k.'); end text(r-0.1,max(x(101:150))+0.05,['\it{r}=',num2str(r)]) end 加密迭代后 clear;clf; axis([0,4,0,4]); hold on for r=0:0.005:3.9 x=[0.1];

for i=2:150 x(i)=r*sin(3.14*x(i-1)); end pause(0.1) for i=101:150 plot(r,x(i),'k.'); end end 运行后得到Feigenbaum图

2.2实验题目二 2.2.1实验问题 某农夫有一个半径10米的圆形牛栏，长满了草。他要将一头牛拴在牛栏边界的桩栏上，但只让牛吃到一半草，问拴牛鼻子的绳子应为多长？ 2.2.2问题分析 如图所示，E为圆ABD的圆心，AB为拴牛的绳子，圆ABD为草场，区域ABCD为牛能到达的区域。问题要求区域ABCD等于圆ABC 的一半，可以设BC等于x,只要求出∠a和∠b就能求出所求面积。先计算扇形ABCD的面积，2a÷π×πx2=2aπ2，再求AB的面积，用扇形ABE的面积减去三角形ABE的面积即可。

MATLAB程序设计与应用(第二版)刘卫国主编_部分实验答案

实验六 2_1 clear; x=linspace(0,2*pi,101); y1=x.^2; y2=cos(2.*x); y3=y1.*y2; plot(x,y1,'b-',x,y2,'r:',x,y3,'g-.'); %y1蓝色实线，y2红色虚线，y3绿色点画线 2_2 subplot(2,2,1); %分四个子图(先画2行2列第1块) plot(x,y1); subplot(2,2,2); plot(x,y2), subplot(2,2,3); plot(x,y3);

2_3 () subplot(3,4,1); %y1的四种图形bar(x,y1); subplot(3,4,2); stairs(x,y1), subplot(3,4,3); stem(x,y1); subplot(3,4,4); fill(x,y1,'b'); subplot(3,4,5); %y2 bar(x,y2); %条形图subplot(3,4,6); stairs(x,y2), %阶梯图subplot(3,4,7); stem(x,y2); %杆图subplot(3,4,8);

fill(x,y2,'b'); %填充图，注意必须加填充颜色 subplot(3,4,9); %y3 bar(x,y3); subplot(3,4,10); stairs(x,y3), subplot(3,4,11); stem(x,y3); subplot(3,4,12); fill(x,y3,'b'); 3 clear; x=-5:0.1:5; if x<=0 y=(x+sqrt(pi)/exp(2)); else y=0.5.*log(x+sqrt(1+x.^2));

数学实验与数学软件(Mathmaticandmatlab)

数学软件与数学实验2013-2014学年度秋季学期期末试卷 专业：统计学 班级：11级2班 学号：20110723 姓名：晏静

一、按要求计算出下列表达式的值 (1)318, 3 162 53 ?? + ? ?? , 21 eπ+, 2.5 tg, 2 log15； (2)给出π的9位和e的10位近似值； (3)求658和4102的最大公约数及35和25的最小公倍数； (4)产生10个0与10之间随机数的一个表; (5)求虚数1453 i i i i +- -的实部,虚部,模,共轭,辐角。 (6)自己运用Table建立两个表，并进行表运算，如连接、并集、交、排序等操作。

二、因式分解 22212321332112322 1 22(1)()()()4;(2)21;x x x x x x x x x x x x x x x +++++---- 解： 三、解方程(组) 1234234124234-2+344-+-3(1)+31-73+3 x x x x x x x x x x x x x -=??=? ? +=??+=-? 65432(2)5232002000.x x x x x x -+--++= 四、求极限 () 20 (1)1sin ;(2);(3)56! ctg x n x n n n Lim x Lim n n →→∞ →∞++

(1) (2) (3) 五、求导数 32 22(1)()=ln(x+1+);(2)()=cos 2,; (3)=log (),Z . x f x x f f x e y x y Z xy x y y ???求的导数已知求求关于的二阶导 (1) (2) (3) 六、求下列定积分与不定积分： ()()()12201+sin ln 1+(1);(2);(3)sin (1+cos ) +1(1+)(2+-) x x dx dx x x x x x x ? ? ?2 2-(4)=0,=1,==.y D D x y y x I x e d σ??设是由直线围成的区域，计算的值 (1) (2)

Matlab数学实验一2015(标准答案版)

Matｌab数学实验一——matlab初体验 一、实验目的及意义 ［1] 熟悉MAＴＬＡB软件的用户环境； ［2] 了解MAＴLAB软件的一般目的命令； ［３] 掌握MATLＡB数组操作与运算函数； 通过该实验的学习,使学生能熟悉matｌab的基础应用，初步应用ＭATLAB软件解决一些简单问题。 二、实验内容 １.认识ｍatｌab的界面和基本操作 2.了解mａtlab的数据输出方式（ｆormat) 3. MATLAＢ软件的数组(矩阵）操作及运算练习; 三、实验任务 根据实验内容和步骤,完成以下具体实验，要求写出实验报告（实验目的→问题→原理→算法与编程→计算结果或图形→心得体会) 完成如下题目,并按照实验报告格式和要求填写实验报告 1．在ｃommanｄwinｄow中分别输入如下值,看它们的值等于多少，并用maｔｌａb的ｈelp中查询这些缺省预定义变量的含义,用中文写出它们的意义。 iｊeps inf nａn pi realｍaｘrealｍin 2.分别输入一个分数、整数、小数等，(如:a=1/9)，观察显示结果，并使用forｍaｔ函数控制数据的显示格式，如：分别输入ｆｏｒmat sｈorｔ、fｏｒmat loｎg、format short e、forｍａt lｏng g、foｒmａt bank、forｍａt heｘ等，然后再在命令窗口中输入ａ,显示a的值的不同形式，并理解这些格式的含义。 3．测试函数clｅａr、clc的含义及所带参数的含义（利用ｍatlａb的ｈelp功能）。 4. 写出在命令窗口中的计算步骤和运行结果。 （1)计算 1.22 10 (ln log) 81 e ππ +- ； >＞(log(pi)+lｏｇ(pi)/lｏg(10)-ｅxp(1.２)）^２/81 >>ａns = ０．０3４8 (2) >> x=2;y＝４； >> z=x^２+exｐ(x+y)-y*lｏg(x)-3 z ＝ 40１.65６2 （3)输入变量 13 5.3, 25 a b ?? ==?? ?? ，在工作空间中使用who，whos，并用ｓave命令将变量存入”D：\exe0 １.ｍat”文件。测试clear命令，然后用load命令将保存的”D：\ｅｘe０１.ｍat”文件载入＞> ａ＝5.3 ａ＝

Matlab程序设计与应用第二版刘卫国课后实验答案

实验一： T1： %%第一小题 z1=2*sin(85*pi/180)/(1+exp(2)) %%第二小题 x=[2,1+2i;-0.45,5]; z2=1/2*log(x+sqrt(1+x.^2)); z2 %%第三小题 a=-3.0:0.1:3.0; z3=1/2*(exp(0.3*a)-exp(-0.3*a)).*sin(a+0.3)+l og((0.3+a)/2) %%第四题 t=0:0.5:2.5 z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t >=2&t<3).*(t.^2-2*t+1) T2： A=[12,34,-4;34,7,87;3,65,7] B=[1,3,-1;2,0,3;3,-2,7] disp ('A+6*B='); disp(A+6*B); disp('A-B+I='); disp(A-B+eye(3)); disp('A*B='); disp(A*B); disp('A.*B='); disp(A.*B); disp('A^3='); disp(A^3); disp('A.^3='); disp(A.^3); disp('A/B='); disp(A/B); disp('B\A='); disp(B\A); disp('[A,B]='); disp([A,B]); disp('[A([1,3],:);B^2]='); disp([A([1,3],:);B^2]); T3: z=1:25; A=reshape(z,5,5)'; B=[3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11]; C=A*B D=C(3:5,2:3) T4-1: a=100:999; b=find(rem(a,21)==0); c=length(b) T4-2: a=input('请输入一个字符串：','s'); b=find(a>='A'&a<='Z'); a(b)=[]; disp(a); 实验二： T1： E=eye(3),R=rand(3,2),O=zeros(2,3),S=diag([1,2 ]); A=[E,R;O,S] disp('A^2='); disp(A^2); disp('[E,R+RS;O,S^2]'); B=[E,R+R*S;O,S^2] T2： H=hilb(5) P=pascal(5) Hh=det(H) Hp=det(P) Th=cond(H) Tp=cond(P) a=abs(Th-1); b=abs(Tp-1); if a>b disp('帕萨卡矩阵P性能更好'); elseif a

数学实验matlab练习题

2015-2016数学实验练习题 一、选择题 1.清除Matlab工作空间（wordspace）变量的命令是（B ） A. clc B. clear C. clf D.delete 2. 清除当前屏幕上显示的所有内容，但不清除工作空间中的数据的命令是（ A ） A. clc B. clear C. clf D.delete 3. 用来清除图形的命令（ C ） A. clc B. clear C. clf D.delete 4. 在MATLAB程序中，使命令行不显示运算结果的符号是（ A ） A. ; B. % C. # D. & 5. 在MATLAB程序中，可以将某行表示为注释行的符号是（ B ） A. ; B. % C. # D. & 6.在循环结构中跳出循环，执行循环后面代码的命令为 ( B ) A. return B. break C. continue D. Keyboard 7.在循环结构中跳出循环，但继续下次循环的命令为（ C ） A. return B. break C. continue D. Keyboard 8. MATLAB中用于声明全局变量的关键字是( C ) A. inf B. syms C. global D. function 9. 用户可以通过下面哪项获得指令的使用说明（ A ） A. help B. load C. demo D. lookfor 10．在MATLAB命令窗口中键入命令S=zoros(3)；可生成一个三行三列的零矩阵，如果省略了变量名S，MATLAB表现计算结果将用下面的哪一变量名做缺省变量名（ A ） A. ans; B. pi; C. NaN; D. Eps. 11. 9/0的结果是（ B ） A. NAN； B. Inf； C. eps； D. 0 12．在MATLAB中程序或语句的执行结果都可以用不同格式显示，将数据结果显示为分数形式，用下面哪一条命令语句（ D ） A. format long； B. format long e； C. format bank； D. fromat rat 13. 下列MATLAB命令中是构造1行3列的(-1,1)均匀分布随机矩阵的命令的是（D）

MATLAB与数学实验练习

一、填空 1、命令clear、clf、clc、who、whos的含义分别是clear用于 清除内存中的所有变量与函数；clf用于清除图形窗口； clc用于清除命令窗口中的所有显示内容；who将内存中 的当前变量以简单的形式列出；whos列出当前内存变量 的名称、大小和类型等信息。 2、若矩阵 2310 , 5112 A B ???? == ? ? ???? ,则A.*B= （对应元素相乘）（2 0； 5 2） ,A*B=（矩阵的乘法）（5 6； 6 2）。 3、生成4行3列的元素全为1的矩阵的命令是 ones(4,3) ；生成4阶全零矩阵的命令是 zeros(4) 。生成对角线元素为[2 4 6]矩阵的命令是v=[2 4 6] A=diag(v,0) 。 4、假定A是一个6阶方阵，选取矩阵A第3行的指令是A(3,:) ，选取第4行第2列元素的指令是A(4,2) 。 5、生成一个从1到50，步长为10的等差数列构成的数组，可以使用的命令是 (from:step:to) 1:10:50 。 6、求x的平方根使用的命令是sqrt(x)，round(pi)的结果 是（取整）3 7、命令subplot(2,2,2)的功能是把图形窗口分为2*2=4 个子图，并把第 2 个子图作为当前图形窗口。 8、matlab中，求解线性方程组时，矩阵行变换化简的命令是

reff 。 9、在用命令 p=polyfit （x ，y ，n ）对数据进行多项式拟合时，参数n 的含义为 n 次多项式 。 10、求函数4y x cos7x =的5阶导数的命令是 diff （y ，5） 。 11、用符号法求解微分方程2xy -3y =x ,y(1)=0,y(5)=0'''的解的 指令是 dsolve(‘x*D2y-3*Dy=x^2’,’y(1)=0,y(5)=0’,’x ’) 12、Matlab 可以输入字母、汉字，但是M 文件中标点符号必须（英文） 状态下输入。 14、求x e 的命令是（exp(x)），求x 的自然对数lnx 的命令（log(x )）。 15、画图时，在x 轴旁边加注文字说明的命令是 (xlabel(‘string ’) )，图名的标注命令是（title(‘string ’)）， 图例标注命令是（legend(‘string ’,’string ’,……)），在鼠标指定 位置上标注的命令是（gtext(‘string ’)），将一个图形窗口分成多 个子图的命令是（subplot(m,n,i)），画空间曲线的命令是 （plot3(x,y,z)）。 16、绘制三维空间曲面图形的命令是（mesh ）和（surf ），生成格 点矩阵的命令是( meshgrid ) 17、用矩形法、复合梯形公式、复合辛普生公式求12 4dx 1x +?的定积分。(详解见P150 7-28) h=0.01;x=0:h:1; y=4./(1+x.^2); format long