重庆大学--数学实验-温罗生

ans =
0.6959 0.6999 0.8200 0.9686 0.7788 0.8754 0.9577 0.6951 0.8444 0.7805
ans =
0.9686
0.9436
0.6718
0.9577
0.8754
ans =
0.0942 0.0336 0.4076 0.1056 0.0908 0.2810 0.2407 0.0680 0.2240 0.0067
ans =
0.9686
0.9436
0.6718
0.9577
0.8754
ans =
2.2825 1.7604 3.1313 2.5417 1.7132 2.5810 3.2978 2.4001 2.3355 2.4515
ans =
0.0942 0.0336 0.4076 0.1056 0.0908 0.2810 0.2407 0.0680 0.2240 0.00670.4229 0.0688 0.5309 0.3251 0.1537 0.4401 0.5181 0.2891 0.2548 0.38680.4709 0.3196 0.6544 0.5313 0.2665 0.4574 0.6377 0.6718 0.3445 0.60220.5985 0.6385 0.7184 0.6110 0.4235 0.5271 0.9436 0.6761 0.6678 0.67530.6959 0.6999 0.8200 0.9686 0.7788 0.8754 0.9577 0.6951 0.8444 0.7805
图形
5.代码
[x,y]=meshgrid(-5:0.1:5);
z=sin(pi*sqrt(x.^2+4*y.^2));
mesh(x,y,z)
图形
6.代码
H=[];
n=1;
for a=1:9
for b=0:9
for c=0:9
h=100*a+10*b+c;
j=a^3+b^3+c^3;
if h==j
H(n)=j;
n=n+1;
end
end
end
end
H
结果：
H =
153 370 371 407
7.代码:(1)f来自nction f=judge(x)
f=1;
if x<4
else for i=2:x-1
k=mod(x,i);
if k==0
f=0;
end
end
end
end
(2)
n=0;
H=[];
for i=2:20000
f=judge(i);
结果
y =
0.4229 0.6999 0.5309 0.9686 0.7788 0.2810 0.5181 0.6761 0.2548 0.78050.0942 0.6385 0.6544 0.5313 0.4235 0.4401 0.9436 0.2891 0.2240 0.67530.5985 0.0336 0.4076 0.3251 0.0908 0.5271 0.6377 0.6718 0.6678 0.00670.4709 0.0688 0.8200 0.1056 0.2665 0.4574 0.9577 0.6951 0.8444 0.60220.6959 0.3196 0.7184 0.6110 0.1537 0.8754 0.2407 0.0680 0.3445 0.3868
2.代码
y=rand(5,10)
max(y),max(y,[],2)
min(y),max(y,[],2)
sum(y)
sort(y)
length(y)
size(y)
说明：rand(x)随机产生矩阵，max(y)求取矩阵y每一行最大值,max(y,[],2)求取矩阵y每一列最大值；其中min（y）类似。sum(y)相应的与max（y）相似的用法，表示求矩阵每一行或每一列的和；这几个函数在用法上大体一致，只是实现的功能有所不同而已。
重庆大学
学生实验报告
实验课程名称数学实验
开课实验室DS1404
学院2014年级电子科学与技术专业班01
学生姓名学号
开课时间2016学年第2学期
总成绩
教师签名
数学与统计学院制
开课学院、实验室：DS1404实验时间：2016年3月10日
课程
名称
数学实验
实验项目
名称
实验项目类型
验证
演示
综合
设计
其他
指导
教师
温老师
成绩
1.代码
x=3.14;
y1=sign(x);
y2=x^3;
y3=sin(x);
y4=cos(x);
y5=tan(x^2);
y6=2*tan(x);
y7=3^x;
y8=log(abs(x-x^x));
y9=log2(x)/log2(3);
y10=fix(x);
y11=asin(x^2-3*pi);
y12=x^(1/6);
y13=sqrt(2*x^2+(4*x/(log(x+exp(x)*sin(x))/2*x^2)));
y=[y1,y2,y3,y4,y5,y6,y7,y8,y9,y10,y11,y12,y13];
y
结果
y =
1.0000 30.9591 0.0016 -1.0000 0.4645 -0.0032 31.4891 3.5025 1.0415 3.0000 0.4498 1.2101 4.6822
ans =
10
3.代码（n阶行列式）
function [u,v]=DET(x,y)
A=rand(x);
u=det(A);
v=0;
for i=1:x
a=A;
a(y,:)=[];
a(:,i)=[];
v=A(y,i)*det(a)*(-1)^(y+i)+v;
end
end
（下面任意选取n阶行列式进行验证，例如四阶行列式）
plot(theta,rou);
t=-10:0.1:10;
y1=(3*t)./(1+t.^3);
y2=(3*t.^2)./(1+t.^3);
subplot(2,2,3);
plot(y1,y2);
subplot(2,2,4);
ezplot('9*(x+2)^2+4*(y-6)^2-36',[-6,2,-4,10])
结果：
[u,v]=DET(4,2)
u = -0.0261
v =-0.0261
4.代码
x=-10:0.1:10;
y=exp(-((x-2).^2)/4);
subplot(2,2,1);
plot(x,y);
theta=-pi:0.1:pi;
rou=sin(3.*theta);
subplot(2,2,2);