重 庆 大 学
学 生 实 验 报 告
实验课程名称 数学实验
开课实验室 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 = 3.5025 1.0415 3.0000 0.4498 1.2101 4.6822 2.代码 y=rand(5,10) max(y),max(y,[],2) min(y),max(y,[],2) 1.0000 30.9591 0.0016 -1.0000 0.4645 -0.0032 31.489sum(y) sort(y) length(y) size(y) 说明：rand(x)随机产生矩阵，max(y)求取矩阵y每一行最大值,max(y,[],2) 求取矩阵y每一列列的和；这几个函数在用法上大体一致，只是实现的功能有所不同而已。 最大值；其中min（y）类似。sum(y)相应的与max（y）相似的用法，表示求矩阵每一行或每一结果 y = 0.4229 0.6999 0.5309 0.9686 0.7788 0.2810 0.5181 0.6760.2548 0.7805 0.0942 0.6385 0.6544 0.5313 0.4235 0.4400.9436 0.2891 0.2240 0.6753 0.5985 0.0336 0.4076 0.3250.0908 0.5271 0.6377 0.6718 0.6678 0.0067 0.4709 0.0680.8200 0.1056 0.2665 0.4574 0.9577 0.6951 0.8440.2407 0.0680 0.3445 0.3868 ans = 0.8444 0.7805 ans = 0.9686 0.9436 0.6718 0.9577 0.8754 ans = 0.2240 0.0067 ans = 0.6022 0.6959 0.3196 0.7184 0.6110 0.1537 0.8750.6959 0.6999 0.8200 0.9686 0.7788 0.8754 0.9577 0.6950.0942 0.0336 0.4076 0.1056 0.0908 0.2810 0.2407 0.0680.9686 0.9436 0.6718 0.9577 0.8754 ans = 2.3355 2.4515 ans = 2.2825 1.7604 3.1313 2.5417 1.7132 2.5810 3.2978 2.4000.0942 0.0336 0.4076 0.1056 0.0908 0.2810 0.2407 0.0680.2240 0.0067 0.4229 0.0688 0.5309 0.3251 0.1530.4401 0.5181 0.2891 0.2548 0.3868 0.4709 0.3196 0.6540.5313 0.2665 0.4574 0.6377 0.6718 0.3445 0.6020.5985 0.6385 0.7184 0.6110 0.4235 0.5271 0.9436 0.6760.8754 0.9577 0.6951 0.8444 0.7805 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; 0.6678 0.6753 0.6959 0.6999 0.8200 0.9686 0.778end end （下面任意选取n阶行列式进行验证，例如四阶行列式） 结果： [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); 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])