实验1 一元函数作图与性质研究
实验目的
1.学习绘制一元函数的图形
2.观察函数的特性（单调性、有界性、奇偶性、周期性），建立数形结合的思想。实验准备
1.基本初等函数及基本性质
2.分段函数,如阶梯函数、单位阶跃函数等。
3.平面曲线的参数表示，如圆，椭圆，直线等。
4.极坐标曲线，如心形线，玫瑰线的等。
实验内容
1.函数图形的绘制
1)直角坐标系下函数的图形并标注
2)参数方程（含极坐标方程）表示的函数图形
3)隐函数的图形
4)分段函数的图形
2.绘制点列的散点图和折线图
3.观察函数的几何特性
MATLAB命令
表1-1 Matlab一元函数绘图命令
- 4 - 第一章 基础实验
4
实验示例
【例1.1】函数的图形
设y 的范围为[-10,10]，画出函数sin =y x 和sec =y x 在区间[2,2]ππ-上的图形，并标注图像。 【程序】：
clf;
hold on
axis([-2*pi 2*pi -10 10])
fplot('sin',[-2*pi 2*pi],'r-') fplot('sec',[-2*pi 2*pi],'b--') h=legend('sinx','secx',1); hold off 或者
clf; hold on
axis([-2*pi 2*pi -10 10]) x=-2*pi:0.1:2*pi; y1=sin(x); y2=sec(x);
plot(x,y1,'r-',x,y2,'b--'); h=legend('sinx','secx',1); hold off
【输出】：见图
1-1
-6-4-20246
-10
-5
5
10
图1-1 函数的图形
【例1.2】参数方程绘图
在同一坐标系中绘制3
3
3cos 3sin ?=??=??x t
y t
与2(sin )2(1cos )=-??=-?x t t y t 在区间[0,2]π上的图形。
实验一 一元函数绘图 5 -
5
【程序】：
clf;
t=0:0.01:2*pi;
x1=3*(cos(t)).^3;y1=3*(sin(t).^3); x2=2*(t-sin(t));y2=2*(1-cos(t));
plot(x1,y1,'r-',x2,y2,'b--','LineWidth',2)
legend('(3cos(t),3sin(t))','6sin(2t),6cos(2t)',1)
【输出】：见图1-2。
-5
051015
图1-2 参数方程绘图
【例1.3】极坐标绘图
画出心形线(1cos ),1,2,02r b at a b t π=-==≤≤的图形。 【程序】：
方法一：利用函数polar a=1;b=2; t=0:0.01:2*pi; r=b*(1-cos(a*t)); polar(t,r);
方法二：将极坐标化成参数方程
cos sin x r t
y r t
=??
=? a=1;b=2;
t=linespace(0,2*pi,1000); r=b*(1-cos(a*t); x=r.*cos(t); y=r.*sin(t); plot(x,y); axis equal; grid;
【输出】：见图1-3。
- 6 - 第一章 基础实验
6
图1-3 极坐标绘图
【例1.4】分段函数绘图
画出函数2sin ,0
(),22cos ,0
【程序】：
%函数定义开始
function y=DemoEx1_4(x) n=length(x);y=zeros(n,1); for i=1:n if x(i)<0
y(i)=sin(x(i)); else
y(i)=x(i)^2*cos(x(i)); end
end %函数定义结束
绘图代码： x=-2:0.01:2; x=DemoEx1_4(x); plot(x,y,’r-‘); grid ;
【输出】：见图1-4。
-2
-101
2
图1-4 分段函数绘图
【例1.5】散点图与折线图
画出表1-3中数据点所对应的散点图。
实验一 一元函数绘图 7 -
7
表1-3 例1.5的数据
【程序】：
x=[0 3 5 7 9 11 12 13 14 15];
y=[0 1.2 1.7 2.0 2.1 2.0 1.8 1.2 1.0 1.6]; plot(x,y,'*') 【输出】：见图1-5
图1-5 散点图
实验练习
1.画出sin cos y x x =+在区间[2,2]ππ-上的图形，设定y 的范围为[1,1]-。
2.画出参数方程1123
11
23cos cos 7sin 17,02sin sin 7cos17x t t t t y t t t
π=++?≤≤?=++?的图形。
3.绘制由隐函数4422
810160x y x y +--+=确定的曲线。 4.绘制函数cos 4
34
2cos sin
,02e
θ
θρθθπ=-+≤≤的图形。
5.绘制分段函数21sin
,0()0
,0
x
x x f x x ?≠?=?
=??的图形。
6.在同一坐标系中画出下列函数
123()sin ,()sin(sin ),()sin(sin(sin ))f x x f x x f x x ===
的图形，并观察它们的变化趋势，并猜测对于任意的x ,
sin
()sin(sin(())n n f x x =
个
有什么变化趋势。
北京交通大学海滨学院考试试题 课程名称：数学实验2010-2011第一学期出题教师：数学组适用专业: 09机械, 物流, 土木, 自动化 班级：学号：姓名： 选做题目序号： 1.一对刚出生的幼兔经过一个月可以长成成兔, 成兔再经过一个月后可以 繁殖出一对幼兔. 如果不计算兔子的死亡数, 请用Matlab程序给出在未来24个月中每个月的兔子对数。 解: 由题意每月的成兔与幼兔的数量如下表所示： 1 2 3 4 5 6 ··· 成兔0 1 1 2 3 5··· 幼兔 1 0 1 1 2 3··· 运用Matlab程序： x=zeros(1,24); x(1)=1;x(2)=1； for i=2:24 x(i+1)=x(i)+x(i-1); end x 结果为x = 1 1 2 3 5 8 13 21 3 4 5 5 89 144 233 377 610 987 1597 2584 4181 6765 1094 6 7711 2865 7 46368 2.定积分的过程可以分为分割、求和、取极限三部分, 以1 x e dx 为例, 利用
已学过的Matlab 命令, 通过作图演示计算积分的过程, 并与使用命令int() 直接积分的结果进行比较. 解：根据求积分的过程，我们先对区间[0,1]进行n 等分， 然后针对函数x e 取和，取和的形式为10 1 i n x i e e dx n ξ=≈ ∑ ? ，其中1[ ,]i i i n n ξ-?。这里取i ξ为区间的右端点，则当10n =时，1 x e dx ?可用10 101 1.805610 i i e ==∑ 来近似计算， 当10n =0时，100 100 1 01 =1.7269100 i x i e e dx =≈ ∑?，当10n =000时，10000 10000 1 1 =1.718410000 i x i e e dx =≈ ∑ ?. 示意图如下图，Matlab 命令如下： x=linspace (0,1,21); y=exp(x); y1=y(1:20); s1=sum(y1)/20 y2=y(2:21); s2=sum(y2)/20 plot(x,y); hold on for i=1:20 fill([x(i),x(i+1),x(i+1),x(i),x(i)],[0,0,y(i),y(i),0],'b') end syms k;symsum(exp(k/10)/10,k,1,10)；%n=10 symsum(exp(k/100)/100,k,1,100)；%n=100 symsum(exp(k/10000)/10000,k,1,10000)；%n=10000
安徽工业大学 大学数学实验课程设计 姓名： 班级： 任课老师：
数学实验 课程设计 问题提出： 某容器盛满水后，低端直径为0d 的小孔开启（图）。根据水力学知识，当水面 高度h 时，水冲小孔中流出的速度v =(g 为重力加速度，0.6为孔口的收缩系数)。 ⑴若容器为倒圆锥形（如图1），现测得容器高和上底面直径均为1.2m ，小孔直径为3cm ，问水从小孔中流完需要多长时间；2min 水面高度是多少。 ⑵若容器为倒葫芦形（如图2），现测得容器高为1.2m ，小孔直径为3cm ，有低端（记作x=0）向上每隔0.1m 测出容器的直径D （m ）如表所示，问水从小孔中流完需要多少时间；2min 时水面的高度是多少。 图1 ： 图2： 问题分析： （1） 倒圆锥形容器流水问题中随时间t 液面高度h 也在变化，同时水的流速也 在变化，再写变化难以用普通的方程进行模拟求解，考虑建立常微分方程竟而代入数值求解。水面的直径等于液面的高度。可以建立容器中水流失的液面高度对时间t 的变化率。 假设t 时，液面的高度h ，此时水的流速流量Q 为：00.6(/4)d π ； 则 在t ?时间内液面下降高度为h ?，可得到关系式：220( )2 4 d dt h dh π = ；
由此可知水下降h ? 时需要的时间：20 40.6 4 h dh t d π π ?= = 根据此关系式知道。 （2） 在第二问中，考虑倒葫芦形容器时因为他的高度h 不同容器直径D 变化 没有规律可循，同第一题相比我们只知道他的一些数值，这就需要我们建立高度h 和容器直径D 之间的关系矩阵，然后再欧拉方程和龙格—库塔方法找出时间t 和液面高度之间的分量关系。 由（1）可同理推知：假设在时间t 时，液面高度为h ，此时流量 为 2 00.6(/4)d π；经过t ?时，液面下降h ?，若我们取的t 是在t(n)和t(n+1) 之间的某一时刻，于是就可在误差范围内得到 (1)()t n t n t +=+?；可以得 到 204 (1)()0.64 h d h dt t n t n d π π =+-=- = ； 建立模型： （1） 在试验中我们不考虑圆锥的缺省对流水的影响，以及其他外界因素和玻璃 的毛细作用，试验中水可以顺利流完。实验中重力加速度g=9.82 /m s ；倒圆锥的液面最初高度为H=1.2m ，液面直径D=1.2m=0.03，小孔的直径为 0d =0.03m ； 接上文中分析结论代入数据：即在T 时间内将1.2m 的液面高度放完， （matlab 不支持一些运算符号，故用matlab 运算格式） dt=-((pi/4)h^2*dh)/(0.6*(pi/4)*d^2*sqrt(gh))=-(h^1.5*dh)/(0.6*d^2*sqrt(g)) h 是由0→1.2m 对t 积分 用matlab 计算上式 编辑文件：a1.m ， d0=0.03; g=9.8; syms h t=(h^1.5)/(0.6*d0^2*sqrt(g)); T=int(t,0,1.2); eval(T) 运行结果： >> a1 ans =
一牛顿环的各环是否等宽?密度是否均匀?解释原因? 因为环是由空气劈上下表面反射的两束光叠加干涉形成的。劈的上表面变化在横向是不均匀的，故光程差也不是均匀变化的。所以各环是不等宽的环的密度也不是均匀的。各环不等宽，半径小的环宽，越到外边越窄，密度是不均匀的,牛顿环的半径公式是：半径r等于根号下(m+1/2)λR，其中m为环的级数。从公式可以看出，半径和环数并不是线性关系，这样环自然不均匀。计算可以知道，越往外环越密。 二牛顿环的干涉圆环是由哪两束相干光干涉产生的? 半凸透镜下表面和下底面上表面的两束反射光 三电桥由哪几部分组成?电桥平衡的条件? 由电源、开关、检流计桥臂电阻组成。 平衡条件是Rx=(R1/R2)R3 四接通电源后,检流计指针始终向一边偏转,试分析出现这种情况的原因? 指针向一侧偏转就说明发生了电子的定向移动了，这个应该没问题。 指针不偏转，有2种情况吧，其1呢是整个电路发生了断路或其他故障，还1种情况则是流过的电流太小，不足于使电表发生偏转或其偏转的角度肉眼根本看不到。 无论如何调节，检流计指针都不动，电路中可能出现故障是调节臂电阻断路或短路。。无论如何调节，检流计指针始终像一边偏而无法平衡，电路中有可能出现故障是有一个臂（非调节臂）的电阻坏了。（断路或短路） 五什么叫铁磁材料的磁滞现象? 铁磁物质经外磁场磁化到饱和以后，把磁场去掉。这些物质仍保留有剩余磁化强度。需要反方向加磁场才能把这剩余磁化强度变为零。这种现象称为铁磁的磁滞现象。也是说，铁磁材料的磁状态，不仅要看它现在所处的磁场条件；而且还要看它过去的状态。 六如何判断铁磁材料属于软.硬材料? 软磁材料的特点是：磁导率大，矫顽力小，磁滞损耗小，磁滞回线呈长条状；硬磁材料的特点是：剩磁大，矫顽力也大 用光栅方程进行测量的条件是什么? 条件是一束平行光垂直射入光栅平面上，光波发生衍射，即可用光栅方程进行计算。如何实现：使用分光计，光线通过平行光管射入，当狭缝位于透镜的焦平面上时，就能使射在狭缝上的光经过透镜后成为平行光 用光栅方程进行测量,当狭缝太窄或者太宽会怎么样?为什么? 缝太窄，入射光的光强太弱，缝太宽，根据光的空间相干性可以知道，条纹的明暗对比度会下降！ 区别是，太窄了，亮纹会越来越暗，暗纹不变，直到一片黑暗！ 太宽，暗条纹会逐渐加强，明纹不变，直到一片光明！
一、实验内容 P206第六题 function f=wuyan2(c) y=[3.9 5.3 7.2 9.6 12.9 17.1 23.2 31.41 38.6 50.2 62.9 76.0 92.0 106.5 123.2 131.7 150.7 179.3 204.0 226.5 251.4 281.4] t=[0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200 210] f=y-c(1)/(1+c(1)/3.9-1)*exp^(-c(2)*t) c0=[1 1] c=lsqnonlin('wuyan2',c0) P206第七题 function f=wuyan1(c) q=[0.4518 0.4862 0.5295 0.5934 0.7171 0.8964 1.0202 1.1963 1.4928 1.6909 1.8548 2.1618 2.6638 3.4634 4.6759 5.8478 6.7885 7.4463 7.8345 8.2068 8.9468 9.7315 10.5172 11.7390 13.6876 ]; k=[0.0911 0.0961 0.1230 0.1430 0.1860 0.2543 0.3121 0.3792 0.4754 0.4410 0.4517 0.5595 0.8080 1.3072 1.7042 2.0019 2.2914 2.4941 2.8406 2.9855 3.2918 3.7214 4.3500 5.5567 7.0477]; l=[4.2361 4.3725 4.5295 4.6436 4.8179 4.9873 5.1282 5.2783 5.4334 5.5329 6.4749 6.5491 6.6152 6.6808 6.7455 6.8065 6.8950 6.9820 7.0637 7.1394 7.2085 7.3025 7.3470 7.4432 7.5200]; f=q-c(1)*k.^c(2).*l.^c(3) c0=[1 1 1] c=lsqnonlin('wuyan1',c0) c = 0.4091 0.6401 1.1446 a=0.4091 α=0.6401 β=1.1446 P239第五题 c=[-20 -30]; A=[1 2;5 4]; b=[20 70]; v1=[0 0]; [x,f,ef,out,lag]=linprog(c,A,b,[],[],v1) z=-f x = 10.0000 5.0000
第一章 函数、极限与连续 由于社会和科学发展的需要，到了17世纪，对物体运动的研究成为自然科学的中心问题.与之相适应，数学在经历了两千多年的发展之后进入了一个被称为“高等数学时期”的新时代，这一时代集中的特点是超越了希腊数学传统的观点，认识到“数”的研究比“形”更重要，以积极的态度开展对“无限”的研究，由常量数学发展为变量数学，微积分的创立更是这一时期最突出的成就之一.微积分研究的基本对象是定义在实数集上的函数. 极限是研究函数的一种基本方法，而连续性则是函数的一种重要属性.因此，本章内容是整个微积分学的基础.本章将简要地介绍高等数学的一些基本概念，其中重点介绍极限的概念、性质和运算性质，以及与极限概念密切相关的，并且在微积分运算中起重要作用的无穷小量的概念和性质.此外，还给出了两个极其重要的极限.随后，运用极限的概念引入函数的连续性概念，它是客观世界中广泛存在的连续变化这一现象的数学描述. 第一节 变量与函数 一、变量及其变化范围的常用表示法 在自然现象或工程技术中，常常会遇到各种各样的量.有一种量，在考察过程中是不断变化的，可以取得各种不同的数值，我们把这一类量叫做变量；另一类量在考察过程中保持不变，它取同样的数值，我们把这一类量叫做常量.变量的变化有跳跃性的，如自然数由小到大变化、数列的变化等，而更多的则是在某个范围内变化，即该变量的取值可以是某个范围内的任何一个数.变量取值范围常用区间来表示.满足不等式a x b ≤≤的实数的全体组成的集合叫做闭区间，记为,a b ????，即 ,{|}a b x a x b =≤≤????； 满足不等式a x b <<的实数的全体组成的集合叫做开区间，记为(,)a b ，即 (,){|}a b x a x b =<<； 满足不等式a x b <≤(或a x b ≤<)的实数的全体组成的集合叫做左(右)开右(左)闭区间，记为 (,a b ?? (或),a b ??)，即 (,{|}a b x a x b =<≤?? (或),{|}a b x a x b =≤
实验一霍尔效应及其应用 【预习思考题】 1．列出计算霍尔系数、载流子浓度n、电导率σ及迁移率μ的计算公式，并注明单位。 霍尔系数，载流子浓度，电导率，迁移率。 2．如已知霍尔样品的工作电流及磁感应强度B的方向，如何判断样品的导电类型？ 以根据右手螺旋定则，从工作电流旋到磁感应强度B确定的方向为正向，若测得的霍尔电压为正，则样品为P型，反之则为N型。 3．本实验为什么要用3个换向开关？ 为了在测量时消除一些霍尔效应的副效应的影响，需要在测量时改变工作电 流及磁感应强度B的方向，因此就需要2个换向开关；除了测量霍尔电压，还要测量A、C间的电位差，这是两个不同的测量位置，又需要1个换向开关。总之，一共需要3个换向开关。 【分析讨论题】 1．若磁感应强度B和霍尔器件平面不完全正交，按式（5.2-5）测出的霍尔系数比实际值大还是小？要准确测定值应怎样进行？ 若磁感应强度B和霍尔器件平面不完全正交，则测出的霍尔系数比实际值偏小。要想准确测定，就需要保证磁感应强度B和霍尔器件平面完全正交，或者设法测量出磁感应强度B和霍尔器件平面的夹角。 2．若已知霍尔器件的性能参数，采用霍尔效应法测量一个未知磁场时，测量误差有哪些来源？ 误差来源有：测量工作电流的电流表的测量误差，测量霍尔器件厚度d的长度测量仪器的测量误差，测量霍尔电压的电压表的测量误差，磁场方向与霍尔器件平面的夹角影响等。 实验二声速的测量 【预习思考题】 1. 如何调节和判断测量系统是否处于共振状态？为什么要在系统处于共振的条件下进行声速测定？ 答：缓慢调节声速测试仪信号源面板上的“信号频率”旋钮，使交流毫伏表指针指示达到最大（或晶体管电压表的示值达到最大），此时系统处于共振状态，显示共振发生的信号指示灯亮，信号源面板上频率显示窗口显示共振频率。在进行声速测定时需要测定驻波波节的位置，当发射换能器S1处于共振状态时，发射的超声波能量最大。若在这样一个最佳状态移动S1至每一个波节处，媒质压缩形变最大，则产生的声压最大，接收换能器S2接收到的声压为最大，转变成电信号，晶体管电压表会显示出最大值。由数显表头读出每一个电压最大值时的位置，即对应的波节位置。因此在系统处于共振的条件下进行声速测定，可以容易和准确地测定波节的位置，提高测量的准确度。 2. 压电陶瓷超声换能器是怎样实现机械信号和电信号之间的相互转换的？ 答：压电陶瓷超声换能器的重要组成部分是压电陶瓷环。压电陶瓷环由多晶结构的压电材料制成。这种材料在受到机械应力，发生机械形变时，会发生极化，同时在极化方向产生电场，这种特性称为压电效应。反之，如果在压电材料上加交
实验2 方程模型及其求解算法 一、实验目的及意义 [1] 复习求解方程及方程组的基本原理和方法； [2] 掌握迭代算法； [3] 熟悉MATLAB软件编程环境；掌握MATLAB编程语句(特别是循环、条件、控制等语句)； [4] 通过范例展现求解实际问题的初步建模过程； 通过该实验的学习，复习和归纳方程求解或方程组求解的各种数值解法（简单迭代法、二分法、牛顿法、割线法等），初步了解数学建模过程。这对于学生深入理解数学概念，掌握数学的思维方法，熟悉处理大量的工程计算问题的方法具有十分重要的意义。 二、实验内容 1．方程求解和方程组的各种数值解法练习 2．直接使用MATLAB命令对方程和方程组进行求解练习 3．针对实际问题，试建立数学模型，并求解。 三、实验步骤 1．开启软件平台——MATLAB，开启MATLAB编辑窗口； 2．根据各种数值解法步骤编写M文件 3．保存文件并运行； 4．观察运行结果(数值或图形)； 5．根据观察到的结果写出实验报告，并浅谈学习心得体会。 四、实验要求与任务 基础实验 1．用图形放大法求解方程x sin(x) = 1. 并观察该方程有多少个根。 画出图形程序： x=-10:0.01:10; y=x.*sin(x)-1; y1=zeros(size(x)); plot(x,y,x,y1) MATLAB运行结果：
-10-8-6-4-20246810 -8-6 -4 -2 2 4 6 8 扩大区间画图程序： x=-50:0.01:50; y=x.*sin(x)-1; y1=zeros(size(x)); plot(x,y,x,y1) MATLAB 运行结果： -50-40-30-20-1001020304050 由上图可知，该方程有偶数个无数的根。
中北大学 信息与通信工程学院 通信工程专业 《电子线路及系统》课程设计任务书2016/2017 学年第一学期 学生姓名：张涛学号： 李子鹏学号： 课程设计题目：数字电子钟的设计 起迄日期：2017年1月4日～2017年7月10日 课程设计地点：科学楼 指导教师：姚爱琴 2017年月日 课程设计任务书
中北大学 信息与通信工程学院 通信工程专业 《电子线路及系统》课程设计开题报告2016/2017 学年第一学期 题目：数字电子钟的设计 学生姓名：张涛学号： 李子鹏学号：
指导教师：姚爱琴 2017 年 1 月 6 日 中北大学 信息与通信工程学院 通信工程专业 《电子线路及系统》课程设计说明书2016/2017 学年第二学期 题目：数字电子钟的设计 学生姓名：张涛学号： 李子鹏学号： 指导教师：姚爱琴 2017 年月日
目录 1 引言 (6) 2 数字电子钟设计方案 (6) 2.1 数字计时器的设计思想 (6) 2.2数字电路设计及元器件参数选择 (6) 2.2.2 时、分、秒计数器 (7) 2.2.3 计数显示电路 (8) 2.2.5 整点报时电路 (10) 2.2.6 总体电路 (10) 2.3 安装与调试 (11) 2.3.1 数字电子钟PCB图 (11) 3 设计单元原理说明 (11) 3.1 555定时器原理 (12) 3.2 计数器原理 (12) 3.3 译码和数码显示电路原理 (12) 3.4 校时电路原理 (12) 4 心得与体会 (12) 1 引言 数字钟是一种用数字电子技术实现时，分，秒计时的装置，具有较高的准确性和直观性等各方面的优势，而得到广泛的应用。此次设计数字电子钟是为了了解数字钟的原理，在设计数字电子钟的过程中，用数字电子技术的理论和制作实践相结合，进一步加深数字电子技术课程知识的理解和应用，同时学会使用Multisim电子设计软件。 2数字电子钟设计方案 2.1 数字计时器的设计思想 要想构成数字钟，首先应选择一个脉冲源——能自动地产生稳定的标准时间脉冲信号。而脉冲源产生的脉冲信号地频率较高，因此，需要进行分频，使得高频脉冲信号变成适合于计时的低频脉冲信号，即“秒脉冲信号”（频率为1Hz）。经过分频器输出的秒脉冲信号到计数器中进行计数。由于计时的规律是：60秒=1分，60分=1小时，24小时=1天，就需要分别设计60进制，24进制计数器，并发出驱动信号。各计数器输出信号经译码器、驱动器到数字显示器，是“时”、“分”、“秒”得以数字显示出来。 值得注意的是：任何记时装置都有误差，因此应考虑校准时间电路。校时电路一般
东华大学M A T L A B数学实验第二版答案(胡良 剑) -CAL-FENGHAI-(2020YEAR-YICAI)_JINGBIAN
数学实验答案 Chapter 1 Page20,ex1 (5) 等于[exp(1),exp(2);exp(3),exp(4)] (7) 3=1*3, 8=2*4 (8) a为各列最小值，b为最小值所在的行号 (10) 1>=4,false, 2>=3,false, 3>=2, ture, 4>=1,ture (11) 答案表明：编址第2元素满足不等式(30>=20)和编址第4元素满足不等式(40>=10) (12) 答案表明：编址第2行第1列元素满足不等式(30>=20)和编址第2行第2列元素满足不等式(40>=10) Page20, ex2 (1)a, b, c的值尽管都是1，但数据类型分别为数值，字符，逻辑，注意a与c 相等，但他们不等于b (2)double(fun)输出的分别是字符a,b,s,(,x,)的ASCII码 Page20,ex3 >> r=2;p=0.5;n=12; >> T=log(r)/n/log(1+0.01*p) Page20,ex4 >> x=-2:0.05:2;f=x.^4-2.^x; >> [fmin,min_index]=min(f) 最小值最小值点编址 >> x(min_index) ans = 0.6500 最小值点 >> [f1,x1_index]=min(abs(f)) 求近似根--绝对值最小的点 f1 = 0.0328 x1_index = 24 >> x(x1_index) ans = -0.8500 >> x(x1_index)=[];f=x.^4-2.^x; 删去绝对值最小的点以求函数绝对值次小的点>> [f2,x2_index]=min(abs(f)) 求另一近似根--函数绝对值次小的点 f2 = 0.0630 x2_index = 65 >> x(x2_index) ans =
数学实验答案 Chapter 1 Page20,ex1 (5) 等于[exp(1),exp(2);exp(3),exp(4)] (7) 3=1*3, 8=2*4 (8) a为各列最小值，b为最小值所在的行号 (10) 1>=4,false, 2>=3,false, 3>=2, ture, 4>=1,ture (11) 答案表明：编址第2元素满足不等式(30>=20)和编址第4元素满足不等式(40>=10) (12) 答案表明：编址第2行第1列元素满足不等式(30>=20)和编址第2行第2列元素满足不等式(40>=10) Page20, ex2 (1)a, b, c的值尽管都是1，但数据类型分别为数值，字符，逻辑，注意a与c相等，但他们不等于b (2)double(fun)输出的分别是字符a,b,s,(,x,)的ASCII码 Page20,ex3 >> r=2;p=0.5;n=12; >> T=log(r)/n/log(1+0.01*p) Page20,ex4 >> x=-2:0.05:2;f=x.^4-2.^x; >> [fmin,min_index]=min(f) 最小值最小值点编址 >> x(min_index) ans = 0.6500 最小值点 >> [f1,x1_index]=min(abs(f)) 求近似根--绝对值最小的点 f1 = 0.0328 x1_index = 24 >> x(x1_index) ans = -0.8500 >> x(x1_index)=[];f=x.^4-2.^x; 删去绝对值最小的点以求函数绝对值次小的点 >> [f2,x2_index]=min(abs(f)) 求另一近似根--函数绝对值次小的点 f2 = 0.0630 x2_index = 65 >> x(x2_index) ans = 1.2500
编号：（）字号 《软件课程设计》报告 班级： 12级信息安全二班 姓名：李江涛 学号： 08123608 指导老师：徐慧 中国矿业大学计算机科学与技术学院 2013年 6 月
软件课程设计任务书 专业年级：信息安全二班 学生姓名：李江涛 任务下达日期：2013 年 4 月日 课程设计日期：2013 年 4 月5日至200年7月 3 日 课程设计题目：面向过程 目录 一第一阶段-----------面向过程 (4) 1 --------------------人民币凑数问题 (4) 1.1 需求分析 (4) 1.2 概要设计 (5) 1.3 详细设计与编码 (5)
1.5 用户使用说明 (6) 1.6 设计体会 (6) 2-------------------- 日期星期转换 (7) 2.1.需求分析： (7) 2.2 概要设计 (7) 2.4.调试分析 (10) 2.5.用户使用说明 (10) 2.6.测试分析 (10) 2.7.设计体会： (10) 二第二阶段------------面向对象 (11) 1--------------------学生管理系统 (11) 1.1----需求分析 (11) 1.2．概要设计 (11) 1.3.详细设计与编码 (11) 1.4 运行结果： (17) 1.5调试分析 (18) 1.6用户使用说明 (18) 1.7测试分析： (18) 1.8 实验体会 (18) 2 面向对象函数模板反向输出 (19) 1--------------------函数模板反向输出 (19) 1.1 需求分析： (19) 1.2函数模板反向输出源代码： (19) 1.4 运行结果： (21) 三第三部分----------可视化 (21) 1--------------------计算器： (21) 用你熟悉的一种可视化编程语言实现如下图所示的计算器。该计算器需要实现基础 的数学运算,如加,减,乘,除。 (21) 1.1重要程序 (21) 1.3运行结果图： (22) 四第四部分----------数据结构 (23) 1--------------------求矩阵的转置 (23) 1.1 需求分析： (23) 1.2 概要设计： (24) 1.3 详细设计与编码： (24) 1.4 运行结果： (27) 1.5 用户使用： (27) 1.6 设计体会： (27) 2--------------------数据结构统计选票 (27) 2.1 需求分析： (28) 2.2 概要设计： (28) 2.3 详细设计与编码： (28) 2.4 运行结果： (30)
实验6 空间曲线与曲面 实验目的 1．学会利用软件命令绘制空间曲线和曲面 2．通过绘制一些常见曲线、曲面去观察空间曲线和曲面的特点 3．绘制多个曲面所围成的区域以及投影区域。 实验准备 1．复习常见空间曲线的方程 2．复习常见空间曲面的方程 实验内容 1．绘制空间曲线 2．绘制空间曲面：直角坐标方程、参数方程 3．旋转曲面的生成 4．空间多个曲面的所围成的公共区域以及投影区域 软件命令 表6-1 Matlab 空间曲线及曲面绘图命令 【例6.1】绘制空间曲线 绘制空间曲线sin ,cos ,x at t y at t z ct ===，在区间09t π≤≤上的图形，这是一条锥面螺旋线，取a=10,c=3。
【程序】： t=0:pi/30:9*pi; a=10; c=3; x=a*t.*sin(t); y=a*t.*cos(t); z=c*t; plot3(x,y,z,’mo ’) 【输出】：见图6-1。 图6-1 空间曲线的绘制 【例6.2】利用多种命令绘制空间曲面 绘制二元函数z = 在区域:99,99D x y -≤≤-≤≤上的图形。 【程序】：参见Exm06Demo02.m 。 【输出】：见图6-2。 图 6-2 绘制空间曲面 【例6.3】绘制Mobius 带 Mobius 带的参数方程为 122122 cos sin cos ,[0,2],[,] sin u u x r u y r u r c v u v a b z v π=??==+∈∈??=?，， 其中,,a b c 为常数，绘制其图形。
【程序】： clear syms u v; c=4.0; a=-2*pi;b=2*pi; c=-1; d=1; x=(c+1/2*v*cos(u/2))*cos(u); y=(c+1/2*v*cos(u/2))*sin(u); z=1/2*v*sin(u/2); ezsurf(x,y,z,[a,b,c,d]) 【输出图形】 图6-2 Mobius 带 【例6.4】 画出上半球面 2222(1)x y z r ++-=与圆锥面2222()r z x y =+所围成的立体的图形及其在xoy 平面与平面y=1上的投影。 【步骤】： 【Step1】：写出它们的参数方程 上半球面参数方程：2sin cos sin sin [0,],[0,2] 1cos x r v u y r v u v u z r v ππ=?? =∈∈??=+?； 圆锥面参数方程：sin cos ,[0,2],[0,1]x y z ρθρθθπρρ=?? =∈∈??=? 【Step2】：绘制上半球面 Clear;clc;r=2/3;a1=0;a2=2*pi;b1=0;b2=pi/2;n1=40;n2=20; %准备上半球面数据 [u,v]=meshgrid(linspace(a1,a2,n1),linspace(b1,b2,n2)); x=r*sin(v).*cos(u);y=r*sin(v).*sin(u);z=1+r*cos(v); 【Step3】：绘制圆锥面 [t,s]=meshgrid(linspace(0,2*pi,20),linspace(0,1,20)); x1=s.*sin(t);y1=s.*cos(t);z1=s;surf(x1,y1,z1); 【Step4】：绘制xoy 平面内的投影：只需要球面的投影即可 z2=zeros(size(u));mesh(x,y,z2); 【Step5】：绘制曲面在y=1内的投影 y3=zeros(size(u))+1; y4=zeros(size(t))+1;% 球面、锥面 mesh(x,y3,z);mesh(x1,y4,z1); 【输出图形】：
二、判断题（“对”在题号前（）中打√×）（10分） （√）1、误差是指测量值与真值之差，即误差=测量值－真值，如此定义的误差反映的是测量值偏离真值的大小和方向，既有大小又有正负符号。 （×）2、残差（偏差）是指测量值与其算术平均值之差，它与误差定义一样。（√）3、精密度是指重复测量所得结果相互接近程度，反映的是随机误差大小的程度。 （√）4、测量不确定度是评价测量质量的一个重要指标，是指测量误差可能出现的范围。 （×）7、分光计设计了两个角游标是为了消除视差。 （×）9、调节气垫导轨水平时发现在滑块运动方向上不水平，应该先调节单脚螺钉再调节双脚螺钉。 （×）10、用一级千分尺测量某一长度（Δ仪=0.004mm），单次测量结果为N=8.000mm，用不确定度评定测量结果为N=（8.000±0.004）mm。 三、简答题（共15分） 1．示波器实验中，（1）CH1（x）输入信号频率为50Hz，CH2（y）输入信号频率为100Hz；（2）CH1（x）输入信号频率为150Hz，CH2（y）输入信号频率为50Hz；画出这两种情况下，示波器上显示的李萨如图形。（8分）
差法处理数据的优点是什么？（7分） 答：自变量应满足等间距变化的要求，且满足分组要求。（4分） 优点：充分利用数据；消除部分定值系统误差 四、计算题（20分，每题10分） 1、用1/50游标卡尺，测得某金属板的长和宽数据如下表所示，求金属板的面 解：（1）金属块长度平均值：)(02.10mm L = 长度不确定度： )(01.03/02.0mm u L == 金属块长度为：mm L 01.002.10±= %10.0=B （2分） （2）金属块宽度平均值：)(05.4mm d = 宽度不确定度： )(01.03/02.0mm u d == 金属块宽度是：mm d 01.005.4±= %20.0=B （2分） （3）面积最佳估计值：258.40mm d L S =?= 不确定度：2222222 221.0mm L d d s L s d L d L S =+=??? ????+??? ????=σσσσσ 相对百分误差：B ＝%100?S s σ=0.25％ （4分） （4）结果表达：21.06.40mm S ±= B ＝0.25％ （2分） 注：注意有效数字位数，有误者酌情扣 5、测量中的千分尺的零点误差属于已定系统误差；米尺刻度不均匀的误差属于未
有关高等数学计算过程中所涉及到的数学公式（集锦） 一、0 101101lim 0n n n m m x m a n m b a x a x a n m b x b x b n m --→∞?=??+++? =??? （系数不为0的情况） 二、重要公式（1）0sin lim 1x x x →= （2）()1 0lim 1x x x e →+= （3 ）)1n a o >= （4 ）1n = （5）lim arctan 2x x π→∞= （6）lim tan 2 x arc x π →-∞=- （7）lim arc cot 0x x →∞ = （8）lim arc cot x x π→-∞ = （9）lim 0x x e →-∞ = （10）lim x x e →+∞ =∞ （11）0 lim 1x x x + →= 三、下列常用等价无穷小关系（0x →） sin x x tan x x a r c s i n x x arctan x x 2 11c o s 2 x x - ()ln 1x x + 1x e x - 1l n x a x a - ()11x x ? +-? 四、导数的四则运算法则 ()u v u v '''±=± ()uv u v uv '''=+ 2u u v uv v v '''-??= ??? 五、基本导数公式 ⑴()0c '= ⑵1 x x μμμ-= ⑶()sin cos x x '= ⑷()cos sin x x '=- ⑸()2 tan sec x x '= ⑹()2 cot csc x x '=- ⑺()sec sec tan x x x '=? ⑻()csc csc cot x x x '=-? ⑼()x x e e '= ⑽()ln x x a a a '= ⑾()1 ln x x '=
1. 统计推断（实验12）—区间估计、假设检验 [mu,sigma,muci,sigmaci]=normfit(x,alpha); %%正态分布检验 [ht,sigt,cit]=ttest(x,mu); %%t 检验 [hz,sigz,ciz,zval]=ztest(x,mu,sigma,alpha,tail); %%z 检验 tail 默认为0 ① P297第2题：（1）分别用两个月的数据验证这种说法的可靠性； 编程：x1=[]; x2=[]; alpha=0.05; [mu1,sigma1,muci1,sigmaci1]=normfit(x1,alpha) %%一月份的均值和标准差以及其置信区间 [mu2,sigma2,muci2,sigmaci2]=normfit(x2,alpha) %%二月份的均值和标准差以及其置信区间 运行结果： （1月）mu1 =115.1500； sigma1 =3.8699； muci1 =113.3388 116.9612； sigmaci1 = 2.9430 5.6523 （2月）mu2 =120.7500； sigma2 =3.7116 muci2 =119.0129 122.4871； sigmaci2 =2.8227 5.4211 （2）分别给出1月和2月汽油价格的置信区间（05.0=α）； 编程：x1=[]; x2=[]; mu=115; alpha=0.05; [h1,sigma1,ci1]=ttest(x1,mu,alpha,0) %%一月份汽油价格的置信区间 [h2,sigma2,ci2]=ttest(x2,mu,alpha,0) %%二月份汽油价格的置信区间 运行结果：（1月）h1 =0； sigma1 =0.8642； ci1 =113.3388 116.9612 （2月）h2 =1； sigma2 =1.3241e-006； ci2 =119.0129 122.4871 （3）如何给出1月和2月汽油价格差的置信区间（05.0=α） 编程：x1=[]; x2=[]; alpha=0.05; [h1,sigma1,ci1]=normfit(x2-x1,alpha) %数据看成同一个加油的数据，其价格差和置信区间 [h2,sigma2,ci2]=ttest(x2,x1,alpha,0) %数据完全随机时，用总体的t 分布检验 运行结果：h1 = 5.6000； sigma1 =5.4715； ci1 =3.0393 8.1607 h2 =1； sigma2 =2.0582e-004； ci2 =3.0393 8.1607 结果分析：根据运行结果，我们可以知道数据完全随机时，用t 分布检验获得的结果更为合理准确。 ②第5题P297：分析：这里一件产品只有合格和不合格之分，用X=0表示合格品，X=1表示废品，可以说总体服从0-1分布，由题意得，合格率为90%，则废品率为10%，)1(2 p p p X -==σμ，方差的期望 双方的置信概率为95%，alpha=1-95%=0.05. 虽然X 不服从正态分布，但根据概率论中心极限定理，党样本容量充分大时，对样本均值 x 有，，,近似的服从)1,0(N ，由此可对总体废品率p 作如下的假设检验： . ,:0100p p H p p H ≠>≤这时应作单侧检验，取)1,0(N 的1-alpha 分位数alpha u -1，设样本的废品 率为 x ， n p p x z /)01(0--= μ，满足 alpha u z -≤1时接受 H ；否则拒绝 ) (10H H 接受 编程：n=50; %样本容量 x=7/n; %样本废品率 p0=1-0.9; %样本废品率期望
第一次练习 教学要求：熟练掌握Matlab 软件的基本命令和操作，会作二维、三维几何图形，能够用Matlab 软件解决微积分、线性代数与解析几何中的计算问题。 补充命令 vpa(x,n) 显示x 的n 位有效数字，教材102页 fplot(‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形 在下面的题目中m 为你的学号的后3位（1-9班）或4位（10班以上） 计算30sin lim x mx mx x →-与3 sin lim x mx mx x →∞- 程序： syms x limit((1001*x-sin(1001*x))/x^3,x,0) 结果： 程序： syms x limit((1001*x-sin(1001*x))/x^3,x,inf) 结果： 0 cos 1000 x mx y e =，求''y 程序： syms x diff(exp(x)*cos(1001*x/1000),2) 结果： -2001/1000000*exp(x)*cos(1001/1000*x)-1001/500*exp(x)*sin(1001/1000*x)
计算 2 2 11 00 x y e dxdy +?? 程序： dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1) 结果： 计算4 2 2 4x dx m x +? 程序： syms x int(x^4/(1000^2+4*x^2)) 结果： (10)cos , x y e mx y =求 程序： syms x diff(exp(x)*cos(1000*x),10) 结果： 给出 0x =的泰勒展式(最高次幂为4). 程序： syms x taylor(sqrt(1001/1000+x),5) 结果： Fibonacci 数列{}n x 的定义是121,1x x ==， 12,(3,4,)n n n x x x n --=+=L 用循环语句编程给出该数列的前20项（要求将结果用向量的形式给出）。 程序： x=[1,1]; for n=3:20 x(n)=x(n-1)+x(n-2); end x 结果： Columns 1 through 10 1 1 2 3 5 8 13 21 3 4 5 5 Columns 11 through 20 89 144 233 377 610 987 1597 2584 4181 6765
上海理工大学 计算机工程学院 实验报告 实验名称红细胞数目统计课程名称数字图像处理 姓名王磊学号0916020226 日期2012-11-27 地点图文信息中心成绩教师韩彦芳
一、设计内容： 主题：《红细胞数目检测》 详细说明：读入红细胞图片，通过中值滤波，开运算，闭运算，以及贴标签等方法获得细胞个数。 二、现实意义： 细胞数目检测在现实生活中的意义主要体现在医学上的作用，可通过细胞数目的检测来查看并估计病人或动物的血液中细胞数，如估测血液中红细胞、白细胞、血小板、淋巴细胞等细胞的数目，同时也可检测癌细胞的数目来查看医疗效果，根据这一系列的指标来对病人或动物进行治疗，是具有极其重要的现实作用的。 三、涉及知识内容： 1、中值滤波 2、开运算 3、闭运算 4、二值化 5、贴标签 四、实例分析及截图效果： （1）代码如下： 1、程序中定义图像变量说明 （1）Image--------------------------------------------------------------原图变量;
（2）Image_BW-------------------------------------------------------值化图象; （3）Image_BW_medfilt-------------------------中值滤波后的二值化图像; （4）Optimized_Image_BW---通过“初次二值化图像”与“中值滤波后的二值化图像”进行“或”运算优化图像效果； （5）Reverse_Image_BW--------------------------优化后二值化图象取反；（6）Filled_Image_BW----------------------已填充背景色的二进制图像；（7）Open_Image_BW--------------------------------------开运算后的图像； 2、实现代码： %-------图片前期处理------------------- %第一步：读取原图,并显示 A = imread('E:\红细胞3.png'); Image=rgb2gray(A); %RGB转化成灰度图 figure,imshow(Image); title('【原图】'); %第二步：进行二值化 Theshold = graythresh(Image); %取得图象的全局域值 Image_BW = im2bw(Image,Theshold); %二值化图象 figure,imshow(Image_BW); title('【初次二值化图像】'); %第三步二值化图像进行中值滤波 Image_BW_medfilt= medfilt2(Image_BW,[13 13]); figure,imshow(Image_BW_medfilt); title('【中值滤波后的二值化图像】'); %第四步：通过“初次二值化图像”与“中值滤波后的二值化图像”进行“或”运算优化图像效果 Optimized_Image_BW = Image_BW_medfilt|Image_BW; figure,imshow(Optimized_Image_BW); title('【进行“或”运算优化图像效果】'); %第五步：优化后二值化图象取反，保证：‘1’-〉‘白色’，‘0’-〉‘黑色’ %方便下面的操作 Reverse_Image_BW = ~Optimized_Image_BW; figure,imshow(Reverse_Image_BW); title('【优化后二值化图象取反】');