北京科技大学数学实验第一次作业

《数学实验》报告

实验名称

学院

专业班级

姓名

学号

2014年11月

一、【实验目的】

1. 比较简单的了解MATLAB软件的基本用途；

2. 了解并掌握MATLAB的工作环境；

3. 掌握MATLAB的基本操作及常用命令；

4. 掌握MATLAB数据显示格式的调节方法；

5. 掌握MATLAB中的常用函数与变量、表达式的定义方法；

6. 掌握MATLAB常用的工作方式M文件的编程工作方式；

7. 掌握MATLAB中矩阵及数组的有关知识；

8. 掌握MATLAB语言中的程序结构。

二、【实验任务】

完成课后第一章第4题（p16），第二章第1、2、3题（p27），第三章第1题（p34）。

三、【实验程序】

P16 第4题

程序：

function jiecheng(n)

y=1;sum=0;

for i=1:n

y=y*i

sum=sum+y;

end

sum

>> jiecheng(20)

sum =

2.5613e+18

P27第1题

程序：

>> format short,exp(1)

ans =

2.7183

>> format long,exp(1)

2.718281828459046

>> format short e,exp(1)

ans =

2.7183e+00

>> format long e,exp(1)

ans =

2.718281828459046e+00 >> format short g,exp(1)

ans =

2.7183

>> format long g,exp(1)

ans =

2.71828182845905 >> format rat,exp(1)

ans =

1457/536

>> format hex,exp(1)

ans =

4005bf0a8b14576a

>> format bank,exp(1)

ans =

2.72

>>

P27第2题

A*B 是A矩阵与B矩阵相乘，A.*B是A矩阵中的每一个数与B矩阵中对应数相乘得到新的一个矩阵。

程序：

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

A =

1 2 3

4 5 6

7 8 9

>> B=[4 6 8;5 5 6;3 2 2]

B =

4 6 8

5 5 6

3 2 2

>> A*B

ans =

23 22 26

59 61 74

95 100 122

>> A.*B

ans =

4 12 24

20 25 36

21 16 18

>>

P27第3题程序：

>> A=[5 2;9 1]

A =

5 2

9 1 >> B=[1 2;9 2]

B =

1 2

9 2 >> A>B

ans =

1 0

0 0 >> A==B

ans =

0 1

1 0 >> A

ans =

0 0

0 1 >> (A==B)&(A

0 0

>> (A==B)&(A>B)

ans =

0 0

0 0

>>

P34第1题

程序：

function PAI(n)

n=1;b=0;a=0;

while (1/(2*n-1))>=0.000001 a=1/(2*n-1);

b=b+((-1)^(n+1))*a;

n=n+1;

end

4*b

>> PAI

ans =

3.14159065358969

>>

四、【实验结果】

P16 第4题

P27第1题

P27第2题

P27

第3题

P34第

1题

五、【实验总结】

最新版北京科技大学第三次数学实验报告

《数学实验》报告 实验名称Matlab三维曲面绘图 学院东凌经济管理学院 专业班级 姓名 学号 2016年3月

一、【实验目的】 1.了解并掌握Matlab三维曲面绘图； 2.进一步掌握绘图程序格式和意义； 3.初步掌握meshgrid, mesh, surf, colordef, colormap, light等使用。 二、【实验任务】 79-7 79-9 三、【实验程序】 79-7 t1=-3:0.1:3; [x1,y1]=meshgrid(t1); z1=x1.^2+y1.^2;

subplot(1,2,1);colordef white;light('position',[20,20,5]);colormap(pin k); mesh(x1,y1,z1),title('x^2+3.*y^2'); subplot(1,2,2);colordef white;light('position',[20,20,5]);colormap(pin k); surf(x1,y1,z1),title('x^2+3.*y^2') 79-9 t=-2:0.1:2; [x,y]=meshgrid(t); z1=5-x.^2-y.^2; subplot(1,3,1),mesh(x,y,z1),title('抛物面') z2=3*ones(size(x)); subplot(1,3,2),mesh(x,y,z2),title('平面') r0=abs(z1-z2)<=0.2; zz=r0.*z2;yy=r0.*y;xx=r0.*x; subplot(1,3,3),plot3(xx,yy,zz,'x'),title('交线') 四、【实验结果】 79-1

MATLAB实验练习题(计算机)-南邮-MATLAB-数学实验大作业答案

“”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> ('(x)-3*x^2',0) = -2*(-1/6*3^(1/2)) -2*(-11/6*3^(1/2)) -2*(1/6*3^(1/2)) 3、求解下列各题： 1）30 sin lim x x x x ->- >> x;

>> (((x))^3) = 1/6 2） (10)cos ,x y e x y =求 >> x; >> ((x)*(x),10) = (-32)*(x)*(x) 3）2 1/2 0(17x e dx ?精确到位有效数字） >> x; >> ((((x^2),0,1/2)),17) =

0.54498710418362222 4）4 2 254x dx x +? >> x; >> (x^4/(25^2)) = 125*(5) - 25*x + x^3/3 5）求由参数方程arctan x y t ??=? =??dy dx 与二阶导 数22 d y dx 。 >> t; >> ((1^2))(t); >> ()() = 1

6）设函数(x)由方程e所确定，求y′(x)。>> x y; *(y)(1); >> ()() = (x + (y)) 7） sin2 x e xdx +∞- ? >> x; >> ()*(2*x); >> (y,0) = 2/5

8） 08x =展开（最高次幂为） >> x (1); taylor(f,0,9) = - (429*x^8)/32768 + (33*x^7)/2048 - (21*x^6)/1024 + (7*x^5)/256 - (5*x^4)/128 + x^3/16 - x^2/8 + 2 + 1 9） 1sin (3)(2)x y e y =求 >> x y; >> ((1)); >> ((y,3),2) =

北京科技大学工科物理实验 下

全息照相1、画出全息照相光路图 P140图 2、全息照相的影响因素 a、系统的稳定性及实验环境是否符合要求 b、参考光与物光的光程差 c、参考光与物光的夹角（30-50）为宜 d、参考光与物光的光强比（3:1-5:1）为宜 e、底片的安放 f、曝光与显影时间的控制（显影1-2分颜色变，定影2-4分，水洗2-3分） 3、全息照相和普通照相的区别 P138 普通照相是用几何光写的方法记录物体上各点的发光强度分布，得到的是二维平面像，像上各点的照度与物体上各点的发光强度一一对应。而全息照相的记录对象是整个物体发出的光波，借助于参考光用干涉的方法记录这个物光波的振幅和位相分布，即记录下物光波与参考光波相干后的全部信息。 4、全息照相与普通照相的本质区别 答：全息照相和常规照相之不同在于，常规照相只是记录了被摄物体表面光线强弱的变化，即只记录了光的振幅；而全息照相则记录了光波的全部信息，除振幅外，还记录了光波的相位。这样就把空间物体光波场的全部信息都贮存记录了下来。然后利用全息照片对特定波长单色照明光的衍射，把原空间景象显现出来。

5、全息照相实验失败率很高，怎样改进 调节物光束与参考光束由分光板至感光板的光程大致相等，用透镜将物光束扩展到整个被摄物都受到均匀的光照，物光束与参考光束的夹角应为30到50度，调整参考光束与物光束的光强比，关闭室内照明灯减少其他光的干扰，控制曝光时间。 6、做全息照相时，能否用一张底片记录两个物体的像 7、全息照相中物象光和参考光的夹角小于多少度 照在感光板上的物光束和参考光束的夹角在30°~40°为宜 8、全息照相的特点 8、全息照相原理 9、全息照相利用光的干涉，将物体的（振幅）和（相位）特性以（干涉条纹）形式表现出来。 10、影响全息照相成败的条件 11、扩光束的作用是什么分光束分得的光，光强较大的是物光还是参考光，为什么 12、全息扩光镜的作用（两个分别写） 13、全息照相中的感光片是否有正反面？ 有正反面，有感光乳胶的面为正面。 传感器 1、传感器的定义 传感器是一种能把特定的被测量信息（包括物理量、化学量、生物学量等）按一定规律转换成某种可用信号输出的器件或装置。 2、传感器通常由哪两部分组成各部分的作用是什么 传感器通常由敏感元件和传感元件组成。

大一数学实验

2017春季数学实验报告 班级：计算机系61 姓名：赵森学号：2160500026（校内赛编号506）班级：计算机系61 姓名：冯丹妮学号：2160500002（校内赛编号327）班级：计算机系63 姓名：郝泽霖学号：2160500054

第一次上机作业 实验8: 练习1： 4．某棉纺厂的原棉需从仓库运送到各车间。各车间原棉需求量、单位产品从各仓库运往各车间的运输费以及各仓库的库存容量如表8.5所列，问如何安排运输任务使得总运费最小？ 设仓库1运往车间1,2,3,的原棉量为x1,x2,x3, 仓库2运往车间1,2,3,的原棉量为x4,x5,x6, 仓库3运往车间1,2,3,的原棉量为x7,x8,x9。 2x1+x2+3x3<=50 2x4+2x5+4x6<=30 3x7+4x8+2x9<=10 X1+x4+x7=40 X2+x5+x8=15 X3+x6+x9=35 程序： c=[2,1,3,2,2,4,3,4,2]; a(1,:)=[1,1,1,0,0,0,0,0,0]; a(2,:)=[0,0,0,1,1,1,0,0,0]; a(3,:)=[0,0,0,0,0,0,1,1,1]; aeq(1,:)=[1,0,0,1,0,0,1,0,0]; aeq(2,:)=[0,1,0,0,1,0,0,1,0]; aeq(3,:)=[0,0,1,0,0,1,0,0,1]; b=[50;30;10]; beq=[40;15;35]; vub=[]; vlb=zeros(9,1); [x,fval]=linprog(c,a,b,aeq,beq,vlb,vub) 结果： x = 10.0000 15.0000 25.0000

c++大作业学生实验报告

学生实验报告 实验课名称： C++程序设计 实验项目名称：综合大作业——学生成绩管理系统专业名称：电子信息工程 班级： 学号： 学生： 同组成员： 教师：

2011 年 6 月 23 日 题目：学生成绩管理系统 一、实验目的： （1）对C++语法、基础知识进行综合的复习。 （2）对C++语法、基础知识和编程技巧进行综合运用，编写具有一定综合应用价值的稍大一些的程序。培养学生分析和解决实际问题的能力，增强学生的自信心，提高学生学习专业课程的兴趣。 （3）熟悉掌握C++的语法和面向对象程序设计方法。 （4）培养学生的逻辑思维能力，编程能力和程序调试能力以及工程项目分析和管理能力。 二、设计任务与要求： （1）只能使用/C++语言，源程序要有适当的注释，使程序容易阅读。 （2）至少采用文本菜单界面（如果能采用图形菜单界面更好）。 （3）要求划分功能模块，各个功能分别使用函数来完成。 三、系统需求分析： 1.需求分析: 为了解决学生成绩管理过程中的一些简单问题，方便对学生成绩的管理 （录入，输出，查找，增加，删除，修改。） 系统功能分析: (1):学生成绩的基本信息：学号、、性别、C++成绩、数学成绩、英语成绩、 总分。 (2):具有录入信息、输出信息、查找信息、增加信息、删除信息、修改信息、 排序等功能。 2.系统功能模块（要求介绍各功能） （1）录入信息(Input)：录入学生的信息。 （2）输出信息(Print)：输出新录入的学生信息。 （3）查找信息(Find)：查找已录入的学生信息。 （4）增加信息(Add)：增加学生信息。 （5）删除信息(Remove)：在查找到所要删除的学生成绩信息后进行删除并输出删除后其余信息。 （6）修改信息(Modify)：在查到所要修改的学生信息后重新输入新的学生信息从而进行修改，然后输出修改后的所有信息。 （7）排序(Sort)：按照学生学号进行排序。 3.模块功能框架图

安培力实验报告(北京科技大学物理实验报告)

北京科技大学实验报告 安培力 实验目的 学会设计简单的实验方案，利用自组仪器合理搭建实验设备；学会用归纳法研究磁场中载流导线的受力规律。实验原理 将一段通电导线置于以均匀磁场中，导线受到的安培力与磁场的方向和强弱、载流方向和强度大小、导体形状和尺寸相关。 实验仪器 天平LGN310、矩形磁极一对、电磁线圈两个、U形铁芯一个、4个尺寸不同的导体、直流电源、电流表、轻质金属导线、开关、导线等。 实验内容及步骤 （1）利用电流秤测量磁场中载流导体的受力。 （2）利用U形铁芯的两个电磁磁极产生较强的匀强磁场。 （3）调节砝码使天平达到平衡。 （4）改变励磁电流、载流线圈电流、载流线圈尺寸测量线圈所受安培力大小。以研究安培力与磁场、载流、导线尺寸的关系。 数据测量结果：

数据分析： 此次实验是为了研究载流导体在磁场中的受力规律的。实验中主要研究安培力大小与励磁电流、载流线圈电流、载流线圈尺寸的关系。实验数据中的“-”号表示反向测量时的电流。并且，实验中以安培力向下方向为正。 需要说明的是，实验中，我们考虑到磁场太小时即使线圈电流很大也不能使得天平有很大偏转。所以我们先在励磁电流为3A的情况下测出不同线圈电流所受的安培力的大小；然后再在线圈电流为3A的情况下测出不同励磁电流（及磁场强度）情况下的安培力大小。这里其实也使用到了控制变量法的基本原理。

用各组数据的前一半数据可画出如下四张图形： 从四张图中可以清楚看到，当磁场、线圈尺寸、线圈匝数不变时，安培力和线圈电流成正比例关系，即：F∝I。 为了方便分析，可以把数据都放到一个图中（如下）： 我们先分析一下A,B,C三组。这三组数据的共同点是线圈匝数相同，所以说，在磁场中线圈长度越长，安培力随电流的变化速率越快。那么两者的定量关系又如何呢？我们可以看到，线圈长度关系为：L(C)=2L(B)=4L(A)。而三者对应的安培力是什么关系呢？我们可以从下图中看出，变化率也是两倍关系。即：F∝L

数学实验作业

练习2﹒1 画出下列常见曲线的图形（其中a=1，b=2,c=3）。 1. 立方抛物线y = 解： x=-4:0.1:4; y=x.^(1/3); plot(x,y) -4 -3-2-101234 0.20.40.60.811.21.4 1.6 2.高斯曲线2 x y e -= 解： fplot('exp(-x^2)',[-4,4])

-4 -3 -2 -1 1 2 3 4 00.10.20.30.40.50.60.70.80.9 1 3、笛卡儿曲线23 3 2 2 33,(3)11at at x y x y axy t t = = +=++ 解：ezplot('x^3+y^3-3*x*y',[-4,4])

-4 -3-2-1 01234 -4-3-2-10123 4x y x 3+y 3-3 x y = 0 或：t=-4:0.1:4; x=3*t./(1+t.^2); y=3*t.^2./(1+t.^2); plot(x,y)

-1.5 -1-0.500.51 1.5 00.5 1 1.5 2 2.5 3 4、蔓叶线233 2 2 2 ,()11at at x x y y t t a x = = = ++- 解：t=-4:0.1:4; x=t.^2./(1+t.^2); y=t.^3,/(1+t.^2); y=t.^3./(1+t.^2); plot(x,y)

00.10.20.30.40.50.60.70.80.91 -4 -3-2-10123 4 或： ezplot('y .^2-x.^3/(1-x)',[-4,4])

北京科技大学参数检测实验报告全

北京科技大学参数检测实验报告全

实验六工业热电偶的校验 摘要：本实验重在了解热电偶的工作原理并通过对热电偶进行校正验证镍铬热电偶的准确性并了解补偿导线的使用方法。 关键词：热电偶校正标准被校补偿导线 1 引言 （1）实验目的 1．了解热电偶的工作原理、构造及使用方法。了解热电势与热端温度的关系。了解对热电偶进行校正的原因及校正方法，能独立地进行校正实验和绘制校正曲线。 2．了解冷端温度对测量的影响及补偿导线的使用方法。 3．通过测量热电势掌握携带式直流电位差计的使用方法。 （2）实验设备 1．铂铑－铂热电偶（标准热电偶）１支 ２．镍铬－镍硅热电偶（被校正热电偶）１支 3．热电偶卧式检定炉（附温度控制器）１台 4．携带式直流电位差计 1台 5．酒精温度计 1支 6．广口保温瓶 1个 7．热浴杯及酒精灯各1个 2 内容 1．了解直流电位差计各旋钮、开关及检流计的作用，掌握直流电位差计的使用方法。 2．热电偶校正 （１）实验开始，给检定炉供电，炉温给定值为400oC。当炉温稳定后，用电位差计分别测量标准热电偶和被校正热电偶的热电势，每个校正点的测量不得少于四次。数据记录于表6－１。 （２）依次校正600oC、 800oC、 1000oC各点。 （３）将测量电势求取平均值并转换成温度，计算误差，根据表6－３判断被热电偶是否合格。绘制校验曲线。 3．热电偶冷端温度对测温的影响及补偿导线的使用方法。 （１）1000oC校正点作完后，保持炉温不变。测量热浴杯中的水温，然后用电位差计分别测量镍铬－镍硅热电偶未加补偿导线和加补偿导线的热电势。数据记录于表6－２中。 （２）用酒精灯加热热浴杯，当水温依次为30oC、 40oC、 50oC时，用电位差计分别测量镍铬－镍硅热电偶未加补偿导线和加补偿导线的热电势。数据记录于表6－２中。 （３）用铂铑－铂热电偶测量炉温，检查实验过程中炉温是否稳定，分析若炉

MATLAB实验练习题(计算机) 南邮 MATLAB 数学实验大作业答案

“MATLAB”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> solve('exp(x)-3*x^2',0) ans = -2*lambertw(-1/6*3^(1/2)) -2*lambertw(-1,-1/6*3^(1/2)) -2*lambertw(1/6*3^(1/2)) 2、求下列方程的根。 1) 5510x x ++= a=solve('x^5+5*x+1',0);a=vpa(a,6)

1.10447+1.05983*i -1.00450+1.06095*i -.199936 -1.00450-1.06095*i 1.10447-1.05983*i 2） 1 sin0 2 x x-=至少三个根 >> fzero('x*sin(x)-1/2', 3) ans = 2.9726 >> fzero('x*sin(x)-1/2',-3) ans = -2.9726 >> fzero('x*sin(x)-1/2',0) ans = -0.7408

3）2sin cos 0x x x -= 所有根 >> fzero('sin(x)*cos(x)-x^2',0) ans = >> fzero('sin(x)*cos(x)-x^2',0.6) ans = 0.7022 3、求解下列各题： 1）30sin lim x x x x ->- >> sym x; >> limit((x-sin(x))/x^3) ans = 1/6 2） (10)cos ,x y e x y =求 >> sym x; >> diff(exp(x)*cos(x),10) ans =

北科大物理实验低真空的获得、测量与用直流溅射法制备金属薄膜

低真空的获得，测量与用直流实验报告 溅射法制备金属薄膜 学院材料班级材料1510 学号 41503000 姓名张问 一、实验目的与实验仪器实验目的： 1)学习真空基本知识和真空的获得与测量技术基础知识。 2)学习用直流溅射法制备薄膜的原理与方法。 3)实际操作一套真空镀膜装置，使用真空泵和真空测量装置，研究该真空系统的抽气特性。 4)用直流溅射法制备一系列不同厚度的金属薄膜，为实验研究金属薄膜厚度对其电阻率影响制备样品。实验仪器： SBC-12 小型直流溅射仪（配有银靶），机械泵，氩气瓶、超声波清洗器、玻璃衬底二、实验原理 （要求与提示：限400 字以内，实验原理图须用手绘后贴图的方式） 1.真空基本知识：该实验首先需制造一个真空条件，在中真空状态，机械泵抽 真空 2.用直流溅射法制备薄膜 溅射就是指荷能粒子轰击表面，使得固体表面原子从表面射出的现象，这些从固体表面 射出的粒子大多呈原子状态，通常称为被溅射原子。溅射粒子轰击靶材，使得表面成为被溅射原子，被溅射原子沉积到衬底上就形成了薄膜。所以这种方法称为溅射法，干燥气体在正常状态下是不导电的，若在气体中安置两个电极并加上电压，少量初始带电粒子与气体院子相互碰撞，使束缚电子脱离气体原子成为电子，这时有电流通过气体，这个现象称为气体放电。 本实验所使用的是直流溅射法，利用了直流电压产生的辉光放电，如图所示，在对系统抽真空后，充入适当压力的惰性气体，作为气体放电的载体，在正负极高压下，气体分子被大量电离，并伴随发出辉光。由放电形成的气体正离子被朝着阴极方向加速，正离子和快速中性粒子获得能量到达靶材，在这些粒子的轰击下，被溅射出来的靶材原子冷凝在阳极上，从而形成了薄膜。

数学实验报告-6

《数学实验》报告 实验名称常微分方程的求解 学院材料科学与工程 专业班级材料1209 姓名曾雪淇 学号 41230265 2014年 5月

一、【实验目的】 掌握常微分方程求解和曲线拟合的方法，通过MATLAB求解一阶甚至是二阶以上的高阶微分方程。 二、【实验任务】 P168习题24，习题27 三、【实验程序】 习题24：dsolve('Dy=x*sin(x)/cos(y)','x') 习题27：function xdot=exf(t,x) u=1-2*t; xdot=[0,1;1,-t]*x+[0 1]'*u; clf; t0=0; tf=pi; x0t=[0.1;0.2]; [t,x]=ode23('exf',[t0,tf],x0t) y=x(:,1); Dy=x(:,2); plot(t,y,'-',t,Dy,'o') 四、【实验结果】 习题24：ans = -asin(-sin(x)+x*cos(x)-C1) 习题27： t = 0.014545454545455 0.087272727272727 0.201440113885487 0.325875614772746 2

0.462108154525786 0.612058884594697 0.777820950596408 0.962141414226468 1.148168188604642 1.276725612086219 1.405283035567796 1.518837016595503 1.670603286779598 1.860122410374634 2.089084425249819 2.356884067351406 2.654570124097287 2.968729389456267 3.141592653589793 x = 0.100000000000000 0.200000000000000 0.103024424647132 0.215787876799993 0.121418223032493 0.288273863806750 0.159807571438023 0.379808018692957 0.211637169341158 0.447918********* 0.275587792496926 0.484712850141869 0.348540604264411 0.481263088285519 3

北京科技大学数学实验第五次讲解学习

北京科技大学数学实 验第五次

精品资料 《数学实验》报告 实验名称 Matlab拟合与插值 2013年12月

一、【实验目的】 1.学习Matlab的一些基础知识，主要多项式及其相关计算等； 2.熟悉Matlab中多项式的拟合，编写一些相关的Matlab命令等； 3.熟悉Matlab中多项式的插值，并编写一些相关的Matlab命令等； 4.完成相关的练习题。 二、【实验任务】 1.在钢线碳含量对于电阻的效应的研究中，得到以下数据.分别用一次、三次、五次多项式拟合曲线来拟合这组数据并画出图形，计算当x=0.45时的电阻值. 碳含量 0.10 0.30 0.40 0.55 0.70 0.80 0.95 x 电阻y 15 18 19 21 22.6 23.8 26 2.在某种添加剂的不同浓度之下对铝合金进行抗拉强度试验，得到数据如下，现分别使用不同的插值方法，对其中没有测量的浓度进行推测，并估算出浓度X=18及26时的抗压强度Y的值. 浓度X 10 15 20 25 30 抗压强度Y 25.2 29.8 31.2 31.7 29.4 3.用不同方法对在(-3,3)上的二维插值效果进行比较.

三、【实验程序】 1.在钢线碳含量对于电阻的效应的研究中，得到以下数据.分别用一次、三次、五次多项式拟合曲线来拟合这组数据并画出图形，计算当x=0.45时的电阻值. M文件 clc; clf; x=[0.1 0.3 0.4 0.55 0.7 0.8 0.95]; y=[15 18 19 21 22.6 23.8 26]; p1=polyfit(x,y,1); p3=polyfit(x,y,3); p5=polyfit(x,y,5); x1=0.1:0.05:1; y1=polyval(p1,x1); y3=polyval(p3,x1); y5=polyval(p5,x1); plot(x,y,'rp',x1,y1,'b-',x1,y3,'g-.',x1,y5,'m--'); legend('拟合点','一次拟合','三次拟合','五次拟合'); disp('以下为当x=0.45时的电阻值：') disp('一阶拟合函数值'),g1=polyval(p1,0.45) disp('三阶拟合函数值'),g3=polyval(p3,0.45) disp('五阶拟合函数值'),g5=polyval(p5,0.45)

matlab与数学实验大作业

《数学实验与MATLAB》 ——综合实验报告 实验名称：不同温度下PDLC薄膜的通透性 与驱动电压的具体关系式的研究学院：计算机与通信工程学院 专业班级： 姓名： 学号： 同组同学： 2014年 6月10日

一、问题引入 聚合物分散液晶(PDLC)是将低分子液晶与预聚物Kuer UV65胶相混合，在一定条件下经聚合反应，形成微米级的液晶微滴均匀地分散在高分子网络中，再利用液晶分子的介电各向异性获得具有电光响应特性的材料，它主要工作在散射态和透明态之间并具有一定的灰度。聚合物分散液晶膜是将液晶和聚合物结合得到的一种综合性能优异的膜材料。该膜材料能够通过驱动电压来控制其通透性，可以用来制作PDLC型液晶显示器等，具有较大的应用范围。已知PDLC薄膜在相同光强度及驱动电压下，不用的温度对应于不同的通透性，不同温度下的阀值电压也不相同。为了尽量得到不同通透性的PDLC薄膜，有必要进行温度对PDLC薄膜的特性的影响的研究。现有不同温度下PDLC 薄膜透过率与驱动电压的一系列数据，试得出不同温度下PDLC薄膜通透性与驱动电压的具体关系式，使得可以迅速得出在不同温度下一定通透性对应的驱动电压。 二、问题分析 想要得到不同温度下PDLC薄膜通透性与驱动电压的具体关系式可以运用MATLAB多项式农合找出最佳函数式，而运用MATLAB多项式插值可以得出在不同温度下一定通透性所对应的驱动电压。 三、实验数据 选择10、20、30摄氏度三个不同温度，其他条件一致。

（1）、10摄氏度 实验程序： x=2:2:40; y=[5.2,5.4,5.8,6.4,7.2,8.2,9.4,10.8,12.2,14.0,16.6,22.0, 30.4,39.8,51.3,55.0,57.5,58.8,59.6,60.2]; p3=polyfit(x,y,3); p5=polyfit(x,y,5); p7=polyfit(x,y,7); disp('三次拟合函数'),f3=poly2str(p3,'x') disp('五次拟合函数'),f5=poly2str(p5,'x') disp('七次拟合函数'),f7=poly2str(p7,'x') x1=0:1:40; y3=polyval(p3,x1); y5=polyval(p5,x1); y7=polyval(p7,x1); plot(x,y,'rp',x1,y3,'--',x1,y5,'k-.',x1,y7); legend('拟合点','三次拟合','五次拟合','七次拟合') 实验结果：

数学实验作业 韩明版

练习6.7 1.有两个煤厂A,B，每月进煤不少于60t,100t，它们担负供应三个居 民区的用煤任务，这三个居民区每月用煤量分别为45t,75t和45t.A 厂离这三个居民区的距离分别为10km,5km,6km,B厂离这三个居民区的距离分别为4km,8km,15km.问这两个煤厂如何分配供煤量能使总运输量（t.km）最小。 解：设甲对三个居民区的供煤量分别为：x1,x2,x3,乙对三个居民区的供煤量分别为x4,x5,x6.由已知有： y=10x1+5x2+6x3+4x4+8x5+15x6 -x1-x2-x3<=-60, -x4-x5-x6<=-100, x1+x4=45,x2+x5=75,x3+x6=40, X1>=0,x2>=0,x3>=0,x4>=0,x5>=0,x6>=0. 输入命令： > c=[10 5 6 4 8 15];A=[-1 -1 -1 0 0 0;0 0 0 -1 -1 -1;0 0 0 0 0 0;0 0 0 0 0 0;0 0 0 0 0 0;0 0 0 0 0 0]; >> b=[-60;-100;0;0;0;0];Aeq=[1 0 0 1 0 0;0 1 0 0 1 0;0 0 1 0 0 1;0 0 0 0 0 0;0 0 0 0 0 0;0 0 0 0 0 0]; >> beq=[45 75 40 0 0 0]; >> lb=ones(6,1); >> [x,fval]=linprog(c,A,b,Aeq,beq,lb) Optimization terminated.

结果为： x = 1.0000 20.0000 39.0000 44.0000 55.0000 1.0000 fval =975.0000 这说明甲乙两个煤厂分别对三个居民区输送1t 20t 39t,44t 55t 1t的煤才能使总运输量最小，且总运输量为975t.km 2．某银行经理计划用一笔资金进行有价证券的投资，可供购进的证券及其信用等级、到期年限、税前收益如下表所示。按照规定，市政证券的收益可以免税，其他证券的收益需按40%的税率纳税。此外还有以下限制： （1）政府及待办机构的证券总共至少购进400万元； （2）所构证券的平均信用等级不超过1.4（信用等级数字越小，信用程度越高）； （3）所构证券的平均到期年限不超过5年。

数学实验报告-2

《数学实验》报告 实验名称 MATLAB绘图 学院材料科学与工程 专业班级材料1209 姓名曾雪淇 学号 41230265 2014年 5月

学会用MATLAB绘制二维曲线、三维曲线，掌握gtext, legend, title,xlabel,ylabel,zlabel,axis 等指令用法,并学会图形的标注。二、【实验任务】 P79 习题1，习题3，习题5 三、【实验程序】 习题一： x=0:pi/10:4*pi; y1=exp(x./3).*sin(3*x); y2=exp(x./3); y3=-exp(x./3); plot(x,y1,'b*',x,y2,'r-.',x,y3,'r-.') 习题二： x1=-pi:pi/10:pi; y1=x1.*cos(x1); x2=pi:pi/10:4*pi; y2=x2.*tan(1./x2).*sin(x2).^3; x3=1:0.1:8; y3=exp(1./x3).*sin(x3); subplot(1,3,1);plot(x1,y1,'r*')，grid on,title(‘y1= x1*cosx1’) subplot(1,3,2) ;plot(x2,y2,’b-‘),grid on,title (‘y2=x2*tan(1/x2)*sinx2^3’) subplot(1,3,3);plot(x3,y3,'g+'),grid on,title (‘y3=exp(1/x3)*sinx3’) gtext(‘y1=x1cos(x1)’),gtext(‘y2=x2tan(1/x2)sin(x2)^3’), gtext(‘y3=exp(1/x3)sin(x3)’) legend(‘y1= x1*cos(x1)’, ‘y2=x2tan(1/x2)sin(x2^)3’ ‘y3=exp(1/x3)sin(x3)’) xlabel(‘x轴’),ylabel(‘y轴’),axis xy 习题三： t=0:pi/10:20*pi; x=t.*cos(pi/6.*t); y=t.*sin(pi/6.*t); z=2*t; plot3(x,y,z,'r*'),grid on title(‘圆锥螺线的图像’) xlabel(‘x轴’),ylabel(‘y轴’),zlabel(‘z轴’)

李萨如图模拟(Matlab大作业)

《数学实验》报告 实验名称李萨如图模拟（Matlab大作业） 2011年11月8日

一、【实验目的】 运用数学知识与MATLAB相结合，运用数学方法，建立数学模型，用MATLAB软件辅助求解模型，解决实际问题。 二、【实验任务】 一个质点沿 X轴和 Y轴的分运动都是简谐运动,分运动的表达式分别为: x=Acos ( w1t+beta ) , y=Acos(w2t+beta ) 。如果二者的频率有简单的整数比, 则相互垂直的简谐运动合成的运动将具有封闭的稳定的运动轨迹, 这种图称为李萨如图。 1，用matlab分别画出同一方向的传播波频率之比为2,3,4/5,1/2,1/3,5/4的图像（未合成）2，用matlab画出同一方向的传播波频率之比为2,3,4/5,1/2,1/3,5/4的合成图像 3，用matlab画出x轴方向和y轴方向传播波频率之比为2,3,4/5,1/2,1/3,5/4的合成图像。（李萨如图） 三、【实验分析及求解】 1，设两个波的振幅为1，他们的beta为pi/5,我们可以根据波的传播公式，y =Acos ( w1t+beta ) 分别画出两个波的传播图像。 2，设两个波的振幅为1，他们的beta为pi/5,我们可以根据波的传播公式，y =Acos ( w1t+beta )， 用matlab画出同一方向的传播波频率之比为2,3,4/5,1/2,1/3,5/4的合成图像。

3，设两个波的振幅为1，他们的beta为pi/5,我们可以根据波的传播公式，画出x轴方向和y 轴方向传播波频率之比为2,3,4/5,1/2,1/3,5/4的合成图像。（李萨如图）。

北京科技大学工科物理实验下

全息照相5.5 1、画出全息照相光路图 P140图5.5-3 2、全息照相的影响因素？ a、系统的稳定性及实验环境是否符合要求 b、参考光与物光的光程差 c、参考光与物光的夹角（30-50）为宜 d、参考光与物光的光强比（3:1-5:1）为宜 e、底片的安放 f、曝光与显影时间的控制（显影1-2分颜色变，定影2-4分，水洗2-3分） 3、全息照相和普通照相的区别 P138 普通照相是用几何光写的方法记录物体上各点的发光强度分布，得到的是二维平面像，像上各点的照度与物体上各点的发光强度一一对应。而全息照相的记录对象是整个物体发出的光波，借助于参考光用干涉的方法记录这个物光波的振幅和位相分布，即记录下物光波与参考光波相干后的全部信息。 4、全息照相与普通照相的本质区别？ 答：全息照相和常规照相之不同在于，常规照相只是记录了被摄物体表面光线强弱的变化，即只记录了光的振幅；而全息照相则记录了光波的全部信息，除振幅外，还记录了光波的相位。这样就把空间物体光波场的全部信息都贮存记录了下来。然后利用全息照片对特定波长单色照明光的衍射，把原空间

景象显现出来。 5、全息照相实验失败率很高，怎样改进？ 调节物光束与参考光束由分光板至感光板的光程大致相等，用透镜将物光束扩展到整个被摄物都受到均匀的光照，物光束与参考光束的夹角应为30到50度，调整参考光束与物光束的光强比，关闭室内照明灯减少其他光的干扰，控制曝光时间。 6、做全息照相时，能否用一张底片记录两个物体的像？ 7、全息照相中物象光和参考光的夹角小于多少度？ 照在感光板上的物光束和参考光束的夹角在30°~40°为宜 8、全息照相的特点？ 8、全息照相原理？ 9、全息照相利用光的干涉，将物体的（振幅）和（相位）特性以（干涉条纹）形式表现出来。 10、影响全息照相成败的条件？ 11、扩光束的作用是什么?分光束分得的光，光强较大的是物光还是参考光，为什么？ 12、全息扩光镜的作用（两个分别写） 13、全息照相中的感光片是否有正反面？ 有正反面，有感光乳胶的面为正面。 传感器7.7 1、传感器的定义 传感器是一种能把特定的被测量信息（包括物理量、化学量、生物学量等）按一定规律转换成某种可用信号输出的器件或装置。 2、传感器通常由哪两部分组成？各部分的作用是什么？

数学实验第七次作业

4. 问题： 某公司将3种不同含硫量的液体原料（分别记为甲、乙、丙）混合生产两种产品（分别记为A,B ）。按照生产工艺的要求，原料甲、乙必须首先导入混合池中混合，混合后的液体再分别与原料丙混合生产A,B 。一直原料甲、乙、丙的含硫量分别是3%，1%，2%，进货价格分别为6千元/t ，16千元/t ，10千元/t ；产品A,B 的含硫量分别不能超过2.5%，1.5%，售价分别为9千元/t ，15千元/t 。根据市场信息，原料甲、乙、丙的供应量都不能超过500t ；产品A,B 的最大市场需求量分别为100t ，200t 。 (1) 应如何安排生产？ (2) 如果产品A 的最大市场需求量增长为600t ，应如何安排生产？ (3) 如果乙的进货价格下降为13千元/t ，应如何安排生产？分别对(1)、(2)两种情况进 行讨论。 模型： （只考虑问题1，问题2,3只需改变一些约束条件） 设生产时使用原料甲、乙分别为12,x x t ，分别取混合后的液体34,x x t 再加入原料丙 56,x x t 生产产品A,B 。 有质量守恒，可得 1234x x x x +=+ 甲乙混合后的液体的含硫量可表示为 12 12 3%x x x x ++，根据含硫量的要求，可得 12 353512 124646 12 3%*2%* 2.5%*()3%*2%* 1.5%*() x x x x x x x x x x x x x x x x +?+≤+?+?? +?+≤+?+? 根据市场的限制，易得 12563546500 500500100200 x x x x x x x x ≤?? ≤?? +≤??+≤??+≤? 当然还有非负约束 123456,,,,,0x x x x x x ≥ 公司的净利润为（单位：千元）： 35461256123456 9()15()61610()6169155z x x x x x x x x x x x x x x =+++---+=--++-+

MTLB实验练习题计算机南邮MATLAB数学实验大作业答案

“M A T L A B ”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> solve('exp(x)-3*x^2',0) ans = -2*lambertw(-1/6*3^(1/2)) -2*lambertw(-1,-1/6*3^(1/2)) -2*lambertw(1/6*3^(1/2)) 2、求下列方程的根。 1) 5510x x ++= a=solve('x^5+5*x+1',0);a=vpa(a,6) a = 1.10447+1.05983*i -1.00450+1.06095*i -. -1.00450-1.06095*i

1.10447-1.05983*i 2） 1 sin0 2 x x-=至少三个根 >> fzero('x*sin(x)-1/2', 3) ans = 2.9726 >> fzero('x*sin(x)-1/2',-3) ans = -2.9726 >> fzero('x*sin(x)-1/2',0) ans = -0.7408 3）2 sin cos0 x x x -=所有根 >> fzero('sin(x)*cos(x)-x^2',0) ans = >> fzero('sin(x)*cos(x)-x^2',0.6)

0.7022 3、求解下列各题： 1）3 0sin lim x x x x ->- >> sym x; >> limit((x-sin(x))/x^3) ans = 1/6 2） (10)cos ,x y e x y =求 >> sym x; >> diff(exp(x)*cos(x),10) ans = (-32)*exp(x)*sin(x) 3）2 1/2 0(17x e dx ?精确到位有效数字） >> sym x; >> vpa((int(exp(x^2),x,0,1/2)),17)

北京科技大学_工科物理实验II题集

操作题 1、受迫振动测振幅 2、(据说，操作题不考全系照相） 3、(据说，K组做牛顿环实验) 4、弗朗克赫兹实验（A组的），...画两个峰值的就行了 5、声速的测量——XX和三个不同组的同学不约而同的考到了这个实验！建议大家好好把公式看两遍，能背就背吧，尤其是那些恶心人的误差公式，真的是要一步一步算的！！！实验特好做，数据特纠结！ 6、7.4实验（非牛顿环部分），听说和咱们实验时做的不一样，就是让测量是个条纹的间距，都是测距离，还算好做。 7、分光仪实验（我们没有做，是别的同学做的） 笔试题： 1、金属镀膜中电压值是多大？在镀膜中的作用？ V=1000V，在接近真空的条件下，高电压使氩气进行“辉光放电”，氩原子电离成氩离子，并在电场力的作用下，加速轰击靶材，靶材原子则会被溅射出来而沉积到玻璃衬底上。 2、本学期测金属薄膜的实验中的方法能否测所有厚度的薄膜？说明理由。 （不行！（这个肯定对吧）。理由：51#：不能。膜太厚的话，将不能分辨出左右（有膜和无膜处）到底相差几个二分之一波长。厚膜需要换一种方法：将空气劈尖转90° 就能观测了。如图所示 3、传感器的定义 4、画出全息照相光路图 5、全息照相的影响因素 6、什么是禁带宽度 禁带宽度是指一个能带宽度.固体中电子的能量是不可以连续取值的，而是一些不连续的能带。要导电就要有自由电子存在。自由电子存在的能带称为导带。被束缚的电子要成为自由电子，就必须获得足够能量从而跃迁到导带，这个能量的最小值就是禁带宽度。禁带非常窄就成为金属了，反之则成为绝缘体。半导体的反向耐压，

正向压降都和禁带宽度有关。 7、锗的适用温度小于45℃ 8、还有个望远镜调法（答案是自准法） 9、游标卡尺上读两个示数的原因（棱镜分光仪吧，消除偏心差） 10、画出受迫振动的幅频相频曲线 11、关于误差分配的 12、最小偏向角与波长关系的 13、改变原子状态一般几种方式，分别是？ 1、以一定频率的光子来实验2、具有一定能量的原子与其碰撞进行能量交换 15、声速的理论值公式是什么？其中的T指的是什么？c0的值是多少？（此题变态程度一般，看了就会，没看就真是不知道。。。那个c0的值不是应该直接给出的么！！！！） 16、在单色仪实验中，如何提高分辨能力？ 17、全息照相和普通照相的区别（巨简单巨简单！！！） 18、6.6设计测量电阻率实验的电路图 19、7.7传感器实验中三个金属片哪个在测量时出现了零值区，为什么？ 铁片，因为金属片的磁导率和电导率越大，线圈和金属片的间距越小，涡流的反作用越强，输出电压越低，甚至输出为零，所以铁片有零值区出现。 20、全息照相中的感光片是否有正反面？ 21、当标准电阻R为0.1级时，其误差为多少？（在6.6实验中实验器材里有） 22、单色仪的光路图 23、涡流传感器的实验原理 由于电流的周期性变化，产生交变磁场H1，金属片靠近时产生磁场H2，电感线圈的等效阻抗发生变化，当各参数、系数保持不变时，阻抗Z只与线圈和金属片的距离 X有关。 24、四引线法适用于什么电阻 测量电阻阻值较低的电阻 25、全息照相实验失败率很高，怎样改进？ 27、四探针法与四引线法有什么异同？ 同：都是用来测量低电阻阻值的方法，都比较有效低消除了接线电阻和接触电阻的影响。异：四探针法比四引线法测得低电阻更为精确，可以测半导体，薄膜的电阻。 28、有和我们做的不一样实验的同学考到，声速的测量，全息照相，高温超导，幅频特性和相频特性曲线，四选二画图。 29、受迫振动实验中驱动力与***（忘记了，大家看书吧~我实在是懒得查了，书上有）为什么查90°相位角。 当驱动力f与震动物体固有f相同时，受迫振动速度幅max，产生速度共振，所以物体振动位移比驱动力之后90°。