重庆大学 数学实验 温罗生

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])
else for i=1:k(1)
D=a;
D(:,i)=b;
Di=det(D);
X(i)=Di/A;
end
end
X=X';
End
结果验证（如验证4阶）
a=rand(4);
b=[1;2;3;4];
X=cramer(a,b)
c=a*X
Ycramer
X =73.6009
6.6966
-69.4126
3.4100
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.0067 0.4229 0.0688 0.5309 0.3251 0.1537 0.4401 0.5181 0.2891 0.2548 0.3868 0.4709 0.3196 0.6544 0.5313 0.2665 0.4574 0.6377 0.6718 0.3445 0.6022 0.5985 0.6385 0.7184 0.6110 0.4235 0.5271 0.9436 0.6761 0.6678 0.6753 0.6959 0.6999 0.8200 0.9686 0.7788 0.8754 0.9577 0.6951 0.8444 0.7805
k=fix(i/3.3);
for l=k:j
h=[i,l,i/l];
if abs(h(3)-pi)<abs(H(3)-pi)
H=h;
end
end
end
f1=H(1);
f2=H(2);
f3=H(3);
end
结果：
[a,b,x]=PI(2)
a =22b =7x =3.1429
[a,b,x]=PI(3)
a =355b =113x =3.1416
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）相似的用法，表示求矩阵每一行或每一列的和；这几个函数在用法上大体一致，只是实现的功能有所不同而已。
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阶行列式进行验证，例如四阶行列式）
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
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
结果
y =
0.4229 0.6999 0.5309 0.9686 0.7788 0.2810 0.5181 0.6761 0.2548 0.7805 0.0942 0.6385 0.6544 0.5313 0.4235 0.4401 0.9436 0.2891 0.2240 0.6753 0.5985 0.0336 0.4076 0.3251 0.0908 0.5271 0.6377 0.6718 0.6678 0.0067 0.4709 0.0688 0.8200 0.1056 0.2665 0.4574 0.9577 0.6951 0.8444 0.6022 0.6959 0.3196 0.7184 0.6110 0.1537 0.8754 0.2407 0.0680 0.3445 0.3868
重庆大学
学生实验报告
实验课程名称数学实验
开课实验室DS1404
学院2014年级电子科学与技术专业班01
学生姓名学号
开课时间2016学年第2学期
总成绩
教师签名
数学与统计学院制
开课学院、实验室：DS1404实验时间：2016年3月10日
课程
名称
数学实验
实验项目
名称
实验项目类型
验证
演示
综合
设计
其他
指导
教师
c =
1.0000
2.0000
3.0000
4.0000
总结与体会
第一次写实验报告是初次将学到的东西进行运用，也将平时老师讲的一些知识点进行连接，实验有助于巩固自己学到的知识。一二题是老师上课讲的一些比较基础的知识；三题需要用到矩阵的知识；四五题则是相应的作图知识，尤其是在四小问中的作图不同于老师上课讲的一些基本的方式，而用了ezplot做出了图形；后面三道题就是M-文件的练习求解，包括如何建立函数式M-文件等等；在做题中我了解了程序设计中的一些基本思想，代码的优化；但是感觉在第一次做起来还是比较吃力，但是我认为将老师上课讲的一些东西，知识点必须要吃透，其次最重要的是一步一步多加练习，总结代码编写的一些经验。
n=n+1;
end
end
end
end
H
结果：
H =
153 370 371 407
7.代码:
(1)
function 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
[a,b,x]=PI(4)
a =1065b =339x =3.1416
说明:[a,b,x]=PI(y),其中a代表分子，b代表分母，y为位数
附加：验证克莱姆法则（n阶）
程序：
function X=cramer(a,b)
A=det(a);
k=size(a);
if A==0
disp('输入有错，解不唯一。');
图形
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;
f=judge(i);
if f==1
n=n+1;
H(n)=i;
end
end
H
L=size(H)
hist(H,10)
等距区间内质数的个数的直方图（以1000为间距）
8.代码：
function [f1,f2,f3]=PI(x)
h=[];
H=[3,1,3];
for i=10^(x-1):10^x-1
j=fix(i/3);
结果：
[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);
温老师
成绩
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));