MATLAB实验代码汇总

2、 3vars global help My_exp sin X+Y _input E-4 Abcd AB_C_ ;3、 A = [ 1 2 3;4 5 6;7 8 9]y = [12 + 2 * (7 -4)]/(3^2)4、 A = 15;B = 20;C = A + B c = a + bA = [1 2 3;4 5 6;7 8 9];B = [9 8 7;6 5 4;3 2 1];A *B A .* BA =10;B =20;C =A/BD =A\Ba = [1 -2 3;4 5 -4;5 -6 7]A = [1,2;3,4] + i * [5,6;7,8]A = [1,2;3,4] + i[5,6;7,8]a = [1 2 3;3 4 2;5 2 3];a^2 a.^2clearX = [1 2;8 9;3 6];X( : )A = [ 1 2 3 ];B = [ 4 5 6 ];C = 3. ^AD = A. ^B5、t = -1:0.01:1;y = (sqrt(3)/2).*exp((-4).*t).*sin(t.*4*sqrt(3)+(pi/3))1、clcclearA = [12 34 -4;34 7 87;3 65 7];B = [1 3 -1;2 0 3;3 -2 7];a=A+6*B,I=eye(3),b=A^2-B+I,c=A*B,d=A.*B,e=A^3,f= A.^3,g=[A,B],h= [(A ([1,3],:));B^2]2、clcclearA = [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 * B ,D = C([3 4 5],[2 3])3、clcclearA = [23 10 -0.778 0; 41 -45 65 5; 32 5 0 32; 6 -9.54 54 3.14];B = A([1:3],:);C = A(:,[1:2]);D = A([2 3 4],[1 2]);E = B * C;a= E<D ,b= E&D,c= E|D ,d=~E|~D,e=(A>=10&A<25)4、clcclear%使用函数实现矩阵左旋90°或右旋90°的功能。

A=magic(3), B=rot90(A),C=rot90(A,3) %B为左旋90°，C为右旋90°5、clcclearA = [1 2 2 2 2 2 2 2 2 2 2],B = cumprod(A),%返回向量A的累乘积向量C = cumsum(B),%返回向量B的累加和向量S = max(C),%求和6、clcclearA = [1 10 11 16 21 2 9 12 17 22 3 8 13 18 23 4 7 14 19 24 5 6 15 20 25];B = reshape(A,5,5)C1 = B(2,[2 3])C2 = B([2 3 4],4)C3 = B([4 5],[1 2 3])实验三1.a=reshape(100:999,30,30);b=rem(a,21) %求余c=find(b==0)length(c) %获取向量长度，若输入参数为矩阵或多维数组，则返回各个维尺寸的最大值2.ch='FrIeNd123'd=find(ch>='A'&ch<='Z')ch(d)=[]3.%若需要在字符串内容中包含单引号，则在键入字符串内容时，连续键入两个单引号即可ch='In GB usage quotation marks are usually singel:''fire!'''5.1A{1,1}=1;A{1,2}='Brenden';A{2,1}=reshape(1:9,3,3);A{2,2}={12,34,2;54,21,3;4,23,67};upper(A{1,2}) %将字符串的字符都转变成为大写字符5.2.1A{1,1}=1;A{1,2}='Brenden';A{2,1}=reshape(1:9,3,3);A{2,2}={12,34,2;54,21,3;4,23,67};A(3)=[] %使用圆括号（）直接访问元胞数组的元胞，获取的数据类型是一个元胞5.2.2A{1,1}=1;A{1,2}='Brenden';A{2,1}=reshape(1:9,3,3);A{2,2}={12,34,2;54,21,3;4,23,67};A{3}=[] %使用大括号{}直接访问元胞数组的元胞，获取的数据类型是元胞的数据类型6.a=rand(4,3) %产生均匀分布的随机数4*3矩阵，数值范围（0,1）b=mean(a) %返回一个行向量，其第i个元素是矩阵a的第i列的算术平均值c=median(a) %返回一个行向量，其第i个元素是a的第i列的中值b1=mean(b)c1=median(c)d=max(a)d1=max(d) %求最大值e=min(a)e1=min(e) %求最小值[i1,j1]=find(a==d1) %最大值元素的位置[i2,j2]=find(a==e1) %最小值元素的位置t=sum(a')x=sum(a)y=sum(x) %求和sort(a) %返回一个对A中的每列元素按升序排列的新矩阵sort(a,2,'descend') %返回一个对A中的每列元素按降序排列的新向量降序实验四1.1s=1;for n=1:100s=s*n;if(s>10^100)break;endenda=sprintf('n=%d',n);b=sprintf('%d!=%e>1e100',n,s);disp(b)1.2n=1;s=1;while(s<=10^100)s=s*n;n=n+1;enda=sprintf('n=%d',n-1);b=sprintf('n!=%e>1e100',s); disp(a)disp(b)2f(1)=1;f(2)=0;f(3)=1;for n=4:50f(n)=f(n-1)-2*f(n-2)+f(n-3); end[MAX,x]=max(f)[MIN,y]=min(f)SUM=sum(f)positive=length(find(f>0)) zero=length(find(f==0)) negative=length(find(f<0))3a(1)=1;a(2)=1;for i=1:10000a(i+2)=a(i)+a(i+1);if a(i)>10000break;endenddisp(['i=',num2str(i)])disp(['a(i)=',num2str(a(i))])1x =1 4 9 16 25 36 49 64 81 100 2ans =100实验五1.1score=input('please input grade=');if(score>100|score<0|rem(score,0.5)~=0)disp('输入的成绩不合理');elseif(score>=90&score<=100)disp('the grade = A !');elseif(score>=80&score<=89)disp('the grade = B !');elseif(score>=70&score<=79)disp('the grade = C !');elseif(score>=60&score<=69)disp('the grade = D !');else if(score<60)disp('the grade = E !');endend1.2score=input('please input grade=');if (score<=100&score>=0&rem(score,0.5)==0) switch fix(score/10)case{9,10}disp('grade=A')case{8}disp('grade=A');case{7}disp('grade=C');case{6}disp('grade=D');otherwisedisp('grade=E');endelsedisp('输入的成绩不合理');end2clearclcx=fix(50*(rand(1)-rand(1)))y=fix(50*(rand(1)-rand(1)))k=input('请输入一个运算符号：','s');switch(k)case{'+'}a=x+y;case{'-'}a=x-y;case{'*'}a=x*y;case{'/'}a=x/y;otherwisedisp('输入的运算符有误');break;endstc=sprintf('a=x%sy=%d',k,a)3clearclcx=[];y=[];x=input('请输入x的值：');for k=1:length(x)if x(k)<0 & x(k)~=-3y(k)=x(k).^2+x(k)-6;elseif x(k)>=0&x(k)<10&x(k)~=2&x(k)~=3y(k)=x(k).^2-5*x(k)+6;else y(k)=x(k).^2-x(k)-1;endendfor i=1:length(y)disp(['y(',num2str(x(i)),')','=',num2str(y(i))]);; end实验六：1.function varargout=chengji(a,b)[x,y]=size(a);[xx,yy]=size(b);if xx==yy&x==y&x==xxvarargout{1}=a*b;varargout{2}=a.*b;else if x==xx&y==yyvarargout{1}=a.*b;disp('只能作点乘')varargout{2}=0;else if x==yy&xx==ydisp('可以做两种矩阵乘')varargout{1}=a*b;varargout{2}=b*a;else if x==yy||xx==ytryvarargout{1}=a*b;disp('只能做矩阵乘')varargout{2}=0;catchvarargout{1}=b*a;disp('只能做矩阵乘')varargout{2}=0;endelseerror('这两个矩阵不能乘')endendendend2.function y=qiuhe(a)y=[];[m,n]=size(a);if m==1s=0;for k=1:ns=s+a(1,k);endy=s;elsefor j=1:ns=0;for i=1:ms=s+a(i,j);endy(1,j)=s;endend3.function s=jiecheng(n)s=1;if n==0s=1;else if n<0error('n必须大于等于0') elsefor i=1:ns=s*i;endendend思考题：function y=temf(a)if length(a)<2y=a;elsey=max(a)-min(a);end实验七：7.1clear all;clc;qname=input('Enter file containing question:','s'); ip=fopen(qname,'rt');if (ip<0)disp('could not open input file')endop=fopen('qq.txt','wt');if (op<0)disp('could not open output file');endq=fgetl(ip);while(ischar(q))fprintf('%s\n',q)a=input('Answer Y(yes) or N(not):','s');while((a~='Y')&&(a~='N'))a=input('Answer Y(yes) or N(not):','s');endfprintf(op,'%s\nAnswer:%s\n',q,a);q=fgetl(ip);endfclose(ip);fclose(op);7.2clear all;clc;f=fopen('a.txt');h=fopen('aa.txt','w');fwrite(h,upper(char(fread(f,inf,'char'))),'char');fclose(f);fclose(h);7.3a=1:5;fid=fopen('fdat.bin','w');fwrite(fid,a,'int16');status=fclose(fid);fid=fopen('fdat.bin','r');status=fseek(fid,6,'bof');b=fread(fid,1,'int16');position=ftell(fid);status=fseek(fid,-4,'cof');c=fread(fid,1,'int16');status=fclose(fid);bcPosition实验八：1.clear,clcx=0:1:4*pi;y=sin(x).*exp(-x/5);xi=0:0.3:4*pi;y1=interp1(x,y,xi,'nearset');y2=interp1(x,y,xi,'linear');y3=interp1(x,y,xi,'spline');y4=interp1(x,y,xi,'cubic');plot(x,y,'-m',xi,y1,'*r',xi,y2,'ok',xi,y3,'^b',xi,y4,'+g') legend('原始数据','最近点插值','线性插值','样条插值','立方插值')title('第一题四种插值方法进行插值','fontname','隶书','fontsize',20)xlabel('x轴','fontname','隶书','fontsize',16)ylabel('y轴','fontname','隶书','fontsize',16)2.clear,clct=[13 15 17 19 20];sex=[0,1]';y=[270 290 330 350 380;250 273 310 330 360];x=[13:20];sex_2=[0,1]';y_2=interp2(t,sex,y,x,sex_2)surf(x,sex_2,y_2)3.x=2007:2013;y=[70 122 144 152 174 196 202];k=polyfit(x,y,1);x1=2007:2016;y1=polyval(k,x1);plot(x,y,'-ro',x1,y1,'-b*');title('某乡镇企业2004-2013年的生产利润')legend('07-13年的利润','07-16年的利润趋势线')xlabel('年份');ylabel('利润');text(2015,polyval(k,2015),'<---2015年利润')text(2016,polyval(k,2016),'<---2016年利润')grid on4.clear all ;clc;load('census.mat');nian=[1990:2010];p=polyfit(cdate,pop,2);z=polyval(p,cdate);zz=polyval(p,nian);my=mean(pop);mz=mean(z);R=sum((pop-my).*(z-mz))./sqrt(sum((pop-my).^2).*sum((z-mz).^2))plot(cdate,pop,'*',cdate,z,'r',nian,zz,'b-')xlabel('年份（十年）','fontname','隶书','fontsize',16) ylabel('人口（百万)','fontname','隶书','fontsize',16)title('第四题美国自1790-1990年的总人口','fontname','隶书','fontsize',20)text(1990,248.73,'<---之后为拟合数据')legend('原数据','拟合曲线')实验九1.1clear all;clc;x1=-2:0.1:2;y1=exp(x1);x2=0:0.1:5;y2=log(x2);plot(x1,y1,'-.r',x2,y2,'--b');grid ontitle('二维图');legend('y=exp(x)','y=log(x)');xlabel('x轴坐标数据');ylabel('y轴坐标数据');1.2clear all;clc;x=0:pi/100:pi/2;y=tan(x);plot(x,y,'-+r');title('正切函数');legend('y=tan(x)');xlabel('x轴数据');ylabel('y轴数据');axis([0,pi/2,0,5]);1.3clear all;clc;t=-pi:pi/20:pi;x=t.*cos(3*t);y=t.*(sin(t)).^2;plot(x,y);legend(strvcat('x=tcos(3t)','y=tsin^2t'));xlabel('x轴数据');ylabel('y轴数据');2.clear all;clc;[x,y]=meshgrid(-1.5*pi:1.5*pi);z=cos(x).*sin(y)./y;subplot(1,2,1);surf(x,y,z)title('曲面图');xlabel('x轴数据');ylabel('y轴数据');zlabel('z轴数据');subplot(1,2,2);mesh(x,y,z);title('线框图');xlabel('x轴数据');ylabel('y轴数据');zlabel('z轴数据');3.clear all;clc;data=[10 2 3 5;5 8 10 3;9 7 6 1;3 5 7 2;4 7 5 3]; a=sum(rand(5,5));theta=2*pi*(0:127)/128;x=cos(theta);y=sin(theta);z=abs(fft(ones(10,1),128))';subplot(2,2,1);bar(data);grid on;title('条状图');subplot(2,2,2);area(data);grid on;title('面积图');subplot(2,2,3);pie(a);grid on;title('饼图');subplot(2,2,4);stem3(x,y,z);grid on; title('火柴杆图');。