MATLAB期末.docx

1、 计算：w=；（2）z=，其中。w=sqrt(2)*(cos(61/180*pi)+0.34*10(-6)(1/3) t=1 2+3iz=0.5*exp(2*t).*log(t+sqrt(1+t.2)2.若x= 2 8 5; 9 7 1, b=2 4 5:a) 利用max和min函数求x的最大值和最小值.b) 求x每行的中间值c) 将b添加为x的第3行x=2 8 5;9 7 1;b=2 4 5;max(x) min(x) ans =9 8 5 ans =2 7 1 max(max(x) min(min(x) ans = 9 ans=1median(x,2) ans = 5;7 x(3,:)=b ans: x = 2 8 5; 9 7 1; 2 4 5注意：median(M),每一列返回一个值,为M该列的从大到小排列的中间值.median(M,dim),dim为1，2。其中1表示按每列返回一个值,为该列从大到小排列的中间值,2表示按每行返回一个值,为该行从大到小排列的中注意:如果行或列的个数为偶数,返回中间两个值的平均值。6、 已知A=2 , 2.1, 2.5 , 2.7 ;20, -20.4, -20.6 , -20.5 ; 执行ceil(A) , round(A) , floor(A), fix(A), 分析他们的异同，及各自功能。ceil为向上取整，即所得矩阵中每一元素对应为矩阵A中每一元素的大于其本身的最小整数。Round为取不大于其本身的最大整数。floor为向下取整。Fix为去掉小数点。生成一个以向量为基础向量的范得蒙矩阵. vander(1;2;3;5)ans= 1 1 1 1 8 4 2 1 27 9 3 1 125 25 5 1num = 28; s = There are int2str(num) characters here s = sprintf(There are %d characters heren, num)s =There are 28characters heres =There are 28 characters here由输出结果我们可以看出两者的区别是28与characters之间是否有空格。其中。代码：a=i+1; b=-7.9; y=2*pi*a2*(1-pi/4)*b-(0.8333-exp(2)/2)*a)即a ,b2 ,abc ,a2b , a_b_1 为有效数字；1a , _abc , a-b为无效数字。name1 = char(Liu Ying, Hu Xu)name1 =Liu YingHu Xu char表示转换为字符串，可以直接运行，会自动分行name11 = Liu Ying; Hu Xu 。name11 =Liu YingHu Xu功能与char一样，可以看到在Hu Xu后有空格是为保持与Liu Ying字符长度一致，否则不能运行。name2 = Liu Ying, Hu Xuname2 =Liu YingHu Xu可以看到Liu YingHu Xu是直接串起来的，没有分行，就像矩阵一样，可以直接运行。name3 = Liu Ying; Hu Xu? Error using = vertcatCAT arguments dimensions are not consistent.无法运行，因为前后字符长度不一样，若将改为则可以运行如下所示。name3 = Liu Ying; Hu Xuname3 = Liu Ying Hu Xustr=answer:; x=1;y=2.5; ss=sprintf(%sn%d +%.1f=%.2fn,str,x,y,x+y) fprintf(%sn%d +%.1f=%.2fn,str,x,y,x+y)ss =answer:1 +2.5=3.50answer:1 +2.5=3.50两者虽然表示方法不一样，其实功能是没什么区别的。6、 已知str= S0619-S1914-S0388-S0348-S0392-S0429-S0436,试使用函数findstr或find及str2num将字符串str中的7个数字619、1914 自动存储到向量V中。str= S0619-S1914-S0388-S0348-S0392-S0429-S0436;id=find(str=S); N=numel(id);V=zeros(N,1);for i=1:N V(i)=str2double(str(id(i)+1:id(i)+4);enddisp(V) 619 1914 388 348 392 429 4367、DNA序列可以看成是由a、c、g、t四个字符按一定顺序排成的字符串，四个字符在DNA中出现的频率可以作为DNA的一个特征。统计如下DNA序列中a、c、g、t各自出现的频率。aggcacggaaaaacgggaataacggaggaggacttggcacggcattacacggaggacgagdna=aggcacggaaaaacgggaataacggaggaggacttggcacggcattacacggaggacgag;N=length(dna);pa=sum(dna=a)/Npc=sum(dna=c)/Npg=sum(dna=g)/Npt=sum(dna=t)/Npa =0.3500pc =0.1833pg =0.3833pt =0.0833已知某3位学生2门课程的成绩如下表。请创建一结构体变量stu，使其4个域名分别为Number、Name、 English、Computer，用于存储表格中数据。NumberNameEnglishComputer1001Zhang san71681002Li si53711003Wang wu6980stu.Number=1001,1002,1003stu.Name=Zhang san,Li siWang wustu.English=71,53,69stu.Computer=68,71,801、 使用subplot在同一图形窗口绘制如下曲线：1）y=x-x3/3, -6x6 2)y= ,-6x6， 3) x2+2y2=64, -8x8 4), -2px2p 。figuresubplot(221)ezplot(x-x3/3,-6 6)subplot(222)ezplot(exp(-x2/2)/(2*pi),-6 6)subplot(223)ezplot(x2+2*y2-64,-8 8)subplot(224)ezplot(t*sin(t),t*cos(t),-2*pi 2*pi)2、 在同一坐标轴内使用plot函数绘制曲线1） y=2x-0.5，-1x1； 2)，并使用ledgend进行标注。figurex1=-1:0.1:1;y1=2*x1-0.5;t=0:0.1:pi;x2=cos(t).*sin(3*t);y2=sin(t).*sin(3*t);plot(x1,y1,x2,y2)legend(y=2x-0.5,x=costsin3t,y=sintsin3t)3 p/34、绘制三维曲线：t=0:0.1:6*pi;x=exp(-t/20).*cos(t); y=exp(-t/20).*sin(t);z=t; plot3(x,y,z) grid on5、分别使用mesh和ezmesh(（或surf和ezsurf）绘制曲面：f(x,y)=-5/(1+x2+y2), |x|3,|y|3.figureezsurf(5/(1+x2+y2),-3 3 -3 3)x=-3:0.1:3;y=-3:0.1:3;x y=meshgrid(x,y);z=5./(1+x.2+y.2);mesh(x,y,z)四、 将窗口按两行两列方式分割成四个子窗口，在各窗口依次按要求画出如下图型：1）用plot在上以绿色的虚线绘制；2）用fplot在上将曲线绘在同一窗口；3）用ezplot画出曲线，；4）用ezpolar绘制曲线，。（20分）代码：figuresubplot(221)x1=-pi:0.1:pi;y1=cos(2*x1)+sin(3*x1);plot(x1,y1,g-)subplot(222)function Y = myfun(x)Y(:,2) = atan(x(:);Y(:,3)=acos(2*cos(x(:)fh = myfun;fplot(fh,-10 10)subplot(223)ezplot(2*t1.2,3*t1.3,0 7)subplot(224)ezpolar(cos(2*t2).*log(1+t2),0 6*pi) 先执行rand(seed,100);A=rand(100)0.5;将矩阵A按“田”字形均分为4块，统计各块中的非零元的个数。（15分）代码：rand(seed,100);A=rand(100)0.5;A1=A(1:50,1:50);A2=A(51:100,1:50);A3=A(1:50,51:100);A4=A(51:100,51:100);a1=find(A1,0);aa1=numel(a1)a2=find(A2,0);aa2=numel(a2)a3=find(A3,0);aa3=numel(a3)a4=find(A4,0);aa4=numel(a4)用if 语句实现下列分段函数，并计算其在x= -1、1、10时的值：function y=func1(x)if x=0 & x10 & x=2 & x=3 y=x*x-5*x+6;else y=x2-x-1;end编好程序后需保存起来，且运行时的当前文件夹需和程序所在位置一致，否则不能运行。运行结果如下： func1(-1)ans = -6 func1(1)ans = 2 func1(10)ans = 892、 利用for循环求1！+2！+3！+ +5！的值。function s=func2(n)s=0;jc=1;for t=1:n jc=jc*t; s=s+jc;end三、编程求满足的最小m值。function m=func3(T)s=0;i=1;while sv=myvander(2 3 4 5)得v=function V=myvander(x)n=length(x);V=ones(n);for r=2:n V(r,:)=V(r-1,:).*x;end1、编写程序实现函数 m,r,c=mymax(A)，求给定矩阵A的最大值m及最大值所在行列坐标r、c。function m,r,c=mymax(A)row,col=size(A)m=A(1,1);r=1;c=1;for i=1:row for j=1:col if A(i,j)m; m=A(i,j); r=i; c=j; end endend运行结果： A=12 13 23 5 6;7 5 9 16 7;28 4 9 6 5;mymax(A)row =3col =5ans =281 利用Lingo求解：minZ=x1+0.1*x2+0.2*x3+0.3*x4+0.8*x5;x1+2*x2+x4=100;2*x3+2*x4+x5=100;3*x1+x2+2*x3+3*x5=100;free(x1);free(x2);free(x3);free(x4);free(x5);1、DNA序列可以看成是由a、c、g、t四个字符按一定顺序排成的字符串，四个字符在DNA中出现的频率可以作为DNA的一个特征。统计文件“DNAData.txt”中各DNA序列中a、c、g、t各自出现的频率。如：对NDA序列:actgca,a出现的频率Pa=a出现的次数/DNA序列长度=2/6=0.33333。将统计结果使用xlswrite保存于result.xls中。 dna=textread(d:DNAData.txt,%s);P=zeros(40,4);for i=1:40 Na=sum(dnai=a); Nc=sum(dnai=c); Nt=sum(dnai=t); Ng=sum(dnai=g); n=Na+Nc+Nt+Ng; P(i,:)=Na Nc Nt Ng/n;enddlmwrite(d:NDAresult.txt,P,t)P=dlmread(d:NDAresult.txt)2、使用textread()函数读取数据“1.1 公汽线路信息.txt”，找出哪些线路通过S0819号站台。将结果使用函数dlmwrite保存于文件result.txt中。 C=textread(1.1 公汽线路信息.txt,%s);N=length(C);js=0;fid=fopen(result.txt,w);%打开文件result.txt，只写for i=1:N T=findstr(Ci,S0819);%查找是否包含S0819 if isempty(T)%若没找到 continue; else %找到 js=js+1; fprintf(fid,%sn,Ci); %将包含S0819的线路写入result.txt文件里 endendfclose(fid)3、使用load 命令导入数据 Indian_pines.mat，执行 imshow(data(:,:,116),) 查看结果，令I= data(:,:,116), 将I中的数据转化为0255之间的无符号8位(uint8)数据，并将其保存为图像 H.tif .load Indian_pines_corrected.matdata=Indian_pines_correctedimshow(data(:,:,151),)I= data(:,:,116);m=min(I(:);M=max(I(:);I1=uint8(1-m)/(M-m)*255);imwrite(I1,H.tif,tiff)figureimshow(I1)sets: SetA/i1.i5/; S