matlab第一章作业

发布 2022-06-30 06:35:28 阅读 4058

第一章作业。

第5题。在主窗口输入:

a=exp(1)

a=piexp(pi)-pi^exp(1) %比较eπ和πe 的大小。

第6题。a=rand(3,3) %产生随机数矩阵。

b=inv(a) %求a 的逆矩阵。

c=b*ad=b.*a

e=sqrt(c)

f=abs(e)

format short %设置数据的形式。

第7题。x=[1:1:100];

y=x.^2+sin(x)+log(x)

第8题。x=0:0.1:6;

a=[x;4*x];

b=sin(a);

plot(x,b(1,:)r',x,b(2,:)b')

第9题。3+(5*i)/(1+i));

sin(4+5*i)

log(3+2*i)

log10(3+2*i)

sinh(3+2*i)

第10题。a=

b=zeros(3,3)

c=cell(3)

第11题。a=magic(3)

a=triu(a,1) %取矩阵a的上三角部分。

b= tril(a, -1) %取矩阵a的下三角部分。

c= sum(a)

d=sum(a,2)

e=median(a)

f=diag(diag(a))

g=inv(a)

h=fliplr(a)

i=flipud(a)

j=rot90(a)

k= sort(a,1)

对各元素进行计算的有a、b、e; 对整列元素进行计算的有c;

对整个矩阵进行计算的有f、g、h、i、j

第12题。a=[1,2,3;4,5,6;7,8,9];

b=[-2;-2;-2];

x=a\b第13题。

a=[1,2,3];

b=diag(a)

c=diag(a,1)

d=diag(a,-1)

第14题。频率比为有理数。

> t=0:0.01:5;

> x=sin(1.*t);y=sin(2.*t);

> comet(x,y)

频率为无理数。

> x=sin(1.*t);y=sin(pi.*t);

> comet(x,y)

也可以写成程序,编程如下:

clear;

t=0:0.005:1000;

w1=2;w2=[1,2,4,3,5,7,sqrt(3),sqrt(5),exp(1)];

x=sin(w1*t);

for n=1:9;

y=sin(w2(n)*t);

subplot(3,3,n)

plot(x,y)

end第15题。

clear all

x=[-1 0 1 -1];

y=[0 1 0 0];

z=[-1 i 1 -1];

subplot(1,2,1)

plot(x,y,'r-+'linewidth',3)

subplot(1,2,2)

plot(z ,'g:o','linewidth',5)

第16题。t=0:0.1:10;

x=8*cos(t);

y=4*sqrt(2)*sin(t);

z=-4*sqrt(2)*sin(t);

comet3(x,y,z)

第17题。fplot('cos(tan(pi*x))'0 2*pi])

第18题。ezplot('exp(y)+cos(x)./x+y=0')

第19题。s=[0:8];t=[0:5*pi];

ezsurf('exp(-s)*cos(t)',exp(-s)*sin(t)',t')

第20题。subplot(1,2,1);

x,y,z]=ellipsoid(0,0,1,2,3,4,50);

surf(x,y,z)

subplot(1,2,2)

t=linspace(0,pi,100);

x,y,z]=cylinder(1+t.^2);

surf(x,y,z)

第21题。x,y]=meshgrid(-5:0.5:5,0:0.2:10);

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

surf(x,y,z)

第22题。subplot(1,2,1)

z=cplxgrid(20);

cplxmap(z,exp(z));

colorbar('vert')

subplot(1,2,2)

z=cplxgrid(20);

cplxroot(5)

colorbar('vert')

第23题。clear;clf;

n=1:4;

x=0:0.01:1;

n x]=meshgrid(n,x)

y1=sin(pi*n.*x);

y2=cos(pi*n.*x);

y3=sin(pi*(n+1/2).*x);

y4=cos(pi*(n+1/2).*x);

y=cat(3,y1,y2,y3,y4每个矩阵单独串接成高维数组。

for k=1:4

subplot(2,2,k)

plot(x,y(:,k),'b')

axis([0,1,-1,1])

endclear;clf;

n=1:4;

x=0:0.01:1;

t=0:0.02*pi:0.8*pi;

n x t]=meshgrid(n,x,t);

for k=1:40

subplot(2,2,1)

y1=sin(pi*n(:,1).*x(:,1)).cos(pi*n(:,1).*t(:,k));

plot(x,y1);

axis([0,1,-1,1])

m1(:,k)=getframe;

subplot(2,2,2)

y2=cos(pi*n(:,1).*x(:,1)).cos(pi*n(:,1).*t(:,k));

plot(x,y2);

axis([0,1,-1,1])

m2(:,k)=getframe;

subplot(2,2,3)

y3=sin(pi*(n(:,1)+1/2).*x(:,1)).cos(pi*n(:,1).*t(:,k));

plot(x,y3);

axis([0,1,-1,1])

m3(:,k)=getframe;

subplot(2,2,4)

y4=cos(pi*(n(:,1)+1/2).*x(:,1)).cos(pi*n(:,1).*t(:,k));

plot(x,y4);

axis([0,1,-1,1])

m4(:,k)=getframe;

endmovie(m1,1);

movie(m2,1);

movie(m3,1);

movie(m4,1);

第24题。t=0:pi/20:16*pi;

x=-4*pi:0.1*pi:4*pi;

m=moviein(160); 在每一步走pi/20的基础上,走完整个全程8pi需要160步

for i=1:160

%保证波是传播的形式过来,没到处振幅为零。

%向右传播的波在x>(-4pi+t)时振幅为零,向左传播的波在x<(4pi-t)时振幅为零。

y1=sin(x+t(i)).x<-4*pi+t(i));运用逻辑判断,括号内条件成立时得到1,不成立得到0.

y2=sin(x-t(i)).x>4*pi-t(i));

y=y1+y2;

y1=sin(x+t(i)) 这样只是波形的传播,而不是波的传播。

y2=sin(x-t(i));

plot(x,y1,x,y2,x,y);

axis([-12,12,-2,2]);

m(:,i)=getframe;

endmovie(m,1,20)

第25题。程序1:两列三角波相向而行,相遇叠加。

x=-5*pi:0.1*pi:5*pi;

t=0:pi/10:10*pi;

m=moviein(100); 在每一步走pi/10的基础上,走完整个全程10pi需要100步

for i=1:100

%保证波是传播的形式过来,没到处振幅为零。

%向右传播的波在x>(-5pi+t)时振幅为零,向左传播的波在x<(5pi-t)时振幅为零。

%运用逻辑判断,括号内条件成立时得到1,不成立得到0.

y1=sawtooth(x+t(i),0.5).*x<-5*pi+t(i));

y2=sawtooth(x-t(i),0.5).*x>5*pi-t(i));

y=y1+y2;

plot(x,y);

axis([-15,15,-2,2]);

m(:,i)=getframe;

endmovie(m,1,20)

程序2:两个三角波相向而行,相遇叠加。

clear all

clcx=-10:.025:10;

for t=-7:0.01:7

y1=(x+1+t).*1 y2=(x+1-t).*1 y=y1+y2;

plot(x,y);

hold on

plot([-7 7],[0 0],'r.')用红色标注a,b两点。

hold off

text([-7 7],[0 0],,horizontalalignment','center','verticalalignment','top','fontunits','normalized','fontsize',.1);

axis([-10,10,-1,2]);

drawnow

end第26题。

程序1:% 牛顿环干涉图样。

clear

r=1.07;n=400;lamda=6328e-6;rr=0.1;

x,y]=meshgrid(linspace(-rr,rr,n));

r=abs(x+i*y);

d=r.^2/r/lamda*pi*2; %两相干光的相位差。

z=cos(d);z=abs(z);

z(:,1)=z/sqrt(2);z(:,2)=z/sqrt(2);

z(:,3)=zeros(n);close all;

h=imshow(z);

t=0;k=1;

set(gcf,'doublebuffer','on');设置当前图形窗口,启动双缓存,一般是为动画演示做准备

title('newton---loop');

matlab第一章作业

1 用maltab命令和m程序两种方式求解三元一次方程组。例题 2a b 3c 5 3a b 5c 5 4a b c 92 用maltab命令和m程序两种方式求解3 3复数矩阵的模和相角。创建一个4 4矩阵,求逆矩阵 产生对角矩阵 求矩阵的秩等。4 创建两个4 4矩阵进行加 减 乘法和左除运算。5 ...

MATLAB第一章作业答案

第一章 matlab概况与基本操作。1.选择题 每题2分,共20分 1 最初的matlab核心程序是采用d语言编写的。2 即将于2011年9月发布的matlab新版本的编号为c。2011ra 2011rb r2011b r2011a 3 在默认设置中,matlab中的注释语句显示的颜色是b。a.黑色...

第一章作业

1 设计操作系统的主要目标是什么?1 方便性。使计算机更易于用户使用。2 有效性。以有效的方式管理计算机系统的资源,合理地组织计算机的工作流程,以防止。对计算机资源的不当或错误使用。这是操作系统可用的关键因素。3 可扩展性。为用户的开发搭建一个平台,允许修改并引进新的功能。2 操作系统的基本功能是什...