MATLAB程序设计与应用课后习题答案.pdf

文档介绍：该【MATLAB程序设计与应用课后习题答案 】是由【cby201601】上传分享，文档一共【21】页，该文档可以免费在线阅读，需要了解更多关于【MATLAB程序设计与应用课后习题答案 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。:.


MATLAB

 
_______________

________________
:.

.!":
CommandWindow
'---f
A=[,,-;,,;,,]
A=B=
B=[l,,-;,,;,-,]
-C=A+*B
D=A-B+eye
E=A*B
F=A.*B
:A+*B>A-B+******@ABICDEFGHG=AA
H=
A*B>A.*BI=A/B
I
J=B\A
K=[A,B]
()AM>M
L=[A([,],:);B*]
NA/BPB\A
()[A,B]>[A([,)UBV]
A=B=C=
---
-
D=
E=F=
-
b-
-
-
G=H=
-
I=J
.-...-..
.-..-..-.
.-..-..-.:.
L=
K=-
--
--
---'
-
.WFGA>BAX,B=-
__
()YZ[\C]
()^FGC_`XaFGbcD]
()deMATLABfghijklm
CommandWindow
A=[,J,,;,,,,;,,,,;,,,,;,,,,]
B=[,,;,-,;,,-;,,;,,]
C=A*B
D=C(end-:end,end-:end)
A=B=
-
-
()
D=
r

():.
ArrayEditorA
BQX)IT|Mack:Base|
I
EE
U
EE
U
ArrayEditorB|T!^ArrayEditor-C
,
|t.|ShBQX
®ia

-
-
ArrayEditor-D
snpqr
oisM
L
jJL
I
I
.stug
()),wxiyz{|}|]
()~}€‚ƒ„…†{AB‡€ˆ]
()I
CommandWindowCommandWindow
count=;ch='sldjADAkdAijfGEIdk,;
fori=:subch=ch(find(ch>>a'&ch<'z'))
ifremdj)==
count=count+l;subch=
end
endsldjkdijfdk
count
count=
:.
X
.~}XFG…Y‰nŠn‹>Œ|]
‰Ž
CommandWindow
A=rand()
a=det(A)
r=rank(A)
t=trace(A)
al=norm(&)
a=norm(A)
ainf=norm(Ainf)
A=
.....
.....
.....
.....
.....
AFG‰nŠn‹‘’:
a=r=t=
-..
Œ|’
al=a=ainf=
...
-
.!"
-
A“”P“”ƒ„…•–A|—˜]
‰Ž:.
CommandWindow
CommandWindow
™A=[-,,;,,;-,,]-
[V,D]=eig(A,'nobalance…)
-
CommandWindowCommandWindow
V=D=
...-.
-.-..
-.
...
.
-///--.'
.¦§}¨©ª«///=.
///._
()ª¬;
()^ª_­ƒ„®¯b°C.,±¬,•²³C´µ>¬¶·
´µ;
() ¸|FGA¹º|•–»¼]
CommandWindow
™A=[l/,/,/;/,/,/;/,/,/];
b=[,,.,.]\
bl=[.,.,.]\
X=inv(A)*bX=XI=
Xl=inv(A)*bl
cl=cond(A)..
c=cond(A).-.
..
c=cond(A,inf)
()b´…A¬B…¶·À´µÁb¬:x´,x´Ã…x´]
()
c=c=
.e+.e+.e+
ÄÅAFG¹º|Æ…ÇȨ©ª«Bb´ɅxÊ^ËÌÆ´:.
µ…Í|ÎÏ©³Ð]
Ñ
.ÒÓÔÕÖff× ªØ’:
()fgÉ|ÙÚÃÉÛ,ÙÚܐÝË%;
()fgÉ|ÞÅÃÉۅßË®U
()AàáâÃÉ® Ë]
ãäáåæf>çÖffÉ| èËf×]
CommandWindow
CommandWindow
m=l;
whilent(()
x=input('(!)"');$%
ifx==
break
end
t=input('$%"’)
ift>
y=*+(t-)**(+.);
elseift<
y=*t-;
else
y=*t;y=
end
-
end
'(()'()()"
'$%"'$%"
x=
y=
é
.êë------=----------,-----,-------->>>ìíî]Èn‘ïnn
:.
Ʌ»ð§ñòóô‘kõö»÷>ƒ„øù]
ƒ„:
-n=|;-n=|;-n=|;
-i=l:n;-i=l:n;-i=:n;
-f=l./i.*;-f=l./i.*;-f=l./i.*;
-y=sum(f);-y=sum(f);-y=sum(f);
-sqrt(*y)-sqrt(*y)-sqrt(*y)
ans=ans=ans=
....
V_a°!"a#b%&'()*
Xn+}h+
x
()+,-./''%Xo=LO
=)>?@=*
()A@-BCDr,FGr'HIJK±JN+a,PQ
b)'R(,)#(,)#(,.)TUVWXHIYZ
[\*
CommandWindowCommandWindowCommandWindow
>>x=;a=;b=l;lx=l;a=;b=;lx=lfa=;b=.;
fori=:fori=:fori=l:
xl=a/(b-hx)fxl=a/(b+x);xl=a/(b-+x);
ifabs(xgx)<=*(g)imabs(xgx)<="(g)ifabs(xgx)<="(-)
breakbreakbreak
endendend
x=x;x=x;x=x;
endendend
xlxlxl
xl=xl=xl=
#...:.
*
CommandWindowCommandWindowCommandWindow
a=l;b=l;a=;b=;a=;b=.;
r=(-b+sqrt(b*+*a))/r=(-b+sqrt(b*+*a))/r=(-b+sqrt(b-+*a))/
r=(-b-sqrt(b*+*a))/r=(-b-sqrt(b-+*a))/r=(-b-sqrt(b"+*a))/
rl=r=rl=
...
r=r=r=
-.-.-.
.pqrstuv)'wxyJz){|}qrstuv)J~)W,
€z)J~z)*AXB=Jz)‚XJ~)W
J~z)*/,
ƒ„…
()~)W'W)*
()†‡ˆ~)WW‰'‚Š~z)‹X*
-k=of
-sm=;
-p=:;
-fori=:
-n=find(rem(p,i)=+p~=i);
-p(n)=[];
-end
-P
-fori=:
-m=i*;k=
-if(find(m==p)>O&find(m==p)<=ength(p))
-k=k+l;
-sum=sum+p(find(m==p));
-end
-endsum=
-k
-sum
–—˜
™#–—…š:.
.›/(")=_+,KrMATLABœ)
'-(X-)+.(X-)+.
,žŸ ¡/xT,x¢¡£¤¥Ÿ¦'/§X¨%©ª£¤*
functionf=fx(x)
fori=:length(x(:end,:))
forj=:length(x(::end))
f(i,j)=/((x(i,j)-).*+.)+/((x(i,j)-).A+.);
end
end|
lA=[l,,;,,;,,];
y=fx(A)
y=
...
...
...
.¯°()
y
±)²±)
³P/(´)=´+l°ln(´+)T,/y'*
¦PT/§´¨=*+*+*+K+µ*§µ+,/y'*
()
functionf=fl(n)ly=fl()/(fl()+fl())
-f=n+.*log(n.*+);
y=
.
*
CommandWindow
|f\mctionf=f(n)
>y=f()/(f()+fy())
-f=f(n-l)+f(n);
y=
functionf=f(n)
-f=n*(n+);
.
–—·:.
›#=.sinx,¸x=
¹ƒ„Rº»¼œ)
COSX
+x
'½¾*
x=linspace(,*pi,);
plot(x,(.+*sin(x)/(+x.*))*cos(x))
¿»¼ÀÁ½¾p
asin(Z?+*<)
ÃUÄÅ)a#b#nW½¾ÆÇ'ÈÉ*
CommandWindow
a=;b=;n=;
t=linspace(,*pi,);
p=a*sin(b+n*t);
polarCt,p)
:.
‡˜ÊËÌUVJa=l,b=,n=l####T'ÍÎ>ÏÐÑPn
%Ò)TÓ¦'̊Ò)rÔÕPn%Ö)TÓ¦'̊€Ö)'q×r
Ô*Å)aؼÀÁÂ'ÙÚÅ)b¢WÌYZÛÜÝÞ*
.»¼½ßÌÆÃYZàáâãä
x=cosScost
y=cosSsintFf*
{z=sinS
-s=linspace(jpi/j);
-t,linspace(j*pi/j);
-[S,Tl^eshgridtsjt);
-x=cos(S).*cos(I);
-y=cos(S).*sin(I);
-z=sin(S);
-surf(x,y,z);
-shadinginterp;
-xlabel(/x-axis/)ylabel(/-axis'))zhbel('z-axis/);
J
-titlefsurf,)|
–—å
.桽ßW综½ßèé)=iOe-°exsin(+):.
êæ¡ëìí½¾vïð@ÌÆñòóôõö÷½¾'øâ#¾
ùX¾úÃ桝ûWçü½¾ýþûÂÿ
-Ix=linspace(,,);
-t=linspace(,,);
-[XjT]=meshgrid(x,t);
-v=*exp(-.*X).*sin(*pi*T-.*X+pi);
-axes('view',[-.,]);
-hs=surface(x,t,v,JFaceColor,,,w*,*EdgeColor,*flat*);
-gridon;
-pause;
-set(he,JFaceColor,,*flaty;
.:MATLAB;<=rand>?***@ABCDEFGH=IJ?,
KLMIJ?=NOP
()EQRSTUV
()WXYZRW[YZ
()X\.=IJ?B?]^?=_G`
-~|ft=rand(,);
-M=mean(A)
-D=std(A)
-E=max(A)
-F=min(A)
-N=length(find(A>.))/:.
()i()
.jklmnonpjq:~:stuvh=wxyz{|}~
€‚
€wxyz{mnƒ„|}~
…th
wxz{tl.......
wyz{†.......
‡:ˆ‰Š‹ŒQGŽqwxy:-:stuvh‘
’=“”z{|°C~
-h=::;
-hl=.::.;
-=[.,.,.,.,.,.,.];
-t=[.,.,.,.,.,.,.];
-=spline(h,hl)
-t=spline(h,t,hl)
CommandWindow
|t=
ColumsthroughColonsthrough
,.......
ColwnnsthroughColumnsthrough
....
.–B—˜™š›~=œ+/+/+|~=+,…:.
žŸ ¡¢:
()ŽP(X)=R(X)+š(X)š@)
()ŽP(x)=¦
()§x¨©ªA=u«YZ…ŽP(x)=Q¬­:
-.-.
A=..
.
()§®©ªA¯°±²…ŽP(x)=Q¬­A=Q³´()
µ¶·
-pl=[l,,,,];
-P=[l,];
-p=[l,,];
-g=[,conv(p,p)];
-p=pl+g
-x=roots(p)
-A=[-l,.,-.;.,,.;,,.];
-xl=polyval()
-x=polyvalm(p,A)|
()()
CommandWindowx=
-.+.i
P=-.-.
-.+.i
-.-.
()()
xl=x=
.e+*.e+*
....-.-.
......
......
¸
.Ž>?¹º»’=?Q¼?:.
XXX
fM=xV,x=,,
x
.:?QU½Ž»¾G
(i)/=f"^cos^sin(r^^=“”Q
-t=:.:*pi;ans=
I=sqrt(cos()+*sin((*t).*)+);
-|trapz(t,I).+.i
.G:¿À·=?QU½ÁÂNUÃÄ
x+y-z+z=-
x=y+zÆÇ=
x+y+z-Ç=
x-y+Ç=ll
ÈÉÁ½:
.
-A=[,,-,;,-,,-;,,,-;,-,,];-.
-b=[-;;;];.
-k=A'b-.
LUGÁ:
«A=[,,-,;,-,,-;,%,-;,-,,];.
-b=[-;;l;ll];-.
_|[UU]=lu(A);.
«x=U(L\b)-.
ÌÁ½:
ÍUÃĖΫÁx
.
.A=[,,-,;,-,,-;,%,-;,-,,];-.
-b=[-;;l;ll];.
-x=ine_solution(A,b)|
-.:.
.STUVWXYZ[\]^
Ï+%+%+=
,Ï+%+E+X=
,x+%+x+x=
function[x,y]=line_solution(A,b)
-[m,n]=size(A);
-y=[];
-ifnorm(b)>
-ifrank(A)==rank([A,b])
-ifrank(A)==n
-disp(,_YZ[`ab^x');
-x=A\b;
-else
-disp('_YZ[`cde^fg^hx,iiVYZ[\jk^lhy');
-x=A\b;
-y=null(A/r);
-end
-else
-disp('YZ[c^’)
-x=[];
-end
-else
-disp('_YZ[`m^x');
-x=zeros(n$);
-ifrank(A)<n
-disp('YZ[`cde^fjk^lhy');
-y=null(A/r,);
-end
-end
ÖA=[,,,;,,,;,,,];
b=[;;];
[x,y]=line_solution(A,b)
_YZ[`ce^fg^hxiiVYZ[\jk^lhy
Warning:Rankdeficient,rank=,tol=.e-.
>Inline_solutionat
X=y=
-..-.
.-..
.
.
.Sn%YZ\%o^
()pqr\sot=,%=,&=ufSYZ[\%o^:.
sinx+y+lnx-=
<x+V-^+l=
x+y+z-=
functionF=myfun(X)
-x=X(l);
-y=X(l);
-z=x|(l);
-F(l)=sin(x)+y.*+log(z)-;
-F()=*x+.Vz.*+Ú
-F()=xy+z-;
ÖX=fsolve('myfun,[,,],optimset('Display'Joff'))
X=
...
.Ž>?¹º»Ût=ÜQ
,
/ÝX+COSX+XlogX/Ý
()/)=à--------ÆdÆJ¹(,)x=W[Q
e
functionf=f(x)
-f=(+cos(x)+x*iog(x))/exp(x);|
x=fmin=
|>>[x,fmin]=fminbnd(?fx?,,)j
.ŽâGUÃ=?QÁ
^%+y=O
clxclx
-y(O)=O
y'(o)=o
.ŽâGUÃÄ=?QÁãäåÁ=æÂ
Y=EX
Ý=ÆAN
è=-«è
y,(O)=O,y(O)=Ly<O)=l:.
functiondy=rigid(t,y)
-dy=zerosé);
-dy(l)=y()*y();
-dyi=-y(l)*yéÚÖ[T,Y]=odeCrigid\[,],[,LJ)Ú
-dy()=-.*y(l)*y();|plot(T,Y(:,)JT,Y(:$),'Y(:$)/+')
ê
.ëìx=,y=,:Cí€îŽz=ïðñ
-x=sym(');
z=
-y=sym(');
-z=(x+l)/(sqrt(-h<)-sqrt(y))|/O-*(/))
.GÁó™
x-y
-symsxy;ans=
-z=x.*-y.*;
-factor(z)I(x-y)*(x+y)*(x*+y*)
.ôõ€î™
()sincos-cosösin
-symsxy;ans=
-z=sin(x)*cos(y)-cos(x)*sin(y);
-simple(z)|sin(x-y)
.ëì:.
bc
ef
hi
÷AŸ øù:
B=ppjA°()B=û©ªãüƒ„
()ýþB©ªÿ

 ()B
.
!"
-symsx;ans=
-f=x*(exp(sin(x)+)-*(exp(tan(x)-)))/(sin(x)).*;
-limit(f)|Inf
.:;
rdx
<I------------------==
J(arcsinx)V+x
-x=sym(=x');ans=
-f=/((asin(x)).**sqrt(-x.*));
-int(f)|-/asin(x)
BC
w