matlab實驗二答案福州大學

A. 跪求MATLAB課後習題答案，求好心人啊！

1，下面為m文件代碼
i=-0.1;
j=1;
if i==0
i=i+0.005;%此步為了避免fzero（f，0）解不出
end
while(i<=0.1)
f=@(x)x*sin(1/x);%創建函數句柄
x(j)=fzero(f,i);%用fzero函數求解i附近的零點
j=j+1;
i=i+0.005;

end
x%此時x為一向量

2，
>> length=@(x)5./cos(x)+10./sin(x)

length =

@(x)5./cos(x)+10./sin(x)

>> [x,fval]= fminbnd(length,0,pi/2)

x =

0.8999

fval =

20.8097

>>length2=@(x)sqrt((5./cos(x)+10./sin(x))^2+10^2)

>> [x,fval]= fminbnd(length2,0,pi/2)

x =

0.8999

fval =

23.0877

B. 誰有matlab程序設計與應用（第二版）劉衛國主編的課後實驗一到實驗十五的全部答案急急急

MATLAB程序設計與應用（第二版）實驗參考答案
%實驗一 MATLAB運算基礎
%第一題
%（1）
z1=2*sin(85*pi/180)/(1+exp(2))
%（2）
x=[2,1+2i;-0.45,5];
z2=0.5*log(x+sqrt(1+x.^2))
%(3)
a=-3.0:0.1:3.0;
z3=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2)
%(4)
t=0:0.5:2.5;
z4=t.^2.*(t>=0&t<1)+(t.^2-1).*(t>=1&t<2)+(t.^2-2*t+1).*(t>=2&t<3)

%第二題
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(size(A))
A*B
A.*B
A^3
A.^3
A/B
B\A
[A,B]
[A([1,3],:);B^2]

%第三題
A=[1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20;21 22 23 24 25]
B=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11]
C=A*B
F=size(C)
D=C(F(1)-2:F(1),F(2)-1:F(2))
whos

%第四題
%(1):
A=100:999;
B=rem(A,21);
C=length(find(B==0))
%(2):
A='lsdhKSDLKklsdkl';
k=find(A>='A'&A<='Z');
A(k)=[]

%實驗二 MATLAB矩陣分析與處理
%第一題
E=eye(3);
R=rand(3,2);
O=zeros(2,3);
S=diag([2,3]);
A=[E,R;O,S];
A^2
B=[E,(R+R*S);O,S^2]
%第二題
H=hilb(5)
P=pascal(5)
Hh=det(H)
Hp=det(P)
Th=cond(H)
Tp=cond(P)
%第三題:
A=fix(10*rand(5))
H=det(A)
Trace=trace(A)
Rank=rank(A)
Norm=norm(A)
%第四題:
A=[-29,6,18;20,5,12;-8,8,5]
[V,D]=eig(A)
%數學意義略
%第五題方法一:
%(1):
A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6];
b=[0.95,0.67,0.52]';
x=inv(A)*b
%(2):
B=[0.95,0.67,0.53]';
x=inv(A)*B
%(3):
cond(A)
%第五題方法二:
A=hilb(4)
A(:,1)=[]
A(4,:)=[]
B=[0.95,0.67,0.52]';
X=inv(A)*B
B1=[0.95,0.67,0.53]';
X1=inv(A)*B1
N=cond(B)
N1=cond(B1)
Na=cond(A) %矩陣A為病態矩陣
%第六題
A=[1,4,9;16,25,36;49,64,81]
B=sqrtm(A)
C=sqrt(A) %sqrtm函數是以矩陣為單位進行計算，sqrt函數是以矩陣中的元素進行計算

%實驗三 選擇程序結構設計
%第一題程序一
x=[-5.0,-3.0,1.0,2.0,2.5,3.0,5.0];
y=[]; %建立存放所有y值的矩陣
for x0=x
if x0<0&x0~=-3
y=[y,x0*x0+x0-6];
elseif x0>=0&x0<5&x0~=2&x0~=3
y=[y,x0*x0-5*x0+6];
else
y=[y,x0*x0-x0-1];
end
end
x %輸出所有x
y %輸出所有y
%第一題程序二
x=[-5,-3,1,2,2.5,3,5];
y=[];
for a=1:7
if x(a)<0&x(a)~=-3
y=[y,(x(a))^2+x(a)-6];
elseif x(a)>=0&x(a)<5&x(a)~=2&x(a)~=3
y=[y,(x(a))^2-5*x(a)+6];
else
y=[y,x(a)*x(a)-x(a)-1];
end
end
%第二題程序一:
x=input('請輸入一個百分製成績：');
if x>100|x<0
disp('您輸入的成績不是百分製成績，請重新輸入。');
else
if x<=100&x>=90
disp('A');
elseif x<=89&x>=80
disp('B');
elseif x<=79&x>=70
disp('C');
elseif x<=69&x>60
disp('D');
else
disp('E');
end
end
%第二題程序二：
s=input('請輸入一個成績（0分到100分之間）：'); %s用於存放成績
while 1 %判斷輸入成績的合理性
if s<0|s>100
disp('輸入的成績需在0到100之間，請重新輸入：')
s=input('請輸入一個成績（0分到100分之間）：');
else
break;
end
end
switch fix(s/10) %對成績做出等級判斷
case {9,10}
disp('A')
case 8
disp('B')
case 7
disp('C')
case 6
disp('D')
otherwise
disp('E')
end
%第三題
n=input('請輸入員工工號：');
h=input('該員工工作時數是：');
if h>120
x=(h-120)*84*(1+0.15)+120*84;
elseif h<60
x=h*84-700;
else
x=h*84;
end
disp([num2str(n),'號員工','的應發工資為',num2str(x)]);
%第四題（還可以用switch語句實現）
a=fix(10+(99-10)*rand(1,2)) %產生兩個隨機整數
x=a(1);
y=a(2);
t=input('請輸入運算符號：','s');
if t=='+'
z=x+y;
elseif t=='-'
z=x-y;
elseif t=='*'
z=x*y;
elseif t=='/'
z=x/y;
end
disp([num2str(x),t,num2str(y),'=',num2str(z)]) %輸出運算結果
%第五題
a=rand(5,6) %產生5x6的隨機矩陣
n=input('請輸入您要輸出矩陣的第幾行：');
if n>5
disp('超出了矩陣的行數，矩陣的最後一行為：')
a(5,:)
else
disp(['矩陣的第',num2str(n),'行為：'])
a(n,:)
end

%實驗四 循環結構程序設計
%第一題程序一
s=0;
n=input('n=?');
for i=1:n
s=s+1/i/i;
end
PI=sqrt(6*s)
pi
%第一題程序二
n=input('n=?');
a=1:n;
b=1./a.^2;
PI=sqrt(6*sum(b))
pi
%第二題
y=0;
n=1;
while(y<3)
y=y+1/(2*n-1);
n=n+1;
end
y=y-1/(2*(n-1)-1)
n=n-2
%第三題
a=input('a=?');
b=input('b=?');
Xn=1;
Xn1=a/(b+Xn);
n=0;
while abs(Xn1-Xn)>1e-5
Xn=Xn1;
Xn1=a/(b+Xn);
n=n+1;
if n==500
break;
end
end
n
Xn1
r1=(-b+sqrt(b*b+4*a))/2
r2=(-b-sqrt(b*b+4*a))/2
%第四題
for i=1:100
if i==1
f(i)=1;
elseif i==2
f(i)=0;
elseif i==3
f(i)=1;
else
f(i)=f(i-1)-2*f(i-2)+f(i-3);
end
end
max(f)
min(f)
sum(f)
length(find(f>0))
length(find(f==0))
length(find(f<0))
%第五題：
s=0;n=0;
for i=2:49
b=i*(i+1)-1;
m=fix(sqrt(b));
for j=2:m
if rem(b,j)==0
break
end
end
if j==m
n=n+1;
s=s+b;
end
end
n
s

%實驗五 函數文件
%第一題
function y=mat1(x) %建立函數文件mat1.m
y=[exp(x),log(x),sin(x),cos(x)];
%在命令窗口調用上述函數文件：
y=mat1(1+i)
%第二題程序一
function [a,b,N,M]=shiyanwu2(m,n,t)
A=[m*cos(t*pi/180),-m,-sin(t*pi/180),0;m*sin(t*pi/180),0,cos(t*pi/180),0;0,n,-sin(t*pi/180),0;0,0,-cos(t*pi/180),1];
B=[0,9.8*m,0,9.8*n];
C=inv(A)*B';
a=C(1);
b=C(2);
N=C(3);
M=C(4);
%在命令窗口調用該函數文件:
m1=input('m1=');
m2=input('m2=');
theta=input('theta=');
[a1,a2,N1,N2]=shiyanwu2(m1,m2,theta)
%第二題程序二
function X=mat2(m1,m2,t)
g=9.8;
A=[m1*cos(t*pi/180),-m1,-sin(t*pi/180),0;m1*sin(t*pi/180),0,cos(t*pi/180),0;0,m2,-sin(t*pi/180),0;0,0,-cos(t*pi/180),1];
B=[0;m1*g;0;m2*g];
X=inv(A)*B;
%在命令窗口調用該函數文件：
X=mat2(1,1,60)
%第三題
function flag=mat3(x)
flag=1;
for i=2:sqrt(x)
if rem(x,i)==0
flag=0;
break;
end
end
%在命令窗口調用該函數文件：
for i=10:99
j=10*rem(i,10)+fix(i/10);
if mat3(i)&mat3(j)
disp(i)
end
end
%第四題
function y=fx(x)
y=1./((x-2).^2+0.1)+1./((x-3).^4+0.01);
%在命令窗口調用該函數文件：
y=fx(2)
a=[1,2;3,4];
y=fx(a)
%第五題
%(1)
function f1=mat5(n)
f1=n+10*log(n*n+5);
%在命令窗口中調用該函數文件：
y=mat5(40)/(mat5(30)+mat5(20))
%(2)方法一
function f2=mat6(n)
f2=0;
for i=1:n
f2=f2+i*(i+1);
end
%在命令窗口中調用該函數文件如：
y=mat6(40)/(mat6(30)+mat6(20))
%(2)方法二
function f2=mat7(n)
i=1:n;
m=i.*(i+1);
f2=sum(m);
end
%在命令窗口中調用該函數文件如：
y=mat7(40)/(mat7(30)+mat7(20))

%實驗六 高層繪圖操作
%第一題：
x=linspace(0,2*pi,101);
y=(0.5+3*sin(x)./(1+x.^2)).*cos(x);
plot(x,y)
%第二題：
%（1）
x=linspace(-2*pi,2*pi,100);
y1=x.^2;
y2=cos(2*x);
y3=y1.*y2;
plot(x,y1,'b-',x,y2,'r:',x,y3,'y--');
text(4,16,'\leftarrow y1=x^2');
text(6*pi/4,-1,'\downarrow y2=cos(2*x)');
text(-1.5*pi,-2.25*pi*pi,'\uparrow y3=y1*y2');
%（2）
x=linspace(-2*pi,2*pi,100);
y1=x.^2;
y2=cos(2*x);
y3=y1.*y2;
subplot(1,3,1);%分區
plot(x,y1);
title('y1=x^2');%設置標題
subplot(1,3,2);
plot(x,y2);
title('y2=cos(2*x)');
subplot(1,3,3);
plot(x,y3);
title('y3=x^2*cos(2*x)');
%（3）
x=linspace(-2*pi,2*pi,20);
y1=x.^2;
subplot(2,2,1);%分區
bar(x,y1);
title('y1=x^2的條形圖');%設置標題
subplot(2,2,2);
stairs(x,y1);
title('y1=x^2的階梯圖');
subplot(2,2,3);
stem(x,y1);
title('y1=x^2的桿圖');
subplot(2,2,4);
fill(x,y1,'r');%如果少了'r'則會出錯
title('y1=x^2的填充圖');
%其他的函數照樣做。
%第三題
x=-5:0.01:5;
y=[];%起始設y為空向量
for x0=x
if x0<=0 %不能寫成x0=<0
y=[y,(x0+sqrt(pi))/exp(2)]; %將x對應的函數值放到y中
else
y=[y,0.5*log(x0+sqrt(1+x0^2))];
end
end
plot(x,y)
%第四題：
a=input('a=');
b=input('b=');
n=input('n=');
t=-2*pi:0.01:2*pi;
r=a*sin(b+n*t);
polar(t,r)
%第五題
x=linspace(-5,5,21);
y=linspace(0,10,31);
[x,y]=meshgrid(x,y);%在[-5,5]*[0,10]的范圍內生成網格坐標
z=cos(x).*cos(y).*exp(-sqrt(x.^2+y.^2)/4);
subplot(2,1,1);
surf(x,y,z);
subplot(2,1,2);
contour3(x,y,z,50);%其中50為高度的等級數，越大越密
%第六題
ezsurf('cos(s)*cos(t)','cos(s)*sin(t)','sin(s)',[0,0.5*pi,0,1.5*pi]); %利用ezsurf隱函數
shading interp %進行插值著色處理

%實驗七 低層繪圖操作
%第一題
h=figure('MenuBar','figure','color','r','WindowButtonDownFcn','disp(''Left Button Pressed'')')
%第二題
x=-2:0.01:2;
y=x.^2.*exp(2*x);
h=line(x,y);
set(h,'color','r','linestyle',':','linewidth',2)
text(1,exp(2),'y=x^2*exp(2*x)')
%第三題
t=0:0.00001:0.001;
[t,x]=meshgrid(t);
v=10*exp(-0.01*x).*sin(2000*pi*t-0.2*x+pi);
axes('view',[-37.5,30]);
h=surface(t,x,v);
title('v=10*exp(-0.01*x).*sin(2000*pi*t-0.2*x+pi)');
xlabel(Ct'),ylabel('x'),zlabel('v')
%第四題
x=0:0.01:2*pi;
y1=sin(x);
y2=cos(x);
y3=tan(x);
y4=cot(x);
subplot(2,2,1);
plot(x,y1);
subplot(2,2,2);
plot(x,y2);
subplot(2,2,3);
plot(x,y3);
subplot(2,2,4);
plot(x,y4);
%第五題
cylinder(5);
light('Position',[0,1,1]);
material shiny

%實驗八 數據處理與多項式運算
%第一題
%(1)
A=rand(1,30000);
b=mean(A)
std(A,0,2)
%(2)
max(A)
min(A)
%(3)
n=0;
for i=1:30000
if A(i)>0.5
n=n+1;
end
end
p=n/30000
%第二題
%(1)
A=45+51*rand(100,5);
[Y,U]=max(A)
[a,b]=min(A)
%(2)
m=mean(A)
s=std(A)
%(3)
sum(A,2)
[Y,U]=max(ans)
[a,b]=min(ans)
%(4)
[zcj,xsxh]=sort(ans)
%第三題
h=6:2:18;
x=6.5:2:17.5;
t1=[18,20,22,25,30,28,24];
t2=[15,19,24,28,34,32,30];
T1=spline(h,t1,x)
T2=spline(h,t2,x)
%第四題
x=1:0.1:101;
y1=log10(x);
p=polyfit(x,y1,5)
y2=polyval(p,x);
plot(x,y1,':',x,y2,'-')
%第五題
%(1)
p1=[1,2,4,0,5];
p2=[1,2];
p3=[1,2,3];
p=p1+[0,conv(p2,p3)] %為使兩向量大小相同，所以補0
%(2)
A=roots(p)
%(3)
A=[-1,1.2,-1.4;0.75,2,3.5;0,5,2.5];
polyval(p,A)
%(4)
polyvalm(p,A)

%實驗十五 綜合實驗
%第一題
%定義函數文件
function zp=fs(z,n)
zp=0;
z=0;
for k=1:n;
z=1+imag(z)-1.4*real(z)^2+i*real(z)*0.3;
zp=[zp,z];
end
%在主窗口調用該函數：
n=30000;
z=0;
zp=fs(z,n);
plot(zp,'p','markersize',2,'color','r')
%第二題
quad('1/sqrt(2*pi)*exp(-x.*x/2)',0,1)
syms x
int(1/sqrt(2*pi)*exp(-x.*x/2),0,1)
%第三題
%定義函數文件
function dx=apollo(t,x)
mu=1/82.45;
mustar=1-mu;
r1=sqrt((x(1)+mu)^2+x(3)^2);
r2=sqrt((x(1)-mustar)^2+x(3)^2);
dx=[x(2);2*x(4)+x(1)-mustar*(x(1)+mu)/r1^3-mu*(x(1)-mustar)/r2^3;
x(4);-2*x(2)+x(3)-mustar*x(3)/r1^3-mu*x(3)/r2^3];
%在主窗口調用該函數：
x0=[1.2;0;0;-1.04935751];
options=odeset('reltol',1e-8);
%該命令的另一種寫法是options=odeset;options.reltol=1e-8;
[t,y]=ode45(@apollo,[0,20],x0,options);
plot(y(:,1),y(:,3))%繪制x和y圖像，也就是衛星的軌跡
title('Appollo衛星運動軌跡')
xlabel('X')
ylabel('Y')
%第四題
%這是一個力矩平衡問題，可列方程：500×1.5×cosα=400×2.0×cos(π/3-α)
%解得α＝26.8021°。
%程序一：迭代法解方程
x1=input('x1=?');
x=pi/3.0-acos(15.0*cos(x1)/16);
while abs(x-x1)>=10e-8
x1=x;
x=pi/3.0-acos(15.0*cos(x1)/16);
if x>pi/3.0
disp('error');
end
end
y=x1*180/pi
%程序二：利用繪圖功能
%分別繪制兩個小孩所產生力矩隨α變化的曲線，兩曲線的交點即是蹺蹺板平衡時的α
alpha=linspace(0,pi/2,80000);
m1=500*1.5*cos(alpha); %第一個小孩產生的力矩
m2=400*2*cos(pi/3-alpha); %第二個小孩產生的力矩
k=find(abs(m1-m2)<1e-2); %找平衡點
alpha0=alpha(k);
m0=500*1.5*cos(alpha0);
plot(alpha,m1,alpha,m2,alpha0,m0,'kp');
alpha=alpha0*180/pi %化為角度
%第五題
clear; close;
fplot('18-t^(2/3)',[0,20]);grid on;hold on;
fplot('5+t+2*t^(2/3)',[0,20],'r');hold off;
%發現t約為4
[t,f,h]=fsolve('18-x^(2/3)-5-x-2*x^(2/3)',4)
%求得t=4.6465
t=linspace(0,t,100); y=18-t.^(2/3)-5-t-2*t.^(2/3);
trapz(t,y)-20
%最大利潤6.3232(百萬元)