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 + b

A = [1 2 3;4 5 6;7 8 9];

B = [9 8 7;6 5 4;3 2 1];

A *

B A .* B

A =10;

B =20;

C =A/B

D =A\B

a = [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.^2

clear

X = [1 2;8 9;3 6];

X( : )

A = [ 1 2 3 ];

B = [ 4 5 6 ];

C = 3. ^A

D = A. ^B

5、t = -1:0.01:1;

y = (sqrt(3)/2).*exp((-4).*t).*sin(t.*4*sqrt(3)+(pi/3))

1、clc

clear

A = [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、clc

clear

A = [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、clc

clear

A = [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=10&A<25)

4、clc

clear

%使用函数实现矩阵左旋90°或右旋90°的功能。

A=magic(3), B=rot90(A),C=rot90(A,3) %B为左旋90°，C为右旋90°5、clc

clear

A = [1 2 2 2 2 2 2 2 2 2 2],

B = cumprod(A),%返回向量A的累乘积向量

C = cumsum(B),%返回向量B的累加和向量

S = max(C),%求和

6、clc

clear

A = [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.1

A{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.1

A{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.2

A{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.1

s=1;

for n=1:100

s=s*n;

if(s>10^100)

break;

end

end

a=sprintf('n=%d',n);

b=sprintf('%d!=%e>1e100',n,s);

disp(b)

1.2

n=1;

s=1;

while(s<=10^100)

s=s*n;

n=n+1;

end

a=sprintf('n=%d',n-1);

b=sprintf('n!=%e>1e100',s); disp(a)

disp(b)

2

f(1)=1;

f(2)=0;

f(3)=1;

for n=4:50

f(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))

3

a(1)=1;

a(2)=1;

for i=1:10000

a(i+2)=a(i)+a(i+1);

if a(i)>10000

break;

end

end

disp(['i=',num2str(i)])

disp(['a(i)=',num2str(a(i))])

1

x =

1 4 9 16 25 36 49 64 81 100 2

ans =

100

实验五

1.1

score=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 !');

end

end

1.2

score=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');

otherwise

disp('grade=E');

end

else

disp('输入的成绩不合理');

end

2

clear

clc

x=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;

otherwise

disp('输入的运算符有误');

break;

end

stc=sprintf('a=x%sy=%d',k,a)

3

clear

clc

x=[];

y=[];

x=input('请输入x的值：');

for k=1:length(x)

if x(k)<0 & x(k)~=-3

y(k)=x(k).^2+x(k)-6;

elseif x(k)>=0&x(k)<10&x(k)~=2&x(k)~=3

y(k)=x(k).^2-5*x(k)+6;

else y(k)=x(k).^2-x(k)-1;

end

end

for 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==xx

varargout{1}=a*b;

varargout{2}=a.*b;

else if x==xx&y==yy

varargout{1}=a.*b;

disp('只能作点乘')

varargout{2}=0;

else if x==yy&xx==y

disp('可以做两种矩阵乘')

varargout{1}=a*b;

varargout{2}=b*a;

else if x==yy||xx==y

try

varargout{1}=a*b;

disp('只能做矩阵乘')

varargout{2}=0;

catch

varargout{1}=b*a;

disp('只能做矩阵乘')

varargout{2}=0;

end

else

error('这两个矩阵不能乘')

end

end

end

end

2.

function y=qiuhe(a)

y=[];

[m,n]=size(a);

if m==1

s=0;

for k=1:n

s=s+a(1,k);

end

y=s;

else

for j=1:n

s=0;

for i=1:m

s=s+a(i,j);

end

y(1,j)=s;

end

end

3.

function s=jiecheng(n)

s=1;

if n==0

s=1;

else if n<0

error('n必须大于等于0') else

for i=1:n

s=s*i;

end

end

end

思考题：

function y=temf(a)

if length(a)<2

y=a;

else

y=max(a)-min(a);

end

实验七：

7.1

clear all;

clc;

qname=input('Enter file containing question:','s'); ip=fopen(qname,'rt');

if (ip<0)

disp('could not open input file')

end

op=fopen('qq.txt','wt');

if (op<0)

disp('could not open output file');

end

q=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');

end

fprintf(op,'%s\nAnswer:%s\n',q,a);

q=fgetl(ip);

end

fclose(ip);

fclose(op);

7.2

clear 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.3

a=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);

b

c

Position

实验八：

1.

clear,clc

x=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,clc

t=[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 on

4.

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.1

clear 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 on

title('二维图');

legend('y=exp(x)','y=log(x)');

xlabel('x轴坐标数据');

ylabel('y轴坐标数据');

1.2

clear 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.3

clear 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('火柴杆图');