matlab实验程序Word文档格式.docx
- 文档编号:20951732
- 上传时间:2023-01-26
- 格式:DOCX
- 页数:15
- 大小:17.81KB
matlab实验程序Word文档格式.docx
《matlab实验程序Word文档格式.docx》由会员分享,可在线阅读,更多相关《matlab实验程序Word文档格式.docx(15页珍藏版)》请在冰豆网上搜索。
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:
end,2:
end);
whos;
T-4-1:
m=100:
999;
n=find(mod(m,21)==0);
length(n)
T-4-2:
ch='
ABC123d4e56Fg9'
;
k=find(ch>
='
A'
&
ch<
Z'
ch(k)=[];
ch
实验二
T-1:
E=eye(3);
R=rand(3,2);
O=zeros(2,3);
S=diag([2,3]);
%S=diag(diag(R))
A=[E,R;
O,S];
B=A^2
C=[E,(R+R*S);
O,S^2]
%B==C
T-2:
H=hilb(5)
P=pascal(5)
Hh=det(H)
Hp=det(P)
Th=cond(H)
Tp=cond(P)
T-3:
%A=fix(10*rand(5))
A=rand(5)
H=det(A)
Trace=trace(A)
Rank=rank(A)
Norm=norm(A)
T-4:
A=[-29,6,18;
20,5,12;
-8,8,5]
[V,D]=eig(A)
A*V
V*D
T-5a:
%
(1):
A=[1/2,1/3,1/4;
1/3,1/4,1/5;
1/4,1/5,1/6];
B1=[0.95,0.67,0.52]'
x1=inv(A)*B1
%
(2):
B2=[0.95,0.67,0.53]'
x2=inv(A)*B2
B3=[0.95,0.67,0.54]'
x3=inv(A)*B3
%(3):
cond(A)%矩阵A为病态矩阵
%C=A;
%C(9)=1/7;
%x3=inv(C)*B2
cond(A)%矩阵A为病态矩阵
T-5b:
A=hilb(4)
A(:
1)=[]
A(4,:
)=[]
B=[0.95,0.67,0.52]'
X=inv(A)*B
B1=[0.95,0.67,0.53]'
X1=inv(A)*B1
N=cond(B)
N1=cond(B1)
Na=cond(A)%矩阵A为病态矩阵
T-6:
A=[1,4,9;
16,25,36;
49,64,81]
B=sqrtm(A)
C=sqrt(A)%sqrtm函数是以矩阵为单位进行计算,sqrt函数是以矩阵中的元素进行计算
实验三
x=input('
请输入x的值'
)
ifx<
0&
x~=-3
y=x^2+x-6;
elseifx>
x<
5&
x~=2&
x~=3
y=x^2-5*x+6;
else
y=x^2-x-1;
end
y
T-1a:
x=[-5.0-3.01.02.02.53.05.0];
y=[];
fori=x
ifi<
i~=-3
y=[y,i^2+i-6];
elseifi>
i<
i~=2&
i~=3
y=[y,i^2-5*i+6];
y=[y,i^2-i-1];
T-2a:
请输入成绩'
ifx>
=90&
=100
disp('
等级为A'
=80&
=89
等级为B'
=70&
=79
等级为C'
=60&
=69
等级为D'
=59
等级为E'
成绩有误'
T-2b:
switch(x)
casenum2cell(90:
100)%case{90,91,92,93,94,95,96,97,98,99,100}
casenum2cell(80:
89)%
casenum2cell(70:
79)
casenum2cell(60:
69)
casenum2cell(0:
59)
otherwise
成绩错误'
n=input('
输入员工的工号:
'
t=input('
该工号员工的工作时数:
ift>
=120
s=t*84+(t-120)*84*0.15;
elseift<
=60
s=t*84-700;
s=t*84;
disp([num2str(n),'
号员工的应发工资为'
num2str(s)]);
T-4:
a=fix(10+90*rand(1,2))%产生两个随机整数
x=a
(1);
y=a
(2);
请输入运算符号:
'
s'
ift=='
+'
z=x+y;
elseift=='
-'
z=x-y;
*'
z=x*y;
/'
z=x/y;
disp([num2str(x),t,num2str(y),'
num2str(z)])%输出运算结果
T-5:
a=rand(5,6)%产生5x6的随机矩阵
请输入您要输出矩阵的第几行:
ifn>
5
超出了矩阵的行数,矩阵的最后一行为:
a(5,:
disp(['
矩阵的第'
num2str(n),'
行为:
])
a(n,:
实验四
T-1a:
s=0;
n=?
fori=1:
n
s=s+1/i/i;
PI=sqrt(6*s)
pi
T-1b:
a=1:
n;
b=1./a.^2;
PI=sqrt(6*sum(b))
n=1;
y=1;
while(y<
3)
n=n+1;
y=y+1/(2*n-1);
y=y-1/(2*n-1)
n=n-1
a=input('
a=?
b=input('
b=?
Xn=1;
Xn1=a/(b+Xn);
n=0;
whileabs(Xn1-Xn)>
1e-5
Xn=Xn1;
Xn1=a/(b+Xn);
ifn==500
break;
%跳出循环
end
Xn1
r1=(-b+sqrt(b*b+4*a))/2
r2=(-b-sqrt(b*b+4*a))/2
100
ifi==1
f(i)=1;
elseifi==2
f(i)=0;
elseifi==3
else
f(i)=f(i-1)-2*f(i-2)+f(i-3);
formatrat%转为有理格式
max(f)
min(f)
sum(f)
length(find(f>
0))
length(find(f==0))
length(find(f<
formatshort
T-4a:
f
(1)=1;
f
(2)=0;
f(3)=1;
fori=4:
T-5:
fori=2:
49%要求[2,50]区间
b=i*(i+1)-1;
%判断b是否是素数
m=fix(sqrt(b));
forj=2:
m
ifrem(b,j)==0%b不是素数
break
elseifj==m%b是素数
s=s+b;
s
实验五
%建立函数文件mat1.m
%functiony=mat1(x)
%y=[exp(x),log(x),sin(x),cos(x)];
%在命令窗口调用上述函数文件:
x=?
y=mat1(x)
%function[a1,a2,n1,n2]=mat2(m1,m2,th)
%t=th*pi/180;
%g=9.8;
%A=[m1*cos(t),-m1,-sin(t),0;
%m1*sin(t),0,cos(t),0;
%0,m2,-sin(t),0;
%0,0,-cos(t),1];
%B=[0;
m1*g;
0;
m2*g];
%C=inv(A)*B;
%a1=C
(1);
%a2=C
(2);
%n1=C(3);
%n2=C(4);
%在命令窗口调用该函数文件:
m1=input('
m1='
m2=input('
m2='
theta=input('
theta='
[a1,a2,N1,N2]=mat2(m1,m2,theta)
%判断x是不是素数,素数y=1;
否则y=0;
%functiony=mat3(x)
%y=1;
%fori=2:
sqrt(x)
%ifrem(x,i)==0
%y=0;
%break;
%end
%在命令窗口调用该函数文件:
fori=10:
99
j=10*rem(i,10)+fix(i/10);
ifmat3(i)&
mat3(j)
disp(i)
%functiony=fx(x)
%y=1./((x-2).^2+0.1)+1./((x-3).^4+0.01);
x='
y=fx(x)
T-5-1:
%functionf1=mat5_1(n)
%f1=n+10*log(n*n+5);
%在命令窗口中调用该函数文件:
y=mat5_1(40)/(mat5_1(30)+mat5_1(20))
T-5-2a:
%functionf2=mat5_2(n)
%f2=0;
%fori=1:
%f2=f2+i*(i+1);
%在命令窗口中调用该函数文件如:
y=mat5_2a(40)/(mat5_2a(30)+mat5_2a(20))
T-5-2b:
%functionf2=mat5_2b(n)
%i=1:
%m=i.*(i+1);
%f2=sum(m);
y=mat5_2b(40)/(mat5_2b(30)+mat5_2b(20))
实验六
x=linspace(0,2*pi,101);
y=(0.5+3*sin(x)./(1+x.^2)).*cos(x);
plot(x,y)
T-2-1:
x=linspace(-2*pi,2*pi,100);
y1=x.^2;
y2=cos(2*x);
y3=y1.*y2;
plot(x,y1,'
b-'
x,y2,'
r:
x,y3,'
g--'
text(4,16,'
\leftarrowy1=x^2'
text(6*pi/4,-1,'
\downarrowy2=cos(2*x)'
text(-1.5*pi,-2.25*pi*pi,'
\uparrowy3=y1*y2'
T-2-2:
subplot(1,3,1);
%分区
plot(x,y1);
title('
y1=x^2'
%设置标题
subplot(1,3,2);
plot(x,y2);
y2=cos(2*x)'
subplot(1,3,3);
plot(x,y3);
y3=x^2*cos(2*x)'
T-2-3:
x=linspace(-2*pi,2*pi,20);
subplot(2,2,1);
bar(x,y1);
y1=x^2的条形图'
subplot(2,2,2);
stairs(x,y1);
y1=x^2的阶梯图'
subplot(2,2,3);
stem(x,y1);
y1=x^2的杆图'
subplot(2,2,4);
fill(x,y1,'
r'
%如果少了'
则会出错
y1=x^2的填充图'
%其他的函数照样做。
x=-5:
0.01:
5;
%起始设y为空向量
forx0=x
ifx0<
=0
y=[y,(x0+sqrt(pi))/exp
(2)];
%将x对应的函数值放到y中
y=[y,0.5*log(x0+sqrt(1+x0^2))];
a='
b='
n='
theta=-2*pi:
2*pi;
rho=a*sin(b+n*theta);
polar(theta,rho)
x=linspace(-5,5,21);
y=linspace(0,10,31);
[x,y]=meshgrid(x,y);
%在[-5,5]*[0,10]的范围内生成网格坐标
z=cos(x).*cos(y).*exp(-sqrt(x.^2+y.^2)/4);
subplot(2,1,1);
surf(x,y,z);
subplot(2,1,2);
contour3(x,y,z,50);
%其中50为高度的等级数,越大越密
ezsurf('
cos(s)*cos(t)'
cos(s)*sin(t)'
sin(s)'
[0,0.5*pi,0,1.5*pi]);
%利用ezsurf隐函数
shadinginterp%进行插值着色处理
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- matlab 实验 程序