MATLAB数学实验第二版答案(胡良剑)(完整教资)
系统教学#
1
数学实验答案
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
1.2500
Page20,ex5
>> 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函数，保存为eg2_1.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
MATLAB程序设计教程(第二版)课后答案(可编辑修改word版)
MATLAB 第二版课后答案unit3-8 unit3实验指导1、 n=input('请输入一个三位数：');a=fix(n/100);b=fix((n-a*100)/10);c=n-a*100-b*10;d=c*100+b*10+a2（1）n=input('请输入成绩');switch ncase num2cell(90:100)p='A';case num2cell(80:89)p='B';case num2cell(70:79)p='C';case num2cell(60:69)p='D';otherwisep='E';endprice=p（2）n=input('请输入成绩');if n>=90&n<=100p='A';elseif n>=80&n<=89p='B';elseif n>=70&n<=79p='C';elseif n>=60&n<=69p='D';elsep='E';endprice=p（3）tryn;catchprice='erroe'end3n=[1,5,56,4,3,476,45,6,3,76,45,6,4,3,6,4,23,76,908,6];b=n(1);for m=2:20if n(m)>aa=n(m);elseif n(m)<bb=n(m);endendmax=amin=b法2n=[1,5,56,4,3,476,45,6,3,76,45,6,4,3,6,4,23,76,908,6];min=min(n)max=max(n)4b=[-3.0:0.1:3.0];for n=1:61a=b(n);y(n)=(exp(0.3*a)-exp(-0.3*a))/2*sin(a+0.3)+log((0.3+a)/2);endy5y1=0;y2=1;n=input('请输入n 的值：');for i=1:ny1=y1+1/i^2;y2=y2*((4*i*i)/((2*i-1)*(2*i+1)));endy1y26A=[1,1,1,1,1,1;2,2,2,2,2,2;3,3,3,3,3,3;4,4,4,4,4,4;5,5,5,5,5,5;6,6,6,6,6,6]; n=input('请输入n 的值:');if n<=5&n>=0disp(A([n],:));elseif n<0disp(lasterr);else disp(A([6],:));disp(lasterr);end7（1）f=[];f(n)=n+10*log(n^2+5);endy=f(40)/(f(30)+f(20))（2）f=[];a=0;for n=1:40f(n)=a+n*(n+1);a=f(n);endy=f(40)/(f(30)+f(20))8y=0;m=input('输入m 的值:');n=input('输入n 值：');for i=1:ny=y+i^m;endy************************************************************ function s=shi8_1(n,m)s=0;for i=1:ns=s+i^m;end************************************************************ shi8_1(100,1)+shi8_1(50,2)+shi8_1(10,1/2)思考练习2N=[1,2,3,4,5];2.*NN./21./N1./N.^23s=fix(100*rand(1,20)*9/10+10)y=sum(s)/20j=0;for i=1:20if s(i)<y&rem(s(i),2)==0j=j+1;A(j)=s(i);else continue;endendA4y1=0;y2=0;n=input('请输入n 的值：'); for i=1:ny1=y1+-(-1)^i/(2*i-1);y2=y2+1/4^i;endy1y2unit4实验指导1（1）x=-10:0.05:10;y=x-x.^3./6;plot(x,y)（2）x=-10:0.5:10;ezplot('x^2+2*y^2-64',[-8,8]); grid on;2t=-pi:pi/10:pi;y=1./(1+exp(-t));subplot(2,2,1);bar(t,y);title('条形图(t,y)');axis([-pi,pi,0,1]);subplot(2,2,2);stairs(t,y,'b');title('阶梯图(t,y)');axis([-pi,pi,0,1]);subplot(2,2,3);stem(t,y,'k');title('杆图(t,y)');axis([-pi,pi,0,1]);subplot(2,2,4);loglog(t,y,'y');title('对数坐标图(t,y)');3（1）t=0:pi/50:2*pi;r=5.*cos(t)+4;polar(t,r);title('\rho=5*cos\theta+4'); （2）t=-pi/3:pi/50:pi/3;r=5.*((sin(t)).^2)./cos(t); polar(t,r);4（1）t=0:pi/50:2*pi;x=exp(-t./20).*cos(t);y=exp(-t./20).*sin(t);z=t;plot3(x,y,z);grid on;（2）[x,y]=meshgrid(-5:5);z=zeros(11)+5;mesh(x,y,z);shading interp;5[x,y,z]=sphere(20);surf(x,y,z);axis off;shading interp;m=moviein(20);for i=1:20axis([-i,i,-i,i,-i,i])m(:,i)=getframe;endmovie(m,4);思考练习2（1）x=-5:0.1:5;y=(1./(2*pi)).*exp((-(x.^2))/2); plot(x,y);（2）t=-2*pi:0.1:2*pi;x=t.*sin(t);y=t.*cos(t);plot(x,y);grid on;3t=0:pi/1000:pi;x=sin(3.*t).*cos(t);y1=sin(3.*t).*sin(t);y2=2.*x-0.5;plot(x,y1,'k',x,y2);hold on;k=find(abs(y1-y2)<1e-2);x1=x(k);y3=2.*x1-0.5;plot(x1,y3,'rp');4x=-2:0.01:2;y=sin(1./x);subplot(2,1,1);plot(x,y);subplot(2,1,2);fplot('sin(1./x)',[-2,2],1e-4);5（1）i=-4*pi:0.1:10;j=12./sqrt(i);polar(i,j);title('{\rho}=12/sqrt(\theta)')（2）a=-pi/6:0.01:pi/6;b=3.*asin(a).*cos(a)./((sin(a)).^3+(cos(a)).^3); polar(a,b);6（1）[u,v]=meshgrid(-4:0.1:4);x=3.*u.*sin(v);y=2.*u.*cos(v);z=4.*u.^2;subplot(2,1,1);mesh(x,y,z);subplot(2,1,2);surf(x,y,z);（2）[x,y]=meshgrid(-3:0.2:3);z=-5./(1+x.^2+y.^2);subplot(1,2,1);mesh(x,y,z);subplot(1,2,2);surf(x,y,z);unit5实验指导1A=randn(10,5)x=mean(A)y=std(A)Max=max(max(A))Min=min(min(A))Sumhang=sum(A,2)SumA=sum(Sumhang)B=sort(A);C=sort(B,2,'descend');C2（1）a=0:15:90;b=a./180.*pi;s=sin(b)c=0:15:75;d=c./180.*pi;t=tan(d)e=input('请输入想计算的值：'); S=sin(e/180*pi)T=tan(e/180*pi)S1=interp1(a,s,e,'spline')T1=interp1(c,t,e,'spline')P1=polyfit(a,s,5);P2=polyfit(c,t,5);S2=polyval(P1,e)T2=polyval(P2,e)（2）n=[1,9,16,25,36,49,64,81,100]; N=sqrt(n);x=input('ji suan zhi : ');interp1(n,N,x,'cubic')3N=64;T=5;t=linspace(0,T,N);h=exp(-t);dt=t(2)-t(1);f=1/dt;X=fft(t);F=X(1:N/2+1);f=f*(0:N/2)/N;plot(f,abs(F),'-*')4P=[2,-3,0,5,13];Q=[1,5,8];p=polyder(P)q=polyder(P,Q)[a,b]=polyder(P,Q)5P1=[1,2,4,0,5];P2=[0,1,2];P3=[1,2,3];P=P1+conv(P2,P3)X=roots(P)A=[-1,1.2,-1.4;0.75,2,3.5;0,5,2.5];p=polyval(P,A)思考练习4A=rand(1,30000);a=mean(A)b=std(A)Max=max(A)Min=min(A)n=0;for i=1:30000if(A(i)>0.5)n=n+1;endendny=n/300005p=[45,74,54,55,14;78,98,45,74,12;87,98,85,52,65][M,S]=max(p)[N,H]=min(p)junzhi=mean(p,1)fangcha=std(p,1,1)zong=sum(p,2);[Max,wei]=max(zong)[Min,wei]=min(zong)[zcj,xsxh]=sort(zong,'descend')6x=[1:10:101];y=[0,1.0414,1.3222,1.4914,1.6128,1.7076,1.7853,1.8513,1.9085,1.9590,2.0043]; [p,s]=polyfit(x,y,5)a=1:5:101;y1=polyval(p,a);plot(x,y,':o',a,y1,'-*')unit6实验指导1A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6];p=[0.95,0.67,0.52]';x=A\pA=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; p=[0.95,0.67,0.53]';x=A\pcond(A)2（1）x1=fzero(@funx1,-1)function fx=funx1(x)fx=x^41+x^3+1;（2）x2=fzero(@funx2,0.5)function fx=funx2(x)fx=x-sin(x)/x;（3）options=optimset('Display','off');x=fsolve(@fun3,[1,1,1]',options)q=fun3(x)function q=fun3(p)x=p(1);y=p(2);z=p(3);q(1)=sin(x)+y^2+log(z)-7;q(2)=3*x+2^y-z^3+1;q(3)=x+y+z-5;3（1）t0=0;tf=5;y0=1;[t,y]=ode23(@fun4,[t0,tf],y0);t'y'function yp=fun4(t,y)yp=-(1.2+sin(10*t))*y;（2）t0=0;tf=5;y0=1;[t,y]=ode23(@fun5,[t0,tf],y0);t'y'function yp=fun5(t,y)yp=cos(t)-y/(1+t^2);4x=fminbnd(@mymin,0,2);-mymin(x)function fx=mymin(x)fx=-(1+x.^2)/(1+x.^4);5options=optimset('Display','off');[x,fval]=fmincon(@fun6,[0,0,0],[],[],a,b,lb,ub)-fvalfunction f=fun6(x)f=-(sqrt(x(1))+(400-x(1))*1.1+(sqrt(x(2))+(400-x(1))*1.1-x(2))*1.1+sqrt(3)+(((400- x(1))*1.1-x(2))*1.1-x(3))*1.1+sqrt(x(x4)));思考练习1（1）A=[2,3,5;3,7,4;1,-7,1];B=[10,3,5]';C1=inv(A)*BC2=A\B[L,U]=lu(A);x=U\(L\B)（2）A=[6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2];B=[-4,13,1,11]';C1=inv(A)*BC2=A\B[L,U]=lu(A);x=U\(L\B)2（1）x1=fzero(@funx1,1.5)function fx=funx1(x)fx=3*x+sin(x)-exp(x);（2）x1=fzero(@funx2,1)function fx=funx2(x)fx=x-1/x+5;（3）options=optimset('Display','off');x=fsolve(@fun3,[3,0]',options)q=fun3(x)function q=fun3(p)x=p(1);y=p(2);q(1)=x^2+y^2-9;q(2)=x+y-1;3（1）t0=0;tf=5;y0=[0,1];[t,y]=ode45(@vdpol,[t0,tf],y0);[t,y]function ydot=vdpol(t,y);ydot(1)=(2-3*y(2)-2*t*y(1))./(1+t^2);ydot(2)=y(1);ydot=ydot';（2）t0=0;tf=5;y0=[1;0;2];[t,y]=ode45(@vdpoll,[t0,tf],y0);[t,y]function ydot=vdpoll(t,y);ydot(1)=cos(t)-y(3)./(3+sin(t))+5*y(1).*cos(2*t)/((t+1).^2)-y(2); ydot(2)=y(1);ydot(3)=y(2);ydot=ydot';4x=fminbnd(@mymin,0,pi);-mymin(x)function fx=mymin(x)fx=-sin(x)-cos(x.^2);5[x,y1]=fminbnd(@mymax,0,1.5);-y1function fx=mymax(x);fx=-(9*x+4*x.^3-12*x.^2);unit7实验指导1（1）format longfx=inline('sin(x)./x');[I,n]=quadl(fx,0,2,1e-10)（2）format longfx=inline('1./((x-0.3).^2+0.01)-1./((x-0.9).^2+0.04)-6');[I,n]=quad(fx,0,1,1e-10)2（1）global ki;ki=0;I=dblquad(@fxy,0,1,0,1)ki（2）f=inline('abs(cos(x+y))','x','y');I=dblquad(f,0,pi,0,pi)3X=0.3:0.2:1.5;F=[0.3895,0.6598,0.9147,1.1611,1.3971,1.6212,1.8325];trapz(X,F)4p=0:pi/5:2*pi;for n=1:3nDX=diff(sin(p),n)end5f=inline('sin(x)./(x+cos(2.*x))');g=inline('(cos(x).*(x+cos(2*x))-sin(x).*(1-2.*sin(2*x)))/(x+cos(2.*x)).^2');x=-pi:0.01:pi;p=polyfit(x,f(x),5);dp=polyder(p);dpx=polyval(dp,x); %求dp 在假设点的函数值dx=diff(f([x,3.01]))/0.01; %直接对f(x)求数值导数gx=g(x); %求函数f 的导函数g 在假设点的导数plot(x,dpx,x,dx,'.',x,gx,'-'); %作图思考练习2format longfx=inline('1./(1+x.^2)');[I,n]=quad(fx,-Inf,Inf,1e-10)[I,n]=quadl(fx,-Inf,Inf,1e-10)x=-100000:100000;y=1./(1+x.^2);trapz(x,y)format short3(1)format longfx=inline('log(1+x)./(1+x.^2)');[I,n]=quad(fx,0,1,1e-10)(2)format longfx=inline('sqrt(cos(t.^2)+4*sin((2*t).^2)+1)'); [I,n]=quad(fx,0,2*pi,1e-10)4f=inline('4.*x.*z.*exp(-z.^2.*y-x.^2)');I=triplequad(f,0,pi,0,pi,0,1)5f=inline('sin(x)');g=inline('cos(x)');x=0:0.01:2*pi;p=polyfit(x,f(x),5);dp=polyder(p);dpx=polyval(dp,x);dx=diff(f([x,2*pi+0.01]))/0.01;gx=g(x);plot(x,dpx,x,dx,'.',x,gx,'-')unit8实验指导1syms x y;s=x^4-y^4;factor(s)factor(5135)2syms x;f=(x-2)/(x^2-4);limit(f,x,2)sym x;f=(sqrt(pi)-sqrt(acos(x)))/sqrt(x+1);limit(f,x,-1,'right')3sym x;f=sin(1/x);diff(f,'x')diff(f,'x',2)sym x;f=(1-cos(2*x))/x;diff(f,'x')diff(f,'x',2)4sym x;f=sqrt(exp(x)+1);int(f,'x')syms x y;f=x/(x+y);int(f,'y')sym x;f=exp(x)*(1+exp(x))^2;int(f,'x',0,log(2))sym x;f=x*log(x);int(f,'x',1,exp(1))5sym x;s=symsum((-1)^(x+1)/x,1,Inf)sym y;z=symsum(y^(2*y-1)/(2*y-1),1,Inf)6sym x;f1=(exp(x)+exp(-x))/2;f2=sqrt(x^3-2*x+1);taylor(f1,x,5,0)taylor(f2,x,6,0)7syms x y a;x=solve('x^3+a*x+1=0','x')x=solve('sin(x)+2*cos(x)-sqrt(x)=0','x')[x y]=solve('log(x/y)=9','exp(x+y)=3','x','y')8syms n;[x,y]=dsolve('x*(D2y)+(1-n)*(Dy)+y=0','y(0)=0','Dy(0)=0','x') 思考练习2syms x B1 B2 a bs1=2*((cos(x))^2)-(sin(x))^2;s2=sin(B1)*cos(B2)-cos(B1)*sin(B2);s3=sqrt((a+sqrt(a^2-b))/2)+sqrt((a-sqrt(a^2-b))/2); s4=(4*x^2+8*x+3)/(2*x+1);h1=simplify(s1)h2=simplify(s2)h3=simplify(s3)h4=simplify(s4)3syms x a;f=abs(x)/x;limit(f,x,0,'left')f=(x+a/x)^x;limit(f,x,inf)4syms x y mf=sqrt(x+sqrt(x+sqrt(x)));m=diff(f,'x')diff(m,'x')syms x y;f=x+y-sqrt(x^2+y^2);z1=diff(f,x)z2=diff(f,y)5syms x y;f=x+y-sqrt(x^2+y^2);z1=diff(f,x)z2=diff(f,y)sym x;f=1/(asin(x)^2*(1-x^2)^(1/2));int(f)6syms xf=1/(1+x);int(f,0,4)sym x;f=x^3*sin(x)^2/(x^6+2*x^4+1);int(f,-1,1)sym x;f=x^3*sin(x)^2/(x^6+2*x^4+1);int(f,-1,1)7syms ns=symsum(1/4^n,1,inf)sym n;s=symsum(((n+1)/n)^(1/2),1,inf)eval(y)8syms xf=tan(x);taylor(f,x,3,0)syms xf=sin(x)^2;taylor(f,x,5,0)9syms xx=solve('log(1+x)-5/(1+sin(x))=2','x')syms x y z[x y z]=solve('4*x^2/(4*x^2+1)=y','4*y^2/(4*y^2+1)=z','4*z^2/(4*z^2+1)=x','x','y','z') 10[x ,y]=dsolve('Dx=3*x+4*y','Dy=5*x-7*y','x(0)=0','y(0)=1','t')。
MATLAB程序设计与应用(第二版)实验答案-095
%在命令窗口调用上述函数文件：
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];
%
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))
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)
%(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)
matlab程序设计与应用第二版习题答案
matlab程序设计与应用第二版习题答案matlab程序设计与应用第二版习题答案【篇一：matlab程序设计与应用(第二版)实验答案】%实验一 matlab运算基础%第1题%（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=0t1)+(t.^2-1).*(t=1t2)+(t.^2-2*t+1).*(t=2t3)%第2题a=[12 34 -4;34 7 87;3 65 7];b=[1 3 -1;2 0 3;3 -2 7];a+6*ba-b+eye(size(a))a*ba.*ba^3a.^3a/bb\a[a,b][a([1,3],:);b^2]%第3题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*bf=size(c)d=c(f(1)-2:f(1),f(2)-1:f(2))whos%第4题%(1):a=100:999;b=rem(a,21);c=length(find(b==0))%(2):a=lsdhksdlkklsdkl;k=find(a=aa=z);a(k)=[]%实验二 matlab矩阵分析与处理 %第1题e=eye(3);r=rand(3,2);o=zeros(2,3);s=diag([2,3]);a=[e,r;o,s];a^2b=[e,(r+r*s);o,s^2]%第2题h=hilb(5)p=pascal(5)hh=det(h)hp=det(p)th=cond(h)tp=cond(p)%第3题a=fix(10*rand(5))h=det(a)trace=trace(a)rank=rank(a)norm=norm(a)%第4题a=[-29,6,18;20,5,12;-8,8,5][v,d]=eig(a)%数学意义略%第5题方法一%(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)%第5题方法二a=hilb(4)a(:,1)=[]a(4,:)=[]b=[0.95,0.67,0.52];x=inv(a)*bb1=[0.95,0.67,0.53];x1=inv(a)*b1n=cond(b)n1=cond(b1)na=cond(a) %矩阵a为病态矩阵%第6题a=[1,4,9;16,25,36;49,64,81]b=sqrtm(a)c=sqrt(a) %sqrtm函数是以矩阵为单位进行计算，sqrt函数是以矩阵中的元素进行计算%实验三选择程序结构设计%第1题程序一x=[-5.0,-3.0,1.0,2.0,2.5,3.0,5.0];y=[]; %建立存放所有y值的矩阵for x0=xif x00x0~=-3y=[y,x0*x0+x0-6];elseif x0=0x05x0~=2x0~=3y=[y,x0*x0-5*x0+6];elsey=[y,x0*x0-x0-1];endendx%输出所有xy%输出所有y%第1题程序二x=[-5,-3,1,2,2.5,3,5];y=[];for a=1:7if x(a)0x(a)~=-3y=[y,(x(a))^2+x(a)-6];elseif x(a)=0x(a)5x(a)~=2x(a)~=3y=[y,(x(a))^2-5*x(a)+6];elsey=[y,x(a)*x(a)-x(a)-1];endend%第2题程序一x=input(请输入一个百分制成绩：);if x100|x0disp(您输入的成绩不是百分制成绩，请重新输入。
MATLAB 程序设计教程第二版实验答案包括第七章
第二章：1.(1) w=sqrt(2)*(1+0.34245*10^(-6))w =1.4142(2) a=3.5;b=5;c=-9.8;x=(2*pi*a+(b+c)/(pi+a*b*c)-exp(2))/(tan(b+c)+a)x =0.9829(3) a=3.32;b=-7.9;y=2*pi*a^2*((1-pi/4)*b-(0.8333-pi/4)*a)y =-128.4271(4) t=[2,1-3i;5,-0.65];>> z=0.5*exp(2*t)*log(t+sqrt(1+t*t))z =1.0e+004 *0.0057 - 0.0007i 0.0049 - 0.0027i1.9884 - 0.3696i 1.7706 - 1.0539i2A=[-1,5,-4;0,7,8;3,61,7];B=[8,3,-1;2,5,3;-3,2,0]; >>(1) A+6*Bans =47 23 -1012 37 26-15 73 7>> A^2-B+eye-18 -216 1823 533 11022 868 526>>(2) A*Bans =14 14 16-10 51 21125 328 180>> A.*Bans =-8 15 40 35 24-9 122 0>> B*Aans =-11 0 -157 228 533 -1 28(3) A/Bans =1.2234 -0.92552.9787 -0.9468 2.3511 -0.9574 4.61703.8723 13.8936>> B\Aans =-0.5106 -8.6170 -1.1277 0.7340 17.5745 1.8085 -0.8830 -21.2128 0.4043>>(4) [A,B]ans =-1 5 -4 8 3 -10 7 8 2 5 33 61 7 -3 2 0>> [A([1,3],:);B^2]ans =-1 5 -43 61 773 37 117 37 13-20 1 9>>3A=[23,10,-0.778,0;41,-45,65,5;32,5,0,32;6,-9.54,54,3.14];(1) b=find(A>=10&A<=25);>> A(b)ans =2310(2) B=A(1:3,:)B =23.0000 10.0000 -0.7780 041.0000 -45.0000 65.0000 5.000032.0000 5.0000 0 32.0000C=A(:,1:2)C =23.0000 10.000041.0000 -45.000032.0000 5.00006.0000 -9.5400>> D=A(2:4,3:4)D =65.0000 5.0000 0 32.0000 54.0000 3.1400>> E=B*CE =1.0e+003 *0.9141 -0.22391.20802.71231.1330 -0.2103(3) E<Dans =0 10 00 1>> E&Dans =1 10 11 1>> E|Dans =1 11 11 1~E|~Dans =0 01 00 0第三章：2．（1）c=input('请输入成绩');if(c>=90&c<=100)disp 'A' ;elseif(c>=80)disp 'B';elseif(c>=70)disp 'C';elseif(c>60)disp 'D';elsedisp 'E';end请输入成绩55Ec=input('please enter:'); switch fix(c/10)case {9,10}disp('A');case {8}disp('B');case {7}disp('C');case {6}disp('D');case {0,1,2,3,4,5} disp('E');end>> asplease enter:55E>> asplease enter:67D（2）while(1)c=input('请输入成绩');if(c>=90&c<=100)disp 'A' ;elseif(c<90&c>=80)disp 'B';elseif(c<80&c>=70)disp 'C';elseif(c<70&c>=60)disp 'D';elseif(c<60&c>=0)disp 'E';elsedisp 'wrong'endend>> third请输入成绩113wrong请输入成绩-47wrong请输入成绩45E请输入成绩89B请输入成绩78C3.a=[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20] b=max(a)c=min(a)b=0c=0for i=1:20if a(i)>bb=a(i);endif a(i)<cc=a(i);endendbc6.A=rand(5,6)n=input('请输入矩阵：');if n>5n=5A(n,n)elseA(n)End6. a=rand(5,6);n=input('请输入n');if n==1disp(a(1,:))endif n==2disp(a(2,:))endif n==3disp(a(3,:))endif n==4disp(a(4,:))endif n==5disp(a(5,:))endif n>5disp(a(5,:))end第四章：1.（1）x=-12:0.001:12;y=x-x.^3/6;plot(x,y);xlabel('x轴');ylabel('y轴');title('Y=X-X^3/3!');（2）ezplot('x^2+2*y^2-64',[-8,8,-8,8])xlabel('x轴');ylabel('y轴');axis equaltitle('x^2+2y^2=64');4.(1)t=0:pi/100:2*pi;x=exp(-t/20).*cos(t);y=exp(-t/20).*sin(t);z=t;plot3(x,y,z)(2)[x,y]=meshgrid(-5:5);z=ones(11);surf(x,y,z);shading interp;第五章1.（1）A=randn(10,5);a=mean(A)a =0.0013 0.2310 0.0253 0.3588 -0.4197b=std(A)b =0.9034 0.8829 1.1898 0.7832 1.0821（2）max(max(A))ans =2.1832>> min(min(A))ans =-2.1707（3）.sum(sum(A))ans =1.9666sum(A,1)ans =0.0128 2.3095 0.2528 3.5884 -4.1969（4）.sort(A)ans =-1.6656 -0.8323 -1.5937 -1.2025 -2.1707-1.1465 -0.5883 -1.4410 -0.3999 -1.6041-0.4326 -0.1867 -1.3362 -0.1567 -1.0565-0.0376 -0.1364 -0.6918 -0.0198 -0.92190.1253 -0.0956 0.2944 0.6686 -0.80510.1746 0.0593 0.5711 0.6900 -0.05920.2877 0.1139 0.7143 0.7119 0.21930.3273 0.7258 0.8580 0.8156 0.25731.1892 1.0668 1.2540 1.1908 0.52871.19092.1832 1.6236 1.2902 1.4151>> sort(ans,2,'descend')ans =-0.8323 -1.2025 -1.5937 -1.6656 -2.1707-0.3999 -0.5883 -1.1465 -1.4410 -1.6041-0.1567 -0.1867 -0.4326 -1.0565 -1.3362-0.0198 -0.0376 -0.1364 -0.6918 -0.92190.6686 0.2944 0.1253 -0.0956 -0.80510.6900 0.5711 0.1746 0.0593 -0.05920.7143 0.7119 0.2877 0.2193 0.11390.8580 0.8156 0.7258 0.3273 0.25731.2540 1.1908 1.1892 1.0668 0.52872.1832 1.6236 1.4151 1.2902 1.19092.（2）n=1:10;a=n.*n;interp1(n,a,2.2,'cubic')ans =4.8133第七章（1）g=inline('sin(x)./x');>> [s,n]=quad(g,0,2)s =1.6054n =18（2）g=inline('1./((x-0.3).*(x-0.3)+0.01)-1./((x-0.9).*(x-0.9)+0.04)-6'); >> [s,n]=quad(g,0,1)s =11.7006n =145。
MATLAB程序设计教程(第二版)课后答案
---MATLAB 第二版课后答案unit3-8unit3实验指导'); 请输入一个三位数：1、n=input('a=fix(n/100); b=fix((n-a*100)/10);c=n-a*100-b*10; d=c*100+b*10+a ）2（1 '); 请输入成绩n=input(' switch ncase num2cell(90:100) p='A';case num2cell(80:89) p='B';case num2cell(70:79)p='C';case num2cell(60:69) p='D'; otherwisep='E'; end price=p请输入成绩）（2 n=input(''); if n>=90&n<=100 p='A';elseif n>=80&n<=89 p='B';elseif n>=70&n<=79p='C';elseif n>=60&n<=69p='D'; elsep='E'; end price=p）3（n;try--------catchprice='erroe' end3n=[1,5,56,4,3,476,45,6,3,76,45,6,4,3,6,4,23,76,908, 6];--------a=n(1)；b=n(1); for m=2:20if n(m)>aa=n(m);elseif n(m)<bb=n(m); end endmax=amin=b 法2n=[1,5,56,4,3,476,45,6,3,76,45,6,4,3,6,4,23,76,908, 6]; min=min(n) max=max(n) 4b=[-3.0:0.1:3.0]; for n=1:61a=b(n);y(n)=(exp(0.3*a)-exp(-0.3*a))/2*sin(a+0.3)+log((0.3+a) /2); end y5 y1=0;y2=1;n=input(' 请输入n 的值：');for i=1:ny1=y1+1/i^2;y2=y2*((4*i*i)/((2*i-1)*(2*i+1))); endy1y2 6A=[1,1,1,1,1,1;2,2,2,2,2,2;3,3,3,3,3,3;4,4,4,4,4,4;5,5 ,5,5,5,5;6,6,6,6,6,6]; n=input(' 请输入n 的值:');if n<=5&n>=0disp(A([n],:));--------elseif n<0 disp(lasterr); else disp(A([6],:));disp(lasterr); end 7（1）f=[];--------for n=1:40f(n)=n+10*log(n^2+5); endy=f(40)/(f(30)+f(20)) （2）f=[];a=0; for n=1:40f(n)=a+n*(n+1); a=f(n);endy=f(40)/(f(30)+f(20)) 8y=0;:'); 输入m 的值m=input(' '); 值：输入n=input('n for i=1:ny=y+i^m; end y********************************************************** ** function s=shi8_1(n,m) s=0;for i=1:ns=s+i^m;end********************************************************** ** shi8_1(100,1)+shi8_1(50,2)+shi8_1(10,1/2) 思考练习2N=[1,2,3,4,5];2.*N N./2 1./N 1./N.^2 3s=fix(100*rand(1,20)*9/10+10)y=sum(s)/20--------j=0; for i=1:20if s(i)<y&rem(s(i),2)==0j=j+1; A(j)=s(i); else continue; end --------endA 4y1=0;y2=0;n=input(' 请输入n 的值：');for i=1:ny1=y1+-(-1)^i/(2*i-1); y2=y2+1/4^i; endy1y2 unit4 实验指导（1 1）x=-10:0.05:10; y=x-x.^3./6;plot(x,y) （2）x=-10:0.5:10;ezplot('x^2+2*y^2-64',[-8,8]);grid on;2t=-pi:pi/10:pi; y=1./(1+exp(-t));subplot(2,2,1);bar(t,y); (t,y)'); title(' 条形图axis([-pi,pi,0,1]);subplot(2,2,2); stairs(t,y,'b'); (t,y)'); 阶梯图title(' axis([-pi,pi,0,1]); subplot(2,2,3); stem(t,y,'k'); (t,y)'); title(' 杆图axis([-pi,pi,0,1]);subplot(2,2,4); loglog(t,y,'y');(t,y)'); title(' 对数坐标图--------3（1）t=0:pi/50:2*pi; r=5.*cos(t)+4;polar(t,r);title('\rho=5*cos\theta+4'); （2）--------t=-pi/3:pi/50:pi/3;r=5.*((sin(t)).^2)./cos(t); polar(t,r); 4（1）t=0:pi/50:2*pi;x=exp(-t./20).*cos(t); y=exp(-t./20).*sin(t); z=t; plot3(x,y,z);grid on; （2）[x,y]=meshgrid(-5:5);z=zeros(11)+5; mesh(x,y,z); shading interp; 5 [x,y,z]=sphere(20);surf(x,y,z); axis off;shading interp; m=moviein(20); for i=1:20 axis([-i,i,-i,i,-i,i]) m(:,i)=getframe;end movie(m,4); 思考练习）（2 1 x=-5:0.1:5; y=(1./(2*pi)).*exp((-(x.^2))/2);plot(x,y); （2）t=-2*pi:0.1:2*pi;x=t.*sin(t);y=t.*cos(t); plot(x,y); grid on;3--------t=0:pi/1000:pi;x=sin(3.*t).*cos(t); y1=sin(3.*t).*sin(t);y2=2.*x-0.5;plot(x,y1,'k',x,y2); hold on;--------k=find(abs(y1-y2)<1e-2);x1=x(k); y3=2.*x1-0.5; plot(x1,y3,'rp'); 4 x=-2:0.01:2;y=sin(1./x); subplot(2,1,1); plot(x,y);subplot(2,1,2);fplot('sin(1./x)',[-2,2],1e-4); 5（1）i=-4*pi:0.1:10; j=12./sqrt(i);polar(i,j);title('{\rho}=12/sqrt(\theta)') （2）a=-pi/6:0.01:pi/6;b=3.*asin(a).*cos(a)./((sin(a)).^3+(cos(a)).^3); polar(a,b); 6（1）[u,v]=meshgrid(-4:0.1:4); x=3.*u.*sin(v); y=2.*u.*cos(v);z=4.*u.^2; subplot(2,1,1); mesh(x,y,z); subplot(2,1,2); surf(x,y,z); （2）[x,y]=meshgrid(-3:0.2:3); z=-5./(1+x.^2+y.^2); subplot(1,2,1); mesh(x,y,z); subplot(1,2,2); surf(x,y,z);unit5 实验指导--------1A=randn(10,5) x=mean(A) y=std(A)Max=max(max(A))Min=min(min(A))--------Sumhang=sum(A,2)SumA=sum(Sumhang) B=sort(A);C=sort(B,2,'descend'); C 2（1）a=0:15:90;b=a./180.*pi; s=sin(b)c=0:15:75;d=c./180.*pi; t=tan(d)');e=input(' 请输入想计算的值：S=sin(e/180*pi) T=tan(e/180*pi)S1=interp1(a,s,e,'spline') T1=interp1(c,t,e,'spline') P1=polyfit(a,s,5); P2=polyfit(c,t,5); S2=polyval(P1,e) T2=polyval(P2,e) ）（2n=[1,9,16,25,36,49,64,81,100];N=sqrt(n);x=input('ji suan zhi : '); interp1(n,N,x,'cubic')3 N=64; T=5;t=linspace(0,T,N); h=exp(-t); dt=t(2)-t(1); f=1/dt; X=fft(t);F=X(1:N/2+1);f=f*(0:N/2)/N;--------plot(f,abs(F),'-*') 4P=[2,-3,0,5,13];Q=[1,5,8]; p=polyder(P) q=polyder(P,Q)[a,b]=polyder(P ,Q)--------5P1=[1,2,4,0,5]; P2=[0,1,2]; P3=[1,2,3];P=P1+conv(P2,P3) X=roots(P)A=[-1,1.2,-1.4;0.75,2,3.5;0,5,2.5]; p=polyval(P ,A)思考练习4A=rand(1,30000); a=mean(A)b=std(A) Max=max(A)Min=min(A) n=0;for i=1:30000if(A(i)>0.5)n=n+1; endendny=n/30000 5p=[45,74,54,55,14;78,98,45,74,12;87,98,85,52,65] [M,S]=max(p) [N,H]=min(p) junzhi=mean(p,1) fangcha=std(p,1,1) zong=sum(p,2);[Max,wei]=max(zong) [Min,wei]=min(zong)[zcj,xsxh]=sort(zong,'descend') 6x=[1:10:101];y=[0,1.0414,1.3222,1.4914,1.6128,1.7076,1.7853,1. 8513,1.9085,1.9590,2.0043];[p,s]=polyfit(x,y,5)--------a=1:5:101; y1=polyval(p,a); plot(x,y,':o',a,y1,'-*') unit6 实验指导1A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6];--------p=[0.95,0.67,0.52]';x=A\pA=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; p=[0.95,0.67,0.53]'; x=A\p cond(A) 2（1）x1=fzero(@funx1,-1)function fx=funx1(x) fx=x^41+x^3+1; （2）x2=fzero(@funx2,0.5)function fx=funx2(x) fx=x-sin(x)/x; （3）options=optimset('Display','off');x=fsolve(@fun3,[1,1,1]',options)q=fun3(x)function q=fun3(p)x=p(1);y=p(2); z=p(3);q(1)=sin(x)+y^2+log(z)-7;q(2)=3*x+2^y-z^3+1;q(3)=x+y+z-5; 3（1）t0=0; tf=5; y0=1;[t,y]=ode23(@fun4,[t0,tf],y0);t'y'function yp=fun4(t,y)yp=-(1.2+sin(10*t))*y; （2）t0=0;--------tf=5; y0=1;[t,y]=ode23(@fun5,[t0,tf],y0);t' y'function yp=fun5(t,y) --------yp=cos(t)-y/(1+t^2);4x=fminbnd(@mymin,0,2);-mymin(x)function fx=mymin(x)fx=-(1+x.^2)/(1+x.^4); 5options=optimset('Display','off');[x,fval]=fmincon(@fun6,[0,0,0],[],[],a,b,lb,ub)-fvalfunction f=fun6(x)f=-(sqrt(x(1))+(400-x(1))*1.1+(sqrt(x(2))+(400-x(1))*1.1 -x(2))*1.1+sqrt(3)+(((400-x(1))*1.1-x(2))*1.1-x(3))*1.1+sqrt(x(x4)));思考练习1（1）A=[2,3,5;3,7,4;1,-7,1]; B=[10,3,5]';C1=inv(A)*B C2=A\B [L,U]=lu(A);x=U\(L\B) （2）A=[6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2];B=[-4,13,1,11]'; C1=inv(A)*BC2=A\B [L,U]=lu(A); x=U\(L\B) 2（1）x1=fzero(@funx1,1.5)function fx=funx1(x) fx=3*x+sin(x)-exp(x); （2）x1=fzero(@funx2,1)function fx=funx2(x) fx=x-1/x+5; （3）--------options=optimset('Display','off');x=fsolve(@fun3,[3,0]',options)q=fun3(x)function q=fun3(p) x=p(1);--------y=p(2);q(1)=x^2+y^2-9;q(2)=x+y-1; 3（1）t0=0;tf=5; y0=[0,1];[t,y]=ode45(@vdpol,[t0,tf],y0);[t,y]function ydot=vdpol(t,y);ydot(1)=(2-3*y(2)-2*t*y(1))./(1+t^2); ydot(2)=y(1); ydot=ydot'; （2）t0=0; tf=5; y0=[1;0;2];[t,y]=ode45(@vdpoll,[t0,tf],y0);[t,y]function ydot=vdpoll(t,y);ydot(1)=cos(t)-y(3)./(3+sin(t))+5*y(1).*cos(2*t)/((t+1).^ 2)-y(2); ydot(2)=y(1);ydot(3)=y(2); ydot=ydot'; 4x=fminbnd(@mymin,0,pi); -mymin(x) function fx=mymin(x) fx=-sin(x)-cos(x.^2); 5[x,y1]=fminbnd(@mymax,0,1.5);-y1function fx=mymax(x);fx=-(9*x+4*x.^3-12*x.^2);unit7 实验指导1（1）--------format longfx=inline('sin(x)./x');[I,n]=quadl(fx,0,2,1e-10) （2）format longfx=inline('1./((x-0.3).^2+0.01)-1./((x-0.9).^2+0.04)-6');--------[I,n]=quad(fx,0,1,1e-10) 2（1）global ki; ki=0;I=dblquad(@fxy,0,1,0,1)ki （2）f=inline('abs(cos(x+y))','x','y'); I=dblquad(f,0,pi,0,pi)3X=0.3:0.2:1.5;F=[0.3895,0.6598,0.9147,1.1611,1.3971,1.6212,1.8 325];trapz(X,F) 4p=0:pi/5:2*pi; for n=1:3nDX=diff(sin(p),n)end5f=inline('sin(x)./(x+cos(2.*x))');g=inline('(cos(x).*(x+cos(2*x))-sin(x).*(1-2.*sin(2*x)))/( x+cos(2.*x)).^2'); x=-pi:0.01:pi; p=polyfit(x,f(x),5);dp=polyder(p); 在假设点的函数% dp 求dpx=polyval(dp,x); 值直接% dx=diff(f([x,3.01]))/ 对f(x) 求数值导数在假设点的导数g 的导函 f % 求函数0.01; gx=g(x); 数作图%plot(x,dpx,x,dx,'.',x,gx,'-'); 思考练习2format longfx=inline('1./(1+x.^2)');[I,n]=quad(fx,-Inf,Inf,1e-10)[I,n]=quadl(fx,-Inf,Inf,1e-10)--------x=-100000:100000;y=1./(1+x.^2); trapz(x,y)format short 3(1)format longfx=inline('log(1+x)./(1+x.^2)');[I,n]=quad(fx,0,1,1e-10)--------(2)format longfx=inline('sqrt(cos(t.^2)+4*sin((2*t).^2)+1)');[I,n]=quad(fx,0,2*pi,1e-10) 4f=inline('4.*x.*z.*exp(-z.^2.*y-x.^2)');I=triplequad(f,0,pi,0,pi,0,1) 5f=inline('sin(x)'); g=inline('cos(x)'); x=0:0.01:2*pi; p=polyfit(x,f(x),5);dp=polyder(p); dpx=polyval(dp,x);dx=diff(f([x,2*pi+0.01]))/0.01; gx=g(x);plot(x,dpx,x,dx,'.',x,gx,'-')unit8 实验指导1 syms x y;s=x^4-y^4; factor(s)factor(5135) 2 syms x; f=(x-2)/(x^2-4);limit(f,x,2)sym x;f=(sqrt(pi)-sqrt(acos(x)))/sqrt(x+1);limit(f,x,-1,'right')3sym x;--------f=sin(1/x); diff(f,'x')diff(f,'x',2)sym x;--------f=(1-cos(2*x))/x;diff(f,'x') diff(f,'x',2) 4 sym x; f=sqrt(exp(x)+1);int(f,'x')syms x y;f=x/(x+y);int(f,'y')sym x;f=exp(x)*(1+exp(x))^2; int(f,'x',0,log(2))sym x; f=x*log(x); int(f,'x',1,exp(1)) 5 sym x; s=symsum((-1)^(x+1)/x,1,Inf)sym y;z=symsum(y^(2*y-1)/(2*y-1),1,Inf)6 sym x;f1=(exp(x)+exp(-x))/2; f2=sqrt(x^3-2*x+1);taylor(f1,x,5,0)taylor(f2,x,6,0) 7syms x y a;x=solve('x^3+a*x+1=0','x')x=solve('sin(x)+2*cos(x)-sqrt(x)=0','x')[x y]=solve('log(x/y)=9','exp(x+y)=3','x','y')--------8 syms n;[x,y]=dsolve('x*(D2y)+(1-n)*(Dy)+y=0','y(0)=0','Dy(0)=0' ,'x') 思考练习2syms x B1 B2 a b--------s1=2*((cos(x))^2)-(sin(x))^2;s2=sin(B1)*cos(B2)-cos(B1)*sin(B2);s3=sqrt((a+sqrt(a^2-b))/2)+sqrt((a-sqrt(a^2-b))/2); s4=(4*x^2+8*x+3)/(2*x+1); h1=simplify(s1)h2=simplify(s2)h3=simplify(s3) h4=simplify(s4) 3syms x a; f=abs(x)/x; limit(f,x,0,'left')f=(x+a/x)^x;limit(f,x,inf) 4syms x y mf=sqrt(x+sqrt(x+sqrt(x))); m=diff(f,'x') diff(m,'x')syms x y;f=x+y-sqrt(x^2+y^2); z1=diff(f,x) z2=diff(f,y)5 syms x y;f=x+y-sqrt(x^2+y^2); z1=diff(f,x)z2=diff(f,y)sym x;f=1/(asin(x)^2*(1-x^2)^(1/2));--------int(f) 6syms xf=1/(1+x); int(f,0,4)--------sym x;f=x^3*sin(x)^2/(x^6+2*x^4+1); int(f,-1,1)sym x;f=x^3*sin(x)^2/(x^6+2*x^4+1); int(f,-1,1)7syms ns=symsum(1/4^n,1,inf)sym n;s=symsum(((n+1)/n)^(1/2),1,inf)eval(y)8syms x f=tan(x);taylor(f,x,3,0)syms x f=sin(x)^2; taylor(f,x,5,0)9syms xx=solve('log(1+x)-5/(1+sin(x))=2','x')syms x y z[x yz]=solve('4*x^2/(4*x^2+1)=y','4*y^2/(4*y^2+1)=z','4*z^2/(4*z^2+1)=x','x','y','z')10[x ,y]=dsolve('Dx=3*x+4*y','Dy=5*x-7*y','x(0)=0','y(0)=1','t')专业资料学习资料教育培训考试建筑装潢资料-----。
数学实验第四章函数和方程习题MATLAB编码(胡良剑版)
%p0,p1是初始值
%delta是给定允许误差
%m是p1-p0的误差估计?
%k是所需
%要的迭代次数?%y=f(p1)
% K=0,p0,p1,feval('f',p0),feval('f',p1)
% for k=1:max1
% if(f1==0)
% root=a;
% end
% if(f2==0)
% root=b;
% end
% if(f1*f2>0)
% disp('两端点函数值乘积大于0!');
% return;
% else
% tol=1;
% fun=diff(sym(f));
% fa=subs(sym(f),findsym(sym(f)),a);
% 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]) %求极大值
%第10题
% t=0:24;
% c=[15 14 14 14 14 15 16 18 20 22 23 25 28 ...
% x(2)=fminsearch(fun2,-2.5);
% x(4)=3;
% feval(fun,x)
%第8题（3）
% fun=inline('abs(x^3-x^2-x-2)');
% fplot(fun,[0 3]);grid on; %作图观察
MATLAB程序设计教程(第二版)课后答案
MATLAB第二版课后答案unit3-8 unit3实验指导1、n=i nput('请输入一个三位数：');a=fix (n/100);b=fix(( n-a*100)/10);c=n-a*100-b*10;d=c*100+b*10+a2 (1)n=input('请输入成绩');switch ncase num2cell(90:100)p='A';case num2cell(80:89)p='B';case num2cell(70:79)p='C';case num2cell(60:69)p='D';otherwisep='E';end price二p(2)n=input('请输入成绩'); if n>=90&n <=100 p='A：elseif n>=80&n v=89p='B';elseif n>=70&n v=79p='C';elseif n>=60&n v=69p='D';elsep='E';endprice=p(3)tryn;catchprice二'erroe'end3n二[1,5,56,4,3,476,45,6,3,76,45,6,4,3,6,4,23,76,908,6]; a=n(1);b=n (1);for m=2:20if n(m )>aa=n( m);elseif n(m)<bb=n( m);endendmax=amin=b法2n二[1,5,56,4,3,476,45,6,3,76,45,6,4,3,6,4,23,76,908,6]; mi n=mi n(n)max二max (n)4b二[-3.0:0.130];for n=1:61a=b( n);y(n )=(exp(0.3*a)-exp(-0.3*a))/2*si n(a+0.3)+log((0.3+a)/2);5y1=0;y2=1;n=input('请输入n的值：');for i=1: ny1=y1+1/i A2;y2二y2*((4*i*i)/((2*i-1)*(2*i+1)));endyiy26A=[1,1,1,1,1,1;2,2,2,2,2,2;3,3,3,3,3,3;4,4,4,4,4,4;5,5,5,5,5,5;6,6,6,6,6,6]; n=input('请输入n的值:');if n<=5&n>=0disp(A([n],:));elseif n<0disp(lasterr);else disp(A([6],:));disp(lasterr);endf=[]；for n=1:40f(n )=n+10*log( n^2+5);endy=f(40)/(f(30)+f(20))(2)f=[];a=0;for n=1:40f(n )=a+n*( n+1);a=f(n);endy=f(40)/(f(30)+f(20))8y=0;m=input('输入m 的值:');n=input('输入n 值：');for i=1: ny=y+i A m;end************************************************************fun cti on s=shi8_1( n,m)s=0;for i=1: ns=s+i A m;end************************************************************ shi8_1(100,1)+shi8_1(50,2)+shi8_1(10,1/2)思考练习2N=[1,2,3,4,5];2.*NN./21./N1./N.A23s=fix(100*ra nd(1,20)*9/10+10)y=sum(s)/20j=0；for i=1:20if s(i)vy&rem(s(i),2)==0j=j+1；A(j)=s(i)；else con ti nue;endendA4y1=0;y2=0;n=input('请输入n的值：'); for i=1: n y1=y1+-(-1)A i/(2*i-1); y2=y2+1/4A i; endyiy2unit4实验指导1 (1)x=-10:0.05:10;y=x-x.A3./6;Plot(x,y)(2)x=-10:0.5:10;ezplot('x A2+2*y A2-64',[-8,8]);grid on;2t=-pi:pi/10:pi;y=1./(1+exp(-t)); subplot(2,2,1);bar(t,y); title('条形图(t,y)');axis([-pi,pi,0,1]); subplot(2,2,2);stairs(t,y,'b');title('阶梯图(t,y)');axis([-pi,pi,0,1]); subplot(2,2,3);stem(t,y,'k');title('杆图(t,y)');axis([-pi,pi,0,1]); subplot(2,2,4);loglog(t,y,'y');title('对数坐标图(t,y)');3 (1)t=0:pi/50:2*pi;r=5.*cos(t)+4;polar(t,r);title('\rho=5*cos\theta+4');(2)t=-pi/3:pi/50:pi/3;r=5.*((si n( t)).A2)./cos(t);polar(t,r);4 (1)t=0:pi/50:2*pi;x=exp(-t./20).*cos(t);y=exp(-t./20).*si n( t);z=t;plot3(x,y,z);grid on;(2)[x,y]=meshgrid(-5:5); z=zeros(11)+5; mesh(x,y,z); shad ing in terp;5[x,y,z]=sphere(20);surf(x,y,z);axis off;shad ing in terp;m=movie in( 20);axis([-i,i,-i,i,-i,i])m(:,i)=getframe;endmovie(m,4);思考练习2 (1)x=-5:0.1:5;y=(1./(2*pi)).*exp((-(x“2))/2); plot(x,y);(2)t=-2*pi:0.1:2*pi;x=t.*si n( t);y=t.*cos(t);grid on;3t=O:pi/1OOO:pi;x=si n( 3.*t).*cos(t);y1=si n(3.*t).*si n( t);y2=2.*x-0.5;plot(x,y1,'k',x,y2);hold on;k二fin d(abs(y1-y2)<1e-2); x1=x(k);y3=2.*x1-0.5;plot(x1,y3,'rp');4x=-2:0.01:2;y=s in (1./x);subplot(2,1,1);plot(x,y);subplot(2,1,2);fplot('si n(1./x)',[-2,2],1e-4);5 (1)i=-4*pi:0.1:10;j=12./sqrt(i);polar(ij);title('{\rho}=12/sqrt(\theta)')(2)a=-pi/6:0.01:pi/6;b=3.*asi n( a).*cos (a) ./((s in (a)).八3+(cos(a)).八3); polar(a,b);6 (1)[u,v]=meshgrid(-4:0.1:4);x=3.*u.*s in (v);y=2.*u.*cos(v);z=4.*u.A2;subplot(2,1,1);mesh(x,y,z);subplot(2,1,2);surf(x,y,z);(2)[x,y]=meshgrid(-3:0.2:3);z=-5./(1+x.A2+y.A2);subplot(1,2,1);mesh(x,y,z);subplot(1,2,2); surf(x,y,z);unit5实验指导1A=ra ndn (10,5)x=mea n(A)y=std(A)Max二max(max(A))Min二min( mi n( A))Sumha ng二sum(A,2) SumA=sum(Sumha ng) B=sort(A); C=sort(B,2,'desce nd'); C2 (1)a=0:15:90;b=a./180.*pi;s=s in(b)c=0:15:75;d二c./180.*pi;t=ta n(d)e=i nput('请输入想计算的值：');S=si n( e/180*pi)T=ta n(e/180*pi)S1=i nterp1(a,s,e,'spli ne')T1二in terp1(c,t,e,'spli ne')P1=polyfit(a,s,5);P2=polyfit(c,t,5);S2=polyval(P1,e)T2=polyval(P2,e)(2)n二[1,9,16,25,36,49,64,81,100]; N=sqrt (n);x=in put('ji sua n zhi :');in terp1( n,N,x,'cubic')3N=64;T=5;t=li nspace(0,T,N);h=exp(-t);dt=t( 2)-t(1);f=1/dt;X=fft(t);F=X(1:N/2+1);f=f*(0:N/2)/N;plot(f,abs(F),'-*')4P=[2,-3,0,5,13];Q=[1,5,8];p=polyder(P)q=polyder(P,Q)[a,b]=polyder(P,Q)5P仁[1,2,4,0,5];P2=[0,1,2];P3=[1,2,3];P=P1+co nv(P2,P3)X=roots(P)A二[-1,1.2,-1.4；0.75,2,3.5;0,5,2.5]; p=polyval(P,A)思考练习A=ra nd(1,30000);a=mea n(A)b二std(A)Max二max(A)Min 二mi n(A)n=0;for i=1:30000if(A(i)>0.5)n=n+1;endendny=n/300005p二[45,74,54,55,14;78,98,45,74,12;87,98,85,52,65] [M,S]=max(p)[N,H]=mi n(p)jun zhi=mea n( p,1)fan gcha=std(p,1,1)zon g=sum(p,2);[Max,wei]=max(z ong)[Mi n, wei]=mi n(zong)[zcj,xsxh]=sort(z on g,'desce nd')6x=[1:10:101];y二[0,1.0414,1.3222,1.4914,1.6128,1.7076,1.7853,1.8513,1.9085,1.9590,2.0043]; [p,s]=polyfit(x,y,5)a=1:5:101;y1=polyval(p,a);plot(x,y,':o',a,y1,'-*')unit6实验指导1A二[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6];p二[0.95,0.67,0.52]';x=A\pA=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6];p二[0.95,0.67,0.53]';x=A\pcon d(A)2 (1)x1二fzero(@fu nx1,-1) fun cti on fx二fun x1(x) fx=x A41+x A3+1;(2)x2=fzero(@fu nx2,0.5)fun cti on fx二fun x2(x)fx=x-si n( x)/x;(3)optio ns=optimset('Display','off); x=fsolve(@fu n3,[1,1,1]',optio ns) q=fu n3(x) fun cti on q二fun 3(p)x=p(1);y=p (2);z=p (3)；q(1)=si n(x)+yA2+log (z)-7;q(2)=3*x+2Ay-zA3+1; q(3)=x+y+z-5;3 (1)t0=0;tf=5;yO=1;[t,y]=ode23(@fu n4 ,[tO,tf],yO); t'y'fun cti on yp二fu n4 (t,y)yp=-(1.2+si n(1O*t))*y;(2)t0=0;tf=5;y0=1;[t,y]=ode23(@fu n5,[t0,tf],y0); t'y'fun cti on yp二fu n5(t,y)yp=cos(t)-y/(1+t A2);4x=fm inbn d(@m ymin, 0,2);-mymi n(x)fun cti on fx=mymi n(x)fx=-(1+x.A2)/(1+x.A4);5optio ns=optimset('Display','off);[x,fval]=fmi neon (@fu n6,[0,0,0],[],[],a,b,lb,ub)-fvalfun eti on f二fun 6(x)f=-(sqrt(x(1))+(400-x(1))*1.1+(sqrt(x (2))+(400-x(1))*1.1-x (2))*1.1+sqrt(3)+(((400- x(1))*1.1-x(2))*1.1-x(3))*1.1+sqrt(x(x4)));思考练习1 (1)A=[2,3,5;3,7,4;1,-7,1];B=[10,3,5]';C1=i nv(A)*BC2=A\B[L,U]=lu(A);x=U\(L\B)(2)A=[6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2]; B二[-4,13,1,11]';C1=i nv(A)*BC2=A\B[L,U]=lu(A);x=U\(L\B)2 (1)x1二fzero(@fu nx1,1.5)fun cti on fx二fun x1(x)fx=3*x+si n(x)-exp(x);(2)x1=fzero(@fu nx2,1)fun cti on fx二fun x2(x) fx=x-1/x+5;(3)optio ns=optimset('Display','off); x=fsolve(@fu n3,[3,0]',optio ns) q=fu n3(x) fun cti on q二fun 3(p)x=p(1);y=p (2)；q(1)=x A2+y A2-9;q(2)=x+y-1;3 (1)t0=0;tf=5;yo=[o,i];[t,y]=ode45(@vdpol,[t0,tf],y0);[t,y]fun cti on ydot二vdpol(t,y);ydot(1)=(2-3* y(2)-2*t*y(1))./(1+t A2);ydot(2)=y(1);ydot=ydot';(2)t0=0;tf=5;y0=[i;0;2];[t,y]=ode45(@vdpoll,[t0,tf],y0);[t,y]fun cti on ydot二vdpoll(t,y);ydot(1)=cos(t)-y(3)./(3+si n( t))+5*y(1).*cos(2*t)/((t+1)A2)-y(2); ydot(2)=y(1);ydot(3)=y(2);ydot=ydot';4x=fmi nbn d(@mymi n,0,pi);-mymi n(x)fun cti on fx=mymi n(x)fx=-si n(x)-cos(x.A2);5[x,y1]=fmi nbn d(@mymax,0,1.5);-yifun cti on fx=mymax(x);fx=-(9*x+4*x.A3-12*x.A2);unit7实验指导1 (1)format longfx二in li ne('s in( x)./x');[I, n ]=quadl(fx,0,2,1e-10)(2)format longfx=i nlin e('1./((x-0.3)A2+0.01)-1./((x-0.9)A2+0.04)-6'); [l,n ]=quad(fx,0,1,1e-10)2 (1)global ki;ki=0;l=dblquad(@fxy,0,1,0,1)ki(2)f=in li ne('abs(cos(x+y))','x','y');匸dblquad(f,O,pi,O,pi)3X=0.3:0.2:1.5;F二[0.3895,0.6598,0.9147,1.1611,1.3971,1.6212,1.8325]; trapz(X,F)4p=0:pi/5:2*pi;for n=1:3DX=diff(s in (p), n)end5f=inlin e('si n(x)./(x+cos(2.*x))');g=i nlin e('(cos(x).*(x+cos(2*x))-si n(x).*(1-2.*si n(2*x)))/(x+cos(2.*x)).八2'); x=-pi:0.01:pi;p=polyfit(x,f(x),5);dp=polyder(p);dpx=polyval(dp,x); %求dp在假设点的函数值dx=diff(f([x,3.01]))/0.01; %直接对f(x)求数值导数gx=g(x); %求函数f的导函数g在假设点的导数plot(x,dpx,x,dx,'.',x,gx,'-'); %作图思考练习2format longfx=i nlin e('1./(1+x.A2)');[I, n]二quad(fx,-l nf,ln f,1e-10)[I, n]=quadl(fx,-I nf,ln f,1e-10) x=-100000:100000;y=1./(1+x.A2);trapz(x,y)format short3(1)format longfx=i nli ne('log(1+x).心+x.A2)');[l,n ]=quad(fx,0,1,1e-10)⑵format longfx=i nli ne('sqrt(cos(t.A2)+4*s in ((2吐).八2)+1)'); [I, n]=quad(fx,0,2*pi,1e-10)4f=in li ne('4.*x.*z.*exp(-z.A2.*y-x.A2)');匸triplequad(f,0,pi,0,pi,0,1)5f=in li ne('s in( x)');g=i nli ne('cos(x)');x=0:0.01:2*pi;p=polyfit(x,f(x),5);dp二polyder(p);dpx=polyval(dp,x);dx二diff(f([x,2*pi+0.01]))/0.01;plot(x,dpx,x,dx,'.',x,gx,'-') unit8实验指导1syms x y;s=x A4-y A4;factor(s)factor(5135)2syms x;f=(x-2)/(xA2-4);Iimit(f,x,2)sym x;f=(sqrt(pi)-sqrt(acos(x)))/sqrt(x+1); limit(f,x,-1,'right')3sym x;f=si n(1/x);diff(f,'x')diff(f,'x',2)f=(1-cos(2*x))/x;diff(f,'x')diff(f,'x',2)4sym x;f=sqrt(exp(x)+1);in t(f,'x')syms x y;f=x/(x+y);in t(f,'y')sym x;f=exp(x)*(1+exp(x))八2;in t(f,'x',0,log(2))sym x;f=x*log(x);in t(f,'x',1,exp(1))5sym x;s=symsum((-1)八(x+1)/x,1,l nf) sym y;z=symsum(y A(2*y-1)/(2*y-1),1,I nf)6sym x;f1=(exp(x)+exp(-x))/2;f2=sqrt(xA3-2*x+1);(f1,x,5,0)(f2,x,6,0)7syms x y a;x=solve('xA3+a*x+1=0','x')x=solve('si n( x)+2*cos(x)-sqrt(x)=0','x')[x y]=solve('log(x/y)=9','exp(x+y)=3','x','y')syms n;[x,y]=dsolve('x*(D2y)+(1-n)*(Dy)+y=O','y(O)=O','Dy(O)=O','x') 思考练习2syms x B1 B2 a bs1=2*((cos(x))A2)-(si n(x))八2;s2=si n(B1)*cos(B2)-cos(B1)*si n(B2);s3二sqrt((a+sqrt(aA2-b))/2)+sqrt((a-sqrt(aA2-b))/2);s4=(4*xA2+8*x+3)/(2*x+1);h1=simplify(s1)h2=simplify(s2)h3=simplify(s3)h4=simplify(s4)3syms x a;f=abs(x)/x;limit(f,x,0,'left')f=(x+a/x)Ax;limit(f,x, inf)4syms x y m f=sqrt(x+sqrt(x+sqrt(x))); m=diff(f,'x') diff(m,'x') syms x y; f=x+y-sqrt(x A2+y A2);z1=diff(f,x)z2=diff(f,y)5syms x y; f=x+y-sqrt(xA2+yA2);z1=diff(f,x)z2=diff(f,y)sym x;f=1/(asi n(xF2*(1-xA2)A(1/2)); int(f)6syms xf=1/(1+x);in t(f,0,4)sym x;f=x A3*si n(xF2/(xA6+2*x%+1);in t(f,-1,1)sym x;f=xA3*si n(xF2/(xA6+2*xA4+1); in t(f,-1,1)7syms ns=symsum(1/4A n,1,i nf)sym n;s=symsum((( n+1)/n)A(1/2),1,i nf) eval(y)8syms xf=ta n( x);(f,x,3,0)syms xf=si n( x)A2;(f,x,5,0)9syms xx=solve('log(1+x)-5/(1+si n(x))=2','x')syms x y z[x yz]=solve('4*x A2/(4*x A2+1)=y','4*y A2/(4*y A2+1)=z','4*z A2/(4*z A2+1)=x','x','y','z')10[x ,y]=dsolve('Dx=3*x+4*y','Dy=5*x-7*y','x(0)=0','y(0)=1','t')。
matlab数学实验胡良剑第四章
第四章%Exercise 1(1)p=[1 1 1];x=roots(p)polyval(p,x) %验算，结果应为零%Exercise 1(2)roots([3 0 -4 0 2 -1])%Exercise 1(3)p=zeros(1,24);p([1 17 18 22])=[5 -6 8 -5]; %这样比直接写24个系数简短。
x=roots(p)polyval(p,x) %验算，结果应为零%Exercise 1(4)p1=[2 3];p2=conv(p1, p1);p3=conv(p1, p2);p3(end)=p3(end)-4; %原p3最后一个分量-4x=roots(p3)polyval(p3,x) %验算，结果应为零%Exercise 2fun=inline('x*log(sqrt(x^2-1)+x)-sqrt(x^2-1)-0.5*x');fzero(fun,2) %注意定义域,初值须大于1%Exercise 3fun=inline('x^4-2^x');fplot(fun,[-2 2]);grid on;fzero(fun,-1),fzero(fun,1),fminbnd(fun,0.5,1.5)%Exercise 4fun=inline('x*sin(1/x)','x');fplot(fun, [-0.1 0.1]);x=zeros(1,10);for i=1:10, x(i)=fzero(fun,(i-0.5)*0.01);end;x=[x,-x]%Exercise 5fun=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])%Exercise 6fun=@(x)[x(1)-0.7*sin(x(1))-0.2*cos(x(2)),x(2)-0.7*cos(x(1))+0.2*sin(x(2))];[a,b,c]=fsolve(fun,[0.5 0.5]) %初值0<x(1)<1, 0<x(2)<1%Exercise 7clear; 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]',[1.5,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]',[1.8,-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]',[3.5,-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])%Exercise 8(1)clear;str='x.^2.*sin(x.^2-x-2)'; %注意数组点运算fun=inline(str);fplot(fun,[-2 2]);grid on; %作图观察x(1)=-2;x(3)=fminbnd(fun,-1,-0.5);x(5)=fminbnd(fun,1,2);fun2=inline(['-',str]);x(2)=fminbnd(fun2,-2,-1);x(4)=fminbnd(fun2,-0.5,0.5);x(6)=2feval(fun,x)%答案: 以上x(1)(3)(5)是局部极小，x(2)(4)(6)是局部极大，从最后一句知道x(1)全局最小，x(2)最大。
【免费下载】MATLAB程序设计教程第二版课后答案
end A 4 y1=0;y2=0; n=input('请输入 n 的值：'); for i=1:n
y1=y1+-(-1)^i/(2*i-1); y2=y2+1/4^i; end y1 y2 unit4 实验指导 1（1） x=-10:0.05:10; y=x-x.^3./6; plot(x,y) （2） x=-10:0.5:10; ezplot('x^2+2*y^2-64',[-8,8]); grid on; 2 t=-pi:pi/10:pi; y=1./(1+exp(-t)); subplot(2,2,1);bar(t,y); title('条形图(t,y)'); axis([-pi,pi,0,1]); subplot(2,2,2); stairs(t,y,'b'); title('阶梯图(t,y)'); axis([-pi,pi,0,1]); subplot(2,2,3); stem(t,y,'k'); title('杆图(t,y)'); axis([-pi,pi,0,1]); subplot(2,2,4); loglog(t,y,'y'); title('对数坐标图(t,y)'); 3（1） t=0:pi/50:2*pi; r=5.*cos(t)+4; polar(t,r); title('\rho=5*cos\theta+4'); （2）
p='A'; case num2cell(80:89)
p='B'; case num2cell(70:79)
p='C'; case num2cell(60:69)
MATLAB程序设计教程（第二版）课后答案
MA TLAB第二版课后答案unit3-8 unit3实验指导1、 n=input('请输入一个三位数：');a=fix(n/100);b=fix((n-a*100)/10);c=n-a*100-b*10;d=c*100+b*10+a2（1）n=input('请输入成绩');switch ncase num2cell(90:100)p='A';case num2cell(80:89)p='B';case num2cell(70:79)p='C';case num2cell(60:69)p='D';otherwisep='E';endprice=p（2）n=input('请输入成绩');if n>=90&n<=100p='A';elseif n>=80&n<=89p='B';elseif n>=70&n<=79p='C';elseif n>=60&n<=69p='D';elsep='E';endprice=p（3）tryn;catchprice='erroe'end3n=[1,5,56,4,3,476,45,6,3,76,45,6,4,3,6,4,23,76,908,6];b=n(1);for m=2:20if n(m)>aa=n(m);elseif n(m)<bb=n(m);endendmax=amin=b法2n=[1,5,56,4,3,476,45,6,3,76,45,6,4,3,6,4,23,76,908,6];min=min(n)max=max(n)4b=[-3.0:0.1:3.0];for n=1:61a=b(n);y(n)=(exp(0.3*a)-exp(-0.3*a))/2*sin(a+0.3)+log((0.3+a)/2);endy5y1=0;y2=1;n=input('请输入n的值：');for i=1:ny1=y1+1/i^2;y2=y2*((4*i*i)/((2*i-1)*(2*i+1)));endy1y26A=[1,1,1,1,1,1;2,2,2,2,2,2;3,3,3,3,3,3;4,4,4,4,4,4;5,5,5,5,5,5;6,6,6,6,6,6]; n=input('请输入n的值:');if n<=5&n>=0disp(A([n],:));elseif n<0disp(lasterr);else disp(A([6],:));disp(lasterr);end7（1）f=[];f(n)=n+10*log(n^2+5);endy=f(40)/(f(30)+f(20))（2）f=[];a=0;for n=1:40f(n)=a+n*(n+1);a=f(n);endy=f(40)/(f(30)+f(20))8y=0;m=input('输入m的值:');n=input('输入n值：');for i=1:ny=y+i^m;endy************************************************************ function s=shi8_1(n,m)s=0;for i=1:ns=s+i^m;end************************************************************ shi8_1(100,1)+shi8_1(50,2)+shi8_1(10,1/2)思考练习2N=[1,2,3,4,5];2.*NN./21./N1./N.^23s=fix(100*rand(1,20)*9/10+10)y=sum(s)/20j=0;for i=1:20if s(i)<y&rem(s(i),2)==0j=j+1;A(j)=s(i);else continue;endendA4y1=0;y2=0;n=input('请输入n的值：'); for i=1:ny1=y1+-(-1)^i/(2*i-1);y2=y2+1/4^i;endy1y2unit4实验指导1（1）x=-10:0.05:10;y=x-x.^3./6;plot(x,y)（2）x=-10:0.5:10;ezplot('x^2+2*y^2-64',[-8,8]); grid on;2t=-pi:pi/10:pi;y=1./(1+exp(-t));subplot(2,2,1);bar(t,y);title('条形图(t,y)');axis([-pi,pi,0,1]);subplot(2,2,2);stairs(t,y,'b');title('阶梯图(t,y)');axis([-pi,pi,0,1]);subplot(2,2,3);stem(t,y,'k');title('杆图(t,y)');axis([-pi,pi,0,1]);subplot(2,2,4);loglog(t,y,'y');title('对数坐标图(t,y)');3（1）t=0:pi/50:2*pi;r=5.*cos(t)+4;polar(t,r);title('\rho=5*cos\theta+4'); （2）t=-pi/3:pi/50:pi/3;r=5.*((sin(t)).^2)./cos(t); polar(t,r);4（1）t=0:pi/50:2*pi;x=exp(-t./20).*cos(t);y=exp(-t./20).*sin(t);z=t;plot3(x,y,z);grid on;（2）[x,y]=meshgrid(-5:5);z=zeros(11)+5;mesh(x,y,z);shading interp;5[x,y,z]=sphere(20);surf(x,y,z);axis off;shading interp;m=moviein(20);for i=1:20axis([-i,i,-i,i,-i,i])m(:,i)=getframe;endmovie(m,4);思考练习2（1）x=-5:0.1:5;y=(1./(2*pi)).*exp((-(x.^2))/2); plot(x,y);（2）t=-2*pi:0.1:2*pi;x=t.*sin(t);y=t.*cos(t);plot(x,y);grid on;3t=0:pi/1000:pi;x=sin(3.*t).*cos(t);y1=sin(3.*t).*sin(t);y2=2.*x-0.5;plot(x,y1,'k',x,y2);hold on;k=find(abs(y1-y2)<1e-2);x1=x(k);y3=2.*x1-0.5;plot(x1,y3,'rp');4x=-2:0.01:2;y=sin(1./x);subplot(2,1,1);plot(x,y);subplot(2,1,2);fplot('sin(1./x)',[-2,2],1e-4);5（1）i=-4*pi:0.1:10;j=12./sqrt(i);polar(i,j);title('{\rho}=12/sqrt(\theta)')（2）a=-pi/6:0.01:pi/6;b=3.*asin(a).*cos(a)./((sin(a)).^3+(cos(a)).^3); polar(a,b);6（1）[u,v]=meshgrid(-4:0.1:4);x=3.*u.*sin(v);y=2.*u.*cos(v);z=4.*u.^2;subplot(2,1,1);mesh(x,y,z);subplot(2,1,2);surf(x,y,z);（2）[x,y]=meshgrid(-3:0.2:3);z=-5./(1+x.^2+y.^2);subplot(1,2,1);mesh(x,y,z);subplot(1,2,2);surf(x,y,z);unit5实验指导1A=randn(10,5)x=mean(A)y=std(A)Max=max(max(A))Min=min(min(A))Sumhang=sum(A,2)SumA=sum(Sumhang)B=sort(A);C=sort(B,2,'descend');C2（1）a=0:15:90;b=a./180.*pi;s=sin(b)c=0:15:75;d=c./180.*pi;t=tan(d)e=input('请输入想计算的值：'); S=sin(e/180*pi)T=tan(e/180*pi)S1=interp1(a,s,e,'spline')T1=interp1(c,t,e,'spline')P1=polyfit(a,s,5);P2=polyfit(c,t,5);S2=polyval(P1,e)T2=polyval(P2,e)（2）n=[1,9,16,25,36,49,64,81,100]; N=sqrt(n);x=input('ji suan zhi : ');interp1(n,N,x,'cubic')3N=64;T=5;t=linspace(0,T,N);h=exp(-t);dt=t(2)-t(1);f=1/dt;X=fft(t);F=X(1:N/2+1);f=f*(0:N/2)/N;plot(f,abs(F),'-*')4P=[2,-3,0,5,13];Q=[1,5,8];p=polyder(P)q=polyder(P,Q)[a,b]=polyder(P,Q)5P1=[1,2,4,0,5];P2=[0,1,2];P3=[1,2,3];P=P1+conv(P2,P3)X=roots(P)A=[-1,1.2,-1.4;0.75,2,3.5;0,5,2.5];p=polyval(P,A)思考练习4A=rand(1,30000);a=mean(A)b=std(A)Max=max(A)Min=min(A)n=0;for i=1:30000if(A(i)>0.5)n=n+1;endendny=n/300005p=[45,74,54,55,14;78,98,45,74,12;87,98,85,52,65][M,S]=max(p)[N,H]=min(p)junzhi=mean(p,1)fangcha=std(p,1,1)zong=sum(p,2);[Max,wei]=max(zong)[Min,wei]=min(zong)[zcj,xsxh]=sort(zong,'descend')6x=[1:10:101];y=[0,1.0414,1.3222,1.4914,1.6128,1.7076,1.7853,1.8513,1.9085,1.9590,2.0043]; [p,s]=polyfit(x,y,5)a=1:5:101;y1=polyval(p,a);plot(x,y,':o',a,y1,'-*')unit6实验指导1A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6];p=[0.95,0.67,0.52]';x=A\pA=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; p=[0.95,0.67,0.53]';x=A\pcond(A)2（1）x1=fzero(@funx1,-1)function fx=funx1(x)fx=x^41+x^3+1;（2）x2=fzero(@funx2,0.5)function fx=funx2(x)fx=x-sin(x)/x;（3）options=optimset('Display','off');x=fsolve(@fun3,[1,1,1]',options)q=fun3(x)function q=fun3(p)x=p(1);y=p(2);z=p(3);q(1)=sin(x)+y^2+log(z)-7;q(2)=3*x+2^y-z^3+1;q(3)=x+y+z-5;3（1）t0=0;tf=5;y0=1;[t,y]=ode23(@fun4,[t0,tf],y0);t'y'function yp=fun4(t,y)yp=-(1.2+sin(10*t))*y;（2）t0=0;tf=5;y0=1;[t,y]=ode23(@fun5,[t0,tf],y0);t'y'function yp=fun5(t,y)yp=cos(t)-y/(1+t^2);4x=fminbnd(@mymin,0,2);-mymin(x)function fx=mymin(x)fx=-(1+x.^2)/(1+x.^4);5options=optimset('Display','off');[x,fval]=fmincon(@fun6,[0,0,0],[],[],a,b,lb,ub)-fvalfunction f=fun6(x)f=-(sqrt(x(1))+(400-x(1))*1.1+(sqrt(x(2))+(400-x(1))*1.1-x(2))*1.1+sqrt(3)+(((400-x(1))*1. 1-x(2))*1.1-x(3))*1.1+sqrt(x(x4)));思考练习1（1）A=[2,3,5;3,7,4;1,-7,1];B=[10,3,5]';C1=inv(A)*BC2=A\B[L,U]=lu(A);x=U\(L\B)（2）A=[6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2];B=[-4,13,1,11]';C1=inv(A)*BC2=A\B[L,U]=lu(A);x=U\(L\B)2（1）x1=fzero(@funx1,1.5)function fx=funx1(x)fx=3*x+sin(x)-exp(x);（2）x1=fzero(@funx2,1)function fx=funx2(x)fx=x-1/x+5;（3）options=optimset('Display','off');x=fsolve(@fun3,[3,0]',options)q=fun3(x)function q=fun3(p)x=p(1);y=p(2);q(1)=x^2+y^2-9;q(2)=x+y-1;3（1）t0=0;tf=5;y0=[0,1];[t,y]=ode45(@vdpol,[t0,tf],y0);[t,y]function ydot=vdpol(t,y);ydot(1)=(2-3*y(2)-2*t*y(1))./(1+t^2);ydot(2)=y(1);ydot=ydot';（2）t0=0;tf=5;y0=[1;0;2];[t,y]=ode45(@vdpoll,[t0,tf],y0);[t,y]function ydot=vdpoll(t,y);ydot(1)=cos(t)-y(3)./(3+sin(t))+5*y(1).*cos(2*t)/((t+1).^2)-y(2); ydot(2)=y(1);ydot(3)=y(2);ydot=ydot';4x=fminbnd(@mymin,0,pi);-mymin(x)function fx=mymin(x)fx=-sin(x)-cos(x.^2);5[x,y1]=fminbnd(@mymax,0,1.5);-y1function fx=mymax(x);fx=-(9*x+4*x.^3-12*x.^2);unit7实验指导1（1）format longfx=inline('sin(x)./x');[I,n]=quadl(fx,0,2,1e-10)（2）format longfx=inline('1./((x-0.3).^2+0.01)-1./((x-0.9).^2+0.04)-6');[I,n]=quad(fx,0,1,1e-10)2（1）global ki;ki=0;I=dblquad(@fxy,0,1,0,1)ki（2）f=inline('abs(cos(x+y))','x','y');I=dblquad(f,0,pi,0,pi)3X=0.3:0.2:1.5;F=[0.3895,0.6598,0.9147,1.1611,1.3971,1.6212,1.8325];trapz(X,F)4p=0:pi/5:2*pi;for n=1:3nDX=diff(sin(p),n)end5f=inline('sin(x)./(x+cos(2.*x))');g=inline('(cos(x).*(x+cos(2*x))-sin(x).*(1-2.*sin(2*x)))/(x+cos(2.*x)).^2');x=-pi:0.01:pi;p=polyfit(x,f(x),5);dp=polyder(p);dpx=polyval(dp,x); %求dp在假设点的函数值dx=diff(f([x,3.01]))/0.01; %直接对f(x)求数值导数gx=g(x); %求函数f的导函数g在假设点的导数plot(x,dpx,x,dx,'.',x,gx,'-'); %作图思考练习2format longfx=inline('1./(1+x.^2)');[I,n]=quad(fx,-Inf,Inf,1e-10)[I,n]=quadl(fx,-Inf,Inf,1e-10)x=-100000:100000;y=1./(1+x.^2);trapz(x,y)format short3(1)format longfx=inline('log(1+x)./(1+x.^2)');[I,n]=quad(fx,0,1,1e-10)(2)format longfx=inline('sqrt(cos(t.^2)+4*sin((2*t).^2)+1)'); [I,n]=quad(fx,0,2*pi,1e-10)4f=inline('4.*x.*z.*exp(-z.^2.*y-x.^2)');I=triplequad(f,0,pi,0,pi,0,1)5f=inline('sin(x)');g=inline('cos(x)');x=0:0.01:2*pi;p=polyfit(x,f(x),5);dp=polyder(p);dpx=polyval(dp,x);dx=diff(f([x,2*pi+0.01]))/0.01;gx=g(x);plot(x,dpx,x,dx,'.',x,gx,'-')unit8实验指导1syms x y;s=x^4-y^4;factor(s)factor(5135)2syms x;f=(x-2)/(x^2-4);limit(f,x,2)sym x;f=(sqrt(pi)-sqrt(acos(x)))/sqrt(x+1);limit(f,x,-1,'right')3sym x;f=sin(1/x);diff(f,'x')diff(f,'x',2)sym x;f=(1-cos(2*x))/x;diff(f,'x')diff(f,'x',2)4sym x;f=sqrt(exp(x)+1);int(f,'x')syms x y;f=x/(x+y);int(f,'y')sym x;f=exp(x)*(1+exp(x))^2;int(f,'x',0,log(2))sym x;f=x*log(x);int(f,'x',1,exp(1))5sym x;s=symsum((-1)^(x+1)/x,1,Inf)sym y;z=symsum(y^(2*y-1)/(2*y-1),1,Inf)6sym x;f1=(exp(x)+exp(-x))/2;f2=sqrt(x^3-2*x+1);taylor(f1,x,5,0)taylor(f2,x,6,0)7syms x y a;x=solve('x^3+a*x+1=0','x')x=solve('sin(x)+2*cos(x)-sqrt(x)=0','x')[x y]=solve('log(x/y)=9','exp(x+y)=3','x','y')8syms n;[x,y]=dsolve('x*(D2y)+(1-n)*(Dy)+y=0','y(0)=0','Dy(0)=0','x') 思考练习2syms x B1 B2 a bs1=2*((cos(x))^2)-(sin(x))^2;s2=sin(B1)*cos(B2)-cos(B1)*sin(B2);s3=sqrt((a+sqrt(a^2-b))/2)+sqrt((a-sqrt(a^2-b))/2); s4=(4*x^2+8*x+3)/(2*x+1);h1=simplify(s1)h2=simplify(s2)h3=simplify(s3)h4=simplify(s4)3syms x a;f=abs(x)/x;limit(f,x,0,'left')f=(x+a/x)^x;limit(f,x,inf)4syms x y mf=sqrt(x+sqrt(x+sqrt(x)));m=diff(f,'x')diff(m,'x')syms x y;f=x+y-sqrt(x^2+y^2);z1=diff(f,x)z2=diff(f,y)5syms x y;f=x+y-sqrt(x^2+y^2);z1=diff(f,x)z2=diff(f,y)sym x;f=1/(asin(x)^2*(1-x^2)^(1/2));int(f)6syms xf=1/(1+x);int(f,0,4)sym x;f=x^3*sin(x)^2/(x^6+2*x^4+1);int(f,-1,1)sym x;f=x^3*sin(x)^2/(x^6+2*x^4+1);int(f,-1,1)7syms ns=symsum(1/4^n,1,inf)sym n;s=symsum(((n+1)/n)^(1/2),1,inf)eval(y)8syms xf=tan(x);taylor(f,x,3,0)syms xf=sin(x)^2;taylor(f,x,5,0)9syms xx=solve('log(1+x)-5/(1+sin(x))=2','x')syms x y z[x y z]=solve('4*x^2/(4*x^2+1)=y','4*y^2/(4*y^2+1)=z','4*z^2/(4*z^2+1)=x','x','y','z') 10[x ,y]=dsolve('Dx=3*x+4*y','Dy=5*x-7*y','x(0)=0','y(0)=1','t')。