资源描述:
《matlab习题——2》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、第二章1。%3X3矩阵 a1=eye(3) a2=ones(3) a3=zeros(3) a4=1-2*rand(3) a5=2*randn(3)+1 %15X8矩阵 b1=eye(15,8) b2=ones(15,8) b3=zeros(15,8) b4=1-2*rand(15,8) b5=2*randn(15,8)+12。%生成a矩阵 x=diag([2,-7,8]); xx=diag([3,5],-1); a=rot90(x+xx) %生成b矩阵 x=diag([258]); x
2、x=rot90(diag([704])); b=x+xx %变成行向量 ah=reshape(a,1,9) bh=reshape(b,1,9)3。a=5-10*rand(50,2); b=round(10*a)/10; c=num2str(b,'%+10.1f')4。alpha=-pi:pi/180:pi; x=sin(alpha); y=cos(alpha); %画曲线 plot(alpha,x,'b-') holdon lot(alpha,y,'r-')5。a=rand(10); b=a>
3、=0.6; total=sum(sum(b))6。a=randn(10); b=and(a<.5,a>-.5); total=sum(sum(b))7。%输入a和b a=input('pleaseinputa:'); b=input('pleaseinputb:'); %分支结构 ifand(a<1,b<=.5) disp('语句1') elseifand(a<1,b>.5) disp('语句2') elseifand(a>1,b<=.5) disp('语句3') else disp('语句
4、4') end8。a=rand(10)+0.2; i=find(a>1); b=a(i)8b。a=rand(10)+0.2; x=a>1; =a.*x; [i,j,b]=find(c); b9。%产生测量阵 a=30*rand(100,3)+5; a(97)=103; a(156)=280; %查找删除 [i,j]=find(a>100); a(i,:)=[]10。%生成10X10矩阵,包含全零行,全零列 a=round(10*rand(10)); a(4,:)=0; a(:,7)
5、=0; %查找全零行,全零列 hang=~any(a'); lie=~any(a); hang_hao=find(hang); lie_hao=find(lie); %删除! a(hang_hao,:)=[] a(:,lie_hao)=[]第三章1.%对数曲线 x=0.1:0.01:100; y=log2(x); semilogx(x,y) title('对数曲线') xlabel('log(x)') ylabel('y') gtext('leftarrow对数曲线') %
6、指数曲线 x=0:0.1:10; y=exp(x); figure(2) plot(x,y) title('指数曲线') xlabel('x') ylabel('y') gtext('指数曲线rightarrow')2.x=0:pi/50:4*pi; k=1:25:201; x(k)=[]; subplot(1,2,1) plot(x,sec(x),'r-') gridon title('正割曲线') subplot(1,2,2) plot(x,csc(x),'g-') gridon
7、 title('余割曲线')3.%对数坐标 x=0.01:0.01:100; y=log10(x); figure(1) subplot(2,1,1) plot(x,y,'k-') gridon title('迪卡尔坐标系中y=log_{10}(x)曲线'),ylabel('y') subplot(2,1,2) semilogx(x,y,'k-'),gridon title('半对数坐标系中y=log_{10}(x)曲线') xlabel('x'),ylabel('y') %极坐标 a=ra
8、ndn(2,2);b=eig(a) c1=abs(b),c2=angle(b) figure(1) subplot(2,1,1) plot(b,'rx'),gridon title('迪卡尔坐标系中特征值的表示') subplot(2,1,2) polar(c2,c1,'rx') gtext('极坐标中的特征值的表