matlab作业

2. tic是计时器开始语句，toc是计时器停止语句。该程序段用于统计中间三条语句耗费的运行时间。

a=rand(500)生成一个500*500的double格式矩阵，b对a矩阵取逆，ans存储的是一个范数，这个范数对应的矩阵是a左乘b得到的矩阵与单位阵相减的结果。计算出来的用时每次都不一样（有时是1.192060秒，有时是0.

237335秒，或者是别的）。

3. syms x s; f=x^5+3*x^4+4*x^3+2*x^2+3*x+6,f=subs(f,x,(s-1)/(s+1))

得出结果是，原来f=x^5 + 3*x^4 + 4*x^3 + 2*x^2 + 3*x + 6，替换结果是f= (3*(s - 1))/s + 1) +2*(s - 1)^2)/(s + 1)^2 + 4*(s - 1)^3)/(s + 1)^3 + 3*(s - 1)^4)/(s + 1)^4 + s - 1)^5/(s + 1)^5 + 6。

4. a=[1:4;4:

-1:1;2:4,1;3 2 4 1]; b=[1+4*j 2+3*j 3+2*j 4+j;4+j 3+2*j 2+3*j 1+4*j;2+3*j 3+2*j 4+j 1+4*j;3+2*j 2+3*j 4+j 1+4*j];

a(5; 6) =5会出错，提示括弧不匹配。a(5,6)=5运行后矩阵变为5行6列，原矩阵部分不动，a(5,6)值为5，其余新扩充出来的部分都是0。

5. a=magic(8), b=a(2:2:length(a),:结果是。a =

b =6. if x>d

y=h;elseif abs(x)<=d

y=h/(d*x);

elsey=-h

end7. 数值方法应为s=0.5*(2^64-1)。使用符号运算求取精确值方法：

syms q

myans=1/2*(q^64-1); myans=subs(myans,q,2);

结果是9.2234e+018。

8. function v=mat_add(varargin)

y=0;for s=1:length(varargin)

y=y+varargin;

endv=y;

end9. function y=fib(k)

if k==1|k==2

y=1;elseif ~isinteger(int8(k))|k<1

disp('输入应该是正整数');

elsey=fib(k-1)+fib(k-2);end

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

for p=2:30001

y(p)=0.3*x(p-1);

x(p)=1+y(p-1)-1.4*x(p-1)*x(p-1);

endplot(x,y,’.xlabel('x'),ylabel('y'),title('henon引力线图')

t=-1:0.01:1;y=sin(1./t);

plot(t,y),title('y=sin(1/t)的等距绘制，横轴间隔0.01'),xlabel('t'),ylabel('y')

t2=[-1:0.01:-0.2,-0.2:0.001:0.2,0.2:0.01:1];y2=sin(1./t2);

plot(t2,y2),title('y=sin(1/t)的等距绘制，t为-0.2~0.2时间距0.

001，其他时候横轴间隔0.01'),xlabel('t'),ylabel('y')

theta=0:0.01:2*pi;

rho=1.0013.*theta.*theta;

subplot(2,2,1),polar(theta,rho),title('fig 1'),xlabel('\theta'),ylabel('ho');

rho=cos(3.5*theta);

subplot(2,2,2),polar(theta,rho),title('fig 2'),xlabel('\theta'),ylabel('ho');

rho=sin(theta)./theta;

subplot(2,2,3),polar(theta,rho),title('fig 3'),xlabel('\theta'),ylabel('ho');

rho=1-cos(7*theta).*cos(7*theta).*cos(7*theta);

subplot(2,2,4),polar(theta,rho),title('fig 4'),xlabel('\theta'),ylabel('ho');

x,y]=meshgrid(-3:0.01:3);

z=x.*y;zz=sin(z);

subplot(1,2,1),surf(x,y,z),title('z=xy的surf图'),xlabel('x'),ylabel('y'),zlabel('z');

subplot(1,2,2),contour3(x,y,z,30),title('z=xy的等值线图（30条）')xlabel('x'),ylabel('y'),zlabel('z');

figure

subplot(1,2,1),surf(x,y,zz),title('z=sin(xy)的图'),xlabel('x'),ylabel('y'),zlabel('z')

subplot(1,2,2),contour3(x,y,zz,30),title('z=sin(xy)的等值线图（30条）')xlabel('x'),ylabel('y'),zlabel('z');

matlab作业

matlab作业

matlab作业

matlab作业

其他用户还读了