西安交大统计学考试试卷 一、单项选择题（每小题2 分，共20 分） 1.在企业统计中，下列统计标志中属于数量标志的是（C） A、文化程度 B、职业 C、月工资 D、行业 2.下列属于相对数的综合指标有（B ） A、国民收入 B、人均国民收入 C、国内生产净值 D、设备台数 3.有三个企业的年利润额分别是5000 万元、8000 万元和3900 万元，则这句话中有（B）个变量？ A、0 个 B、两个 C、1 个 D、3 个 4.下列变量中属于连续型变量的是（A ） A、身高 B、产品件数 C、企业人数 D、产品品种 5.下列各项中，属于时点指标的有（A ） A、库存额 B、总收入 C、平均收入 D、人均收入 6.典型调查是（B ）确定调查单位的 A、随机 B、主观 C、随意 D 盲目 7.总体标准差未知时总体均值的假设检验要用到（A ）： A、Z 统计量 B、t 统计量 C、统计量 D、X 统计量 8.把样本总体中全部单位数的集合称为（A ） A、样本 B、小总体 C、样本容量 D、总体容量 9.概率的取值范围是p（D ） A、大于1 B、大于－1 C、小于1 D、在0 与1 之间 10.算术平均数的离差之和等于（A ） A、零 B、1 C、－1 D、2 二、多项选择题（每小题2 分，共10 分。每题全部答对才给分，否则不计分） 1.数据的计量尺度包括（ABCD ）： A、定类尺度 B、定序尺度 C、定距尺度 D、定比尺度 E、测量尺度 2.下列属于连续型变量的有（BE ）： A、工人人数 B、商品销售额 C、商品库存额 D、商品库存量 E、总产值 3.测量变量离中趋势的指标有（ABE ） A、极差 B、平均差 C、几何平均数 D、众数 E、标准差 4.在工业企业的设备调查中（BDE ） A、工业企业是调查对象 B、工业企业的所有设备是调查对象 C、每台设备是 填报单位D、每台设备是调查单位E、每个工业企业是填报单位 5.下列平均数中，容易受数列中极端值影响的平均数有（ABC ） A、算术平均数 B、调和平均数 C、几何平均数 D、中位数 E、众数 三、判断题（在正确答案后写“对”，在错误答案后写“错”。每小题1 分，共10 分） 1、“性别”是品质标志。（对）
流体力学知识要点 第一章 流体及其主要物理性质 1. 流体的连续介质模型 a) 流体的定义：任何微小的剪切力都会导致连续变形的物质 b) 质点：含有足够多分子数，并且具有确定宏观统计特征的分子集合。 c) 连续介质模型：（欧拉）假定组成流体的最小物理实体是流体质点而不是流体分子， 即：流体是由无穷多个、无穷小的、紧密毗邻、连绵不断的流体质点所组成的一种绝无间隙的连续介质。 2. 流体的主要物理性质 a) 流体的密度：表征流体在空间某点质量的密集程度 i. 密度：'lim V V m V （'V 特征体积，此时具有统计平均特性和确定性） ii. 比容：1v b) 压缩性：当作用在一定量流体上的压强增加时，其体积将减小, 用单位压强所引起 的体积变化率表示 i. 压缩性系数b : /b dV V dp ii. 体积弹性模量E ：1 /b dp Vdp E dV V dV (Pa) v dp E d (1/)(1/)/V dp Vdp dp dp m dp dV d dV d d m 对气体： (等温 E p ；等熵 E p ，一般 1.4 ) 对液体，无明确比例 可压缩流体和不可压缩流体 液体的体积弹性模量值大，液体平衡和运动的绝大多数问题可以用不 可压缩流体解决。 气体的体积弹性模量值小，气体平衡和运动的大多数问题需要按可压 缩流体来解决。 c) 流体的粘性：是流体抵抗剪切变形或相对运动的一种固有属性，表现为流体内摩擦 i. 粘性内摩擦力产生的原因： 分子间吸引力（内聚力）产生阻力 分子不规则运动的动量交换产生的阻力 ii. 牛顿粘性实验
U U F A F A h h 牛顿内摩擦定律： /U F A h （μ 动力粘性系数，Pa ·s ） du d dy dt （d dt 角变形率） iii. 粘性系数 动力粘性系数 Pa ·s 运动粘性系数 2 /m s iv. 影响粘性的因素 压强：0p p e 正相关 温度：液体温度大粘度小 气体温度大粘度大 v. 理想流体：不具有粘性（对应粘性流体，一切实际流体都具有粘性） vi. 牛顿流体：满足牛顿内摩擦定律的流体（对应非牛顿流体，不满足牛顿内摩擦 定律） 3. 作用在流体上的力 （ 表面力 质量力） a) 表面力：作用在所取的流体分离体表面上的力。即分离体以外的流体通过接触面作 用在分离体上的力（压力，粘性力） b) 质量力：外力场作用在流体质点上的非接触力，在流体质量均匀情况下又称体积力。 质量力与外力场的强度和流体的分布有关，与它周围的微元体积无关。（重力） 4. 理想流体中的压力与方向无关 a) ,,p p x y z （即理想流体中任一点流体静压强的大小与其作用的面在空间的 方位无关，只是该点坐标的函数） b) 流体静压强的方向沿作用面的内法线方向。 第二章 流体静力学 1. 流体静压强及其特性 a) 流体静压强的方向沿作用面的内法线方向。 b) 静止流体中任一点上不论来自何方的静压强均相等，所以在静止流体中流体静压强 是空间坐标的连续函数。 2. 静止流体平衡微分方程式（欧拉平衡微分方程式） a) 1 0g a p dy du dt dy dudt dt d /
1.计算以下和式：01421181 84858616n n S n n n n ∞ =?? =--- ?++++??∑ ，要求： (1)若保留11个有效数字，给出计算结果，并评价计算的算法； (2)若要保留30个有效数字，则又将如何进行计算。 （1）题目分析 该题是对无穷级数求和，因此在使用matlab 进行累加时需要一个累加的终止条件。这里令?? ? ??+-+-+-+= 681581482184161n n n n a n n ，则 ()()1.016 1 6855844864816114851384128698161 681581482184161148113811282984161111<< ? ??? ????? ??++++++???? ????? ??++++++=??? ????? ??+-+-+-+??? ????? ??+-+-+-+=+++n n n n n n n n n n n n n n n n a a n n n n n n 故近似取其误差为1+≈k a ε，并且有m -1m -111021 21 ?=?=≈+βεk a ， （2）算法依据 使用matlab 编程时用digits 函数和vpa 函数来控制位数。 （3）Matlab 运行程序 %%保留11位有效数字 k1=11; s1=0;%用于存储这一步计算值 for n=0:50 a=(1/16^n)*(4/(8*n+1)-2/(8*n+4)-1/(8*n+5)-1/(8*n+6)); n1=n-1; if a<=0.5*10^(1-k1) break end end; for i=0:1:n1 t=(1/16^i)*(4/(8*i+1)-2/(8*i+4)-1/(8*i+5)-1/(8*i+6)); s1=s1+t; end s11=vpa(s1,k1); disp('保留11位有效数字的结果为：');disp(s11); disp('此时n 值为：');disp(n1); %%保留30位有效数字 clear all; k2=30;
一、 选择题(略) 二、 判断题(略) 三、 简答题 1.等压面是水平面的条件是什么？ :①连续介质 ② 同一介质 ③ 单一重力作用下. 2. 同一容器中装两种液体，且21ρρ?，在容器侧壁装了两根测压管。试问：图中所标明的测压管中液面位置对吗？为什么？
C (c) 盛有不同种类溶液的连通器 D C D 水 油 B B (b) 连通器被隔断 A A (a) 连通容器 解：不对，（右测压管液面要低一些，从点压强的大小分析） 3. 图中三种不同情况，试问：A-A 、B-B 、C-C 、D-D 中哪个是等压面？哪个不是等压面？为什么？ ：（ a ）A-A 是 （b ）B-B 不是 （c ）C-C 不是， D-D 是。 四、作图题(略) 五、计算题(解题思路与答案) 1. 已知某点绝对压强为80kN/m 2 ，当地大气压强p a =98kN/m 2 。试将该点绝对压强、相对压强和真空压强用水柱及水银柱表示。 解： 用水柱高表示 （1）该点绝对压强：8.16mH 2o （2）该点相对压强：-1.84mH 2o （3）该点真空压强：1.84mH 2o
用水银柱高表示 （1）该点绝对压强：599.1mm H g （2）该点相对压强：-135.4 mm H g （3）该点真空压强：135.4 mm H g 2. 一封闭水箱自由表面上气体压强p 0=25kN/m 2 ，h 1=5m ，h 2＝2m 。求A 、B 两点的静水压强。 解：由压强基本公式gh p p ρ+= 0求解 A p = 7.551 mH 2o (74 kN/m 2 ) B p = 4.551 mH 2o (44.6 kN/m 2 ) 3 如图所示为一复式水银测压计，已知m 3.21=?，m 2.12=?， m 5.23=?，m 4.14=?，m 5.15 =?(改为3.5m)。试求水箱液面上的绝 对压强0p ＝？
实验课题四曲面图与统计图 第一大题：编程作下列曲面绘图： 用平面曲线r=2+cos(t)+sin(t)，t∈(0,π)绘制旋转曲面 t=0:0.02*pi:pi; r=2+cos(t)+sin(t); cylinder(r,30) title('旋转曲面'); shading interp 用直角坐标绘制双曲抛物面曲面网线图，z2=xy (-3 axis off 用直角坐标绘制修饰过的光滑曲面曲面：z 4=sin(x )－cos(y ) x 与y 的取值在(-π,π) [x,y]=meshgrid(-pi:0.02*pi:pi); z4=sin(x)-cos(y); surf(x,y,z4); title('picture 4'); shading interp axis off 用连续函数绘图方法绘制曲面)2 s in (6522x y x z ++=,x ∈[-2pi,2pi], y ∈[-2pi,2pi],并作图形修饰。 ezsurf(@(x,y)(x^2+y^2+6*sin(2*x)),[-2*pi 2*pi -2*pi 2*pi]) title('picture 5'); shading interp axis off 第二大题：按要求作下列问题的统计图： x21是1—10的10维自然数构成的向量，y21是随机产生的10维整数向量，画出条形图。(提示bar(x,y)) x21=1:10; y21=randn(10,1); bar(x21,y21) 随机生成50维向量y22，画出分5组的数据直方图。（提示hist(y,n)） 西安交通大学2009年攻读硕士学位研究生入学考试试题 科目代码：818 科目名称：高等代数 一 （20分）计算行列式： 000 00 0000 00000n D αβαβαβαβαβαβαβαβ +++=+ + 二 （20分）已知12(0,1,0),(3,2,2)T T αα==-，是线性方程组 1231231 2321341x x x x x x ax bx cx d -+=-??++=??++=? 的两个解，求此方程组的全部解． 三 （20）当t 取什么值时，下面二次型是正定的： 222123123121323(,,)42106f x x x x x x tx x x x x x =+++++ 四（15分）设3阶实对称矩阵A 有特征值1231,1λλλ=-==，A 的属于特征值-1的特征向量1(0,1,1)T ξ=，矩阵32B A A E =-+，其中E 为3阶单位阵（下同），问： （1） 1ξ是否为B 的特征向量？求B 的所有特征值和特征向量； （2） 求矩阵B ． 五（15分）设，1200000,,,,00,,,00a c x W a a b c R W y x y z R c b z z ????????????????=∈=∈???????????????????????? （1） 求12W W +； （2） 记12W W W =+，试求空间3W 使得33()M R W W =⊕（其中3()M R 为实数域 上3阶矩阵全体），并说明理由． 六（15分）设向量组12,,,r ααα线性无关，而12,,,,,r αααβγ线性相关．证明： 要么β与γ中至少有一个可被12,,,r ααα线性表出，要么12,,,,r αααβ与12,,,,r αααγ等价． 七（15分）设A 为(1)n n ?+阶常数矩阵，X 为(1)n n +?阶未知数矩阵．试证明矩阵方程AX E =有解的充要条件为()r A n =． 八（10）若12,αα是数域F 上的二维线性空间2()V F 的基，σ和τ是2()V F 上的线性变换，且满足 112212121212,,(),()σαβσαβτααββτααββ==+=+-=- 试证：στ=． 九（10）设A 和B 是两个n 阶实正交矩阵，并且det()det()A B =-．证明 ()r A B n +<． 十（10分）证明A 可与一个对角矩阵相似的充要条件是：对于A 的任意特征值i λ，方程组 2()0i E A X λ-=与()0i E A X λ-= 是同解的，其中11(,,,)n n X x x x =．需要更多试题请https://www.doczj.com/doc/0510289380.html,/exam.taoba -//maths :http 高等代数试题分数分布： 行列式：20分（1）； 线性方程组：35分（2）； 矩阵：15分（1）； 二次型：20分（1）； 线性空间：15分（1）； 欧几里得空间：10分（1） 线性变换：35分（3） 第一章流体及其主要物理性质 一概念 流体：在任何微小剪切力持续作用下连续变形的物体。 连续介质模型：假定组成流体的最小物质实体是流体质点，流体是由无限多个流体质点组成，质点之间不存在间隙。 适用条件：分子平均自由程远小于流动问题特征尺寸。 不适用条件：稀薄气体，激波层内等。 粘性：流体抵抗剪切变形（相对运动）的一种属性。 流体层间无相对运动时不表现粘性 牛顿平板实验： 两板间流体速度： 剪切力，即： 或： 式中与板间流体的种类、流体的温度、压强有关，成为液体的动力粘性系数，简称粘性系数。 流体做任意层状流动： 牛顿内摩擦定律的数学表达式 式中是角变形率或角变形速度。 凡符合牛顿内摩擦定律的流体称为牛顿流体。 产生粘性机理： 同一种流体的动力粘性系数与流体的温度有很大的关系，而受压强的影响很小。液体与气体的产生粘性机理不一样，液体的粘性主要取决于分子间的距离和分子间的吸引力，故温度升 高粘性下降；气体的粘性主要取决于分子气体热运动所产生的动量交换，故温度升高，其粘性增大。 在国际单位中： μ反应流体真实粘性的大小 运动粘性系数： 物理单位是： 粘性系数等于零的流体称为理想流体或无粘流体。 工程上常用体积弹性模量衡量流体的可压缩性，体积弹性模量定义为： 体积弹性模量的量纲和压强相同，是或。 流体的压缩性越大，则越大，即越小；反之，可压缩性越小，则越 小，即越大。 体积弹性模量又可以表示为： 等温体积弹性模量： 等熵体积弹性模量： 由： 当很小或者很大，由或者二者兼得，则此时流体的密度相对变化量就很小。如果忽 略流体密度的变化，不考虑流体的可压缩性的影响，这种简化的模型称为不可压缩流体，其密度可视为常量；反之，考虑密度为变量或压缩性影响的流体，称为可压缩流体。 不可压缩流体：液体低速流动的气体 可压缩流体：气体水下爆炸和水锤现象的液体 实验课题五线性代数 第一大题：创建矩阵： 1.1 用元素输入法创建矩阵 ??? ???? ??-=34063689 864275311A ?????? ? ? ?--=96 5 214760384 32532A A1=[1 3 5 7;2 4 6 8;9 8 6 3;-6 0 4 3] A2=[3 5 -2 3;4 8 3 0;6 7 4 -1;2 5 6 9] 1.2 创建符号元素矩阵 ???? ? ?=54 3 2 15432 13y y y y y x x x x x A ??? ? ??+=)cos(1)sin(42x x x x A A3=sym('[x1 x2 x3 x4 x5;y1 y2 y3 y4 y5]') A4=sym('[sin(x) x^2;1+x cos(x)]') 1.3 生成4阶随机整数矩阵B B=rand(4) 1.4 由向量t=[2 3 4 2 5 3]生成范德蒙矩阵F t=[2 3 4 2 5 3]; F=vander(t) 1.5 输入4阶幻方阵C C=magic(4) 1.6 用函数创建矩阵：4阶零矩阵Q ; 4阶单位矩阵E ; 4阶全壹矩阵N Q=zeros(4) E=eye(4) N=ones(4) 1.7 用前面题目中生成的矩阵构造8×12阶大矩阵： ???? ? ?=16A C N Q E B A A6=[B E Q;N C A1] 第二大题：向量计算： 2.1计算：a21是A1的列最大元素构成的向量，并列出所在位置。提示：[a21,i]=max(A1) a22是A1的列最小元素构成的向量，并列出所在位置. a23是A1的列平均值构成的向.， a24是A1的列中值数构成的向量. a25是A1的列元素的标准差构成的向量. a26是A1的列元素和构成的向量. [a21,i]=max(A1) [a22,j]=min(A1) a23=mean(A1) a24=median(A1) a25=std(A1) a26=sum(A1) 2.2计算a27=A1+A2；a28=A1×A2 a27=A1+A2 a28=A1.*A2 2.3取矩阵A2的一、三行与二、三列的交叉元素做子矩阵A29. A29=A2([1,3],[2,3]) 第三大题：矩阵运算 3.1生成6阶随机整数矩阵A A=fix(15*rand(6)) 3.2作A31等于A的转置;作A32等于A的行列式;作A33等于A的秩。 A31=A' A32=det(A) A33=rank(A) 3.3判断A是否可逆.若A可逆，作A34等于A的逆，否则输出‘A不可逆’。 if det(A)==0 disp('A不可逆'); else A34=inv(A) end 西安交通大学 数据结构与算法课程实验 实验名称：数据结构与算法课程专题实验 所属学院：电信学院 专业班级：计算机32班 小组成员： 指导老师：赵仲孟教授 实验一背包问题的求解 1.问题描述 假设有一个能装入总体积为T的背包和n件体积分别为w1,w2,…w n的物品，能否从n件物品中挑选若干件恰好装满背包，即使w1+w2+…+w m=T，要求找出所有满足上述条件的解。 例如：当T=10，各件物品的体积{1,8,4,3,5,2}时，可找到下列4组解： (1,4,3,2) (1,4,5) (8,2) (3,5,2)。 2.实现提示 可利用回溯法的设计思想来解决背包问题。首先，将物品排成一列，然后，顺序选取物品装入背包，若已选取第i件物品后未满，则继续选取第i+1件，若该件物品“太大”不能装入，则弃之，继续选取下一件，直至背包装满为止。 如果在剩余的物品中找不到合适的物品以填满背包，则说明“刚刚”装入的物品“不合适”，应将它取出“弃之一边”，继续再从“它之后”的物品中选取，如此重复，直到求得满足条件的解，或者无解。 由于回溯求解的规则是“后进先出”，自然要用到“栈”。 3.问题分析 1、设计基础 后进先出，用到栈结构。 2、分析设计课题的要求，要求编程实现以下功能： a．从n件物品中挑选若干件恰好装满背包 b. 要求找出所有满足上述条件的解，例如：当T=10，各件物品的体积{1，8，4， 3，5，2}时，可找到下列4组解：（1，4，3，2）、（1，4，5）、（8，2）、（3，5，2）3，要使物品价值最高，即p1*x1+p2*x1+...+pi*xi(其1<=i<=n，x取0或1，取1表示选取物品i) 取得最大值。在该问题中需要决定x1 .. xn的值。假设按i = 1，2，...，n 的次序来确定xi 的值。如果置x1 = 0，则问题转变为相对于其余物品(即物品2，3，.，n)，背包容量仍为c 的背包问题。若置x1 = 1，问题就变为关于最大背包容量为c-w1 的问题。现设r={c，c-w1} 为剩余的背包容量。在第一次决策之后，剩下的问题便是考虑背包容量为r 时的决策。不管x1 是0或是1，[x2 ，.，xn ] 必须是第一次决策之后的一个最优方案。也就是说在此问题中，最优决策序列由最优决策子序列组成。这样就满足了动态规划的程序设计条件。 4.问题实现 代码1： #include"iostream" using namespace std; class Link{ public: int m; Link *next; Link(int a=0,Link *b=NULL){ m=a; next=b; } }; class LStack{ private: Link *top; 西南交通大学限修课数学实验题目及答案六 实验课题六一元微积分 第一大题函数运算 1.用程序集m 文件中定义函数： 键盘输入自变量x ,由下列函数 求函数值：f 1 (12) f 1 (-32) function y=f1(x) if x>0 y=4*x^3+5*sqrt(x)-7 else y=x^2+sin(x) end end 2. 用函数m 文件定义函数f 2 ???<+≥+=06)5sin(0 3232x x x x x e f x 求f 2(-6) f 2(11) function y=f2(x) if x<0 y=sin(5*x)+6*x^3 else y=exp(2*x)+3*x ???≤+>-+=0 )sin(0 754123x x x x x x f 313-+=x x f end end 3.已知 求 其反函 数 syms x f3=(1+x)/(x-3); g=finverse(f3) %g =(3*x + 1)/(x - 1) 4.已知： 92847 653423234-++=+-+=x x x g x x x f 做函数运算：u1 = f 4+ g 4 ; u2 = f 4 – g 4 ; u3 = f 4 * g 4 ; u4 = f 4 / g 4 u5=)(4)(4x g x f ,u6=()()x g f 44 syms x f4=3*x^4+5*x^3-6*x^2+7 g4=8*x^3+2*x^2+x-9 u1=f4+g4 u2=f4-g4 u3=f4*g4 u4=f4/g4 u5=f4^g4 u6=compose(f4,g4) %u1 =3*x^4 + 13*x^3 - 4*x^2 + x - 2 %u2 =3*x^4 - 3*x^3 - 8*x^2 - x + 16 %u3 =(3*x^4 + 5*x^3 - 6*x^2 + 7)*(8*x^3 + 2*x^2 + x - 9) %u4 =(3*x^4 + 5*x^3 - 6*x^2 + 7)/(8*x^3 + 2*x^2 + x - 9) %u5 =(3*x^4 + 5*x^3 - 6*x^2 + 7)^(8*x^3 + 2*x^2 + x - 9) %u6 =5*(8*x^3 + 2*x^2 + x - 9)^3 - 6*(8*x^3 + 2*x^2 + x - 9)^2 + 3*(8*x^3 + 高等代数与解析几何同济答案 【篇一：大学所有课程课后答案】 资料打开方法：按住 ctrl键，在你需要的资料上用鼠标左键单击 资料搜索方法：ctrl+f 输入关键词查找你要的资料 【数学】 o o o o o o o o o o o o o o o o o 习题答案 o o o o o o o o o o o o o o o o o o o o 【计算机/网络/信息】 o o o o 【经济/金融/营销/管理/电子商务】 o o o o o o o o o o o o 【 o o o o o 【篇二：各门课程课后答案】 式]《会计学原理》同步练习题答案 [word格式]《成本会计》习题及答案（自学推荐，23页） [word格式]《成本会计》配套习题集参考答案 [word格式]《实用成本会计》习题答案 [word格式]《会计电算化》教材习题答案（09年） [jpg格式]会计从业《基础会计》课后答案 [word格式]《现代西方经济学（微观经济学）》笔记与课后习题详解（第3版，宋承先） [word格式]《宏观经济学》习题答案（第七版，多恩布什） [word格式]《国际贸易》课后习题答案（海闻 p.林德特王新奎） [pdf格式]《西方经济学》习题答案（第三版，高鸿业）可直接打印 [word格式]《金融工程》课后题答案（郑振龙版） [word格式]《宏观经济学》课后答案（布兰查德版） [jpg格式]《投资学》课后习题答案（英文版，牛逼版） [pdf格式]《投资学》课后习题答案（博迪，第四版） [word格式]《微观经济学》课后答案（高鸿业版） [word格式]《公司理财》课后答案（英文版，第六版） [word格式]《国际经济学》教师手册及课后习题答案（克鲁格曼，第六版） 计算方法上机报告 姓名： 学号： 班级： 目录 题目一------------------------------------------------------------------------------------------ - 4 - 1.1题目内容 ---------------------------------------------------------------------------- - 4 - 1.2算法思想 ---------------------------------------------------------------------------- - 4 - 1.3Matlab源程序----------------------------------------------------------------------- - 5 - 1.4计算结果及总结 ------------------------------------------------------------------- - 5 - 题目二------------------------------------------------------------------------------------------ - 7 - 2.1题目内容 ---------------------------------------------------------------------------- - 7 - 2.2算法思想 ---------------------------------------------------------------------------- - 7 - 2.3 Matlab源程序---------------------------------------------------------------------- - 8 - 2.4计算结果及总结 ------------------------------------------------------------------- - 9 - 题目三----------------------------------------------------------------------------------------- - 11 - 3.1题目内容 --------------------------------------------------------------------------- - 11 - 3.2算法思想 --------------------------------------------------------------------------- - 11 - 3.3Matlab源程序---------------------------------------------------------------------- - 13 - 3.4计算结果及总结 ------------------------------------------------------------------ - 14 - 题目四----------------------------------------------------------------------------------------- - 15 - 4.1题目内容 --------------------------------------------------------------------------- - 15 - 4.2算法思想 --------------------------------------------------------------------------- - 15 - 4.3Matlab源程序---------------------------------------------------------------------- - 15 - 4.4计算结果及总结 ------------------------------------------------------------------ - 16 - 题目五----------------------------------------------------------------------------------------- - 18 - 计算方法（B）实验报告 姓名： 学号： 学院： 专业： 实验一 三对角方程组Tx f =的求解 一、 实验目的 掌握三对角方程组Tx f =求解的方法。 二、 实验内容 求三对角方程组Tx f =的解，其中： 4 -1 -1 4 -1 -1 4 1 -1 4T ????????=?? ?? ???? ， 3223f ?? ? ? ?= ? ? ??? 三、 算法组织 设系数矩阵为三对角矩阵 11222333111 b c a b c a b c a b c b n n n n T ---???????? =?????? ?????? 则方程组Tx f =称为三对角方程组。 设矩阵T 非奇异，T 可分解为T=LU ，其中L 为下三角矩阵，U 为单位上三角矩阵，记 1 1 212 313 1 1 1111 ,11n n n n n r l r l r L U l r l μμμμμ---???? ? ? ? ? ? ?== ? ? ? ? ? ? ? ? ? ?? ? ? ? 可先依次求出,L U 中的元素后，令Ux y =，先求解下三角方程组Ly f =得出 y ，再求解上三角方程组Ux y =。 追赶法的算法组织如下： 1.输入三对角矩阵T 和右端向量f ； 2.将Tx f =压缩为四个一维数组{}{}{}{}i i i i a b c d 、、、，{}{}{}i i i a b c 、、是T 的三对角线性方程组的三个对角，{}i d 是右端向量。将分解矩阵压缩为三个一维数组 {}{}{}i i i l r μ、、。 3.对T 做Crout 分解（也可以用Doolittle 分解）导出追赶法的计算步骤如下： 1111,b r c μ== for 2i n = 111, , ,i i i i i i i i i i i i i l a b a r r c y d l y μμ---==-==- end 4.回代求解x /n n n x y μ= for 11i n =- 1()/i i i i i x y c x μ+=- end 5. 停止，输出结果。 四、 MATLAB 程序 MATLAB 程序见附件1. 五、 结果及分析 实验结果为： (1.0000 1.0000 1.0000 1.0000)T x = 数学实验报告 制作成员班级学号 2011年6月12日 培养容器温度变化率模型 一、实验目的 利用matlab软件估测培养容器温度变化率 二、实验问题 现在大棚技术越来越好，能够将温度控制在一定温度范围内。为利用这种优势，实验室现在需要培植某种适于在8.16℃到10.74℃下能够快速长大的甜菜品种。为达到实验所需温度，又尽可能地节约成本，研究所决定使用如下方式控制培养容器的温度：1，每天加热一次或两次，每次约两小时； 2，当温度降至8.16℃时，加热装置开始工作；当温度达到10.74℃时，加热装置停止工作。 已知实验的时间是冬天，实验室为了其它实验的需要已经将实验室的温度大致稳定在0℃。下表记录的是该培养容器某一天的温度 时间(h)温度（℃）时间（h）温度（℃）09.68 1.849.31 0.929.45 2.959.13 3.878.981 4.989.65 4.988.811 5.909.41 5.908.691 6.839.18 7.008.5217.938.92 7.938.3919.048.66 8.978.2219.968.43 9.89加热装置工作20.848.22 10.93加热装置工作22.02加热装置工作10.9510.8222.96加热装置工作12.0310.5023.8810.59 12.9510.2124.9910.35 13.889.9425.9110.18 三、建立数学模型 1，分析：由物理学中的傅利叶传热定律知温度变化率只取决于温度 差，与温度本身无关。因为培养容器最低温度和最高温度分别是：8.16℃和10.74℃；即最低温度差和最高温度差分别是：8.16℃和10.74℃。而且，16.8/74.10≈1.1467，约为1，故可以忽略温度对温度变化率的影响2， 将温度变化率看成是时间的连续函数，为计算简单，不妨将温度变化率定义成单位时间温度变化的多少，即温度对时间连续变化的绝对值（温度是下降的），得到结果后再乘以一系数即可。 四、问题求解和程序设计流程1）温度变化率的估计方法 根据上表的数据，利用matlab 做出温度-时间散点图如下： 下面计算温度变化率与时间的关系。由图选择将数据分三段，然后对每一段数据做如下处理：设某段数据为{(0x ,0y ),(1x ,1y ),(2x , 2y ),…,(n x ,n y )},相邻数据中点的平均温度变化率采取公式： 温度变化率=（左端点的温度-右端点的温度）/区间长度算得即：v( 2 1i i x x ++)=(1+-i i y y )/(i i x x - +1). 每段首尾点的温度变化率采用下面的公式计算：v(0x )=(30y -41y +2y )/(2x -0x )v(n x )=(3n y -41+n y +2+n y )/(n x -2-n x ) 计算方法A 上机大作业 1. 共轭梯度法求解线性方程组 算法原理：由定理3.4.1可知系数矩阵A 是对称正定矩阵的线性方程组Ax=b 的解与求解二次函数1()2 T T f x x Ax b x =-极小点具有等价性，所以可以利用共轭梯度法求解1()2 T T f x x Ax b x = -的极小点来达到求解Ax=b 的目的。 共轭梯度法在形式上具有迭代法的特征，在给定初始值情况下，根据迭代公式： (1)()()k k k k x x d α+=+ 产生的迭代序列(1)(2)(3)x x x ，，，... 在无舍入误差假定下，最多经过n 次迭代，就可求得()f x 的最小值，也就是方程Ax=b 的解。 首先导出最佳步长k α的计算式。 假设迭代点()k x 和搜索方向()k d 已经给定，便可以通过()()()() k k f x d φαα=+的极小化 ()()min ()()k k f x d φαα=+ 来求得，根据多元复合函数的求导法则得: ()()()'()()k k T k f x d d φαα=?+ 令'()0φα=，得到: ()() ()()k T k k k T k r d d Ad α=，其中()()k k r b Ax =- 然后确定搜索方向()k d 。给定初始向量(0)x 后，由于负梯度方向是函数下降最快的方向，故第一次迭代取搜索方向(0) (0)(0)(0)()d r f x b Ax ==-?=-。令 (1)(0)00x x d α=+ 其中(0)(0)0(0)(0) T T r d d Ad α=。第二次迭代时，从(1) x 出发的搜索方向不再取(1)r ，而是选取(1) (1)(0)0d r d β=+，使得(1)d 与(0)d 是关于矩阵A 的共轭向量，由此可 求得参数0β: 实验报告（三） 完成人：L.W.Yohann 注：本次实验主要学习了用MATLAB循环结构、选择结构进行编程，在学习完成后小组对65页的上机练习题进行了 程序编辑和运行。 1.使用for循环求和. 解：在编辑窗口输入： clear;clc; n=20;s=0; for i=1:n s=s+((i^2+3*i)/(2*i+1)); fprintf('i=%.0f,s=%.5f\n',i,s) end 并保存，命名为lab1； 在命令窗口中输入lab1,得： i=1,s=1.33333 i=2,s=3.33333 i=3,s=5.90476 i=4,s=9.01587 i=5,s=12.65224 i=6,s=16.80608 i=7,s=21.47275 i=8,s=26.64922 i=9,s=32.33343 i=10,s=38.52391 i=11,s=45.21956 i=12,s=52.41956 i=13,s=60.12326 i=14,s=68.33016 i=15,s=77.03984 i=16,s=86.25196 i=17,s=95.96624 i=18,s=106.18246 i=19,s=116.90041 i=20,s=128.11992 2.编写程序，通过键盘输入一组数，找出其中的最大数和最 小数. 3.解：在编辑窗口输入： a=input('请输入一组数x（用中括号括起来）：'); n=length(a); m=a(1);M=a(1); for i=2:n if a(i) 实验二、三：quartusⅡ原理图设计1.实验原理图 2.实验仿真波形 实验四：Verilog描述组合逻辑电路1.一位数值比较器 1.1源代码 module compare(a_gt,a_eq,a_lt,a,b); input a,b; output a_gt,a_eq,a_lt; assign a_gt=a&~b; assign a_eq=a&b|~a&~b; assign a_lt=~a&b; endmodule 1.2代码生成原理图 2.七段译码器 2.1源代码 module decode4_7(codeout,indec); input[3:0] indec; output[6:0] codeout; reg[6:0] codeout; always@(indec) begin case(indec) 4'd0:codeout=7'b1111110; 4'd1:codeout=7'b0110000; 4'd2:codeout=7'b1101101; 4'd3:codeout=7'b1111001; 4'd4:codeout=7'b0110011; 4'd5:codeout=7'b1011011; 4'd6:codeout=7'b1011111; 4'd7:codeout=7'b1110000; 4'd8:codeout=7'b1111111; 4'd9:codeout=7'b1111011; default: codeout=7'b1001111; endcase end endmodule 2.2代码生成原理图 3.总原理图 4.实验仿真波形图 西安交通大学入学测试机考 专升本高数（一）模拟题1、题目Z1-1（2）（） A．A B．B C．C D．D 标准答案：B 2、题目1-1（2）（） A．A B．B C．C D．D 标准答案：C 3、题目1－2（2）（） A．A B．B C．C D．D 标准答案：A 4、题目1－3（2）（） A．A B．B C．C D．D 标准答案：B 5、题目6－1：（2）（） A．A B．B C．C D．D 标准答案：B 6、题目1－4（2）（） A．A B．B C．C D．D 标准答案：D 7、题目1－5（2）（） A．A B．B C．C D．D 标准答案：C 8、题目1－6（2）（） A．A B．B C．C D．D 标准答案：D 9、题目1－7（2）（） A．A B．B C．C D．D 标准答案：B 10、题目1－8（2）（） A．A B．B C．C D．D 标准答案：A 11、题目1－9（2）（） A．A B．B C．C D．D 标准答案：B 12、题目1－10（2）（） A．A B．B C．C D．D 标准答案：B 13、题目6－2：（2）（） A．A B．B C．C D．D 标准答案：C 14、题目2－1（2）（） A．A B．B C．C D．D 标准答案：D 15、题目2－2（2）（） A．A B．B C．C D．D 标准答案：C 16、题目2－3（2）（） A．A B．B C．C D．D 标准答案：D 17、题目6－3：（2）（） A．A B．B C．C D．D 标准答案：C 18、题目2－4（2）（） A．A B．B C．C D．D 标准答案：B 19、题目6－4：（2）（） A．A B．B C．C D．D 标准答案：C 20、题目2－5（2）（） A．A B．B C．C D．D 标准答案：A 21、题目2－6（2）（）2009西安交通大学高等代数考研真题
西安交通大学-流体力学-期末总复习
西南交通大学限修课数学实验题目及答案五
数据结构与算法分析专题实验-西安交大-赵仲孟
西南交通大学限修课数学实验题目及答案六
高等代数与解析几何同济答案
西安交通大学计算方法B大作业
西交大计算方法上机报告
matlab数学实验报告5
西交计算方法A上机大作业
西安交通大学数学实验报告(用MATLAB绘制二维、三维图形)(MATLAB循环结构、选择结构)
(精选)西南交大数电实验报告
西安交通大学入学测试机考《高等数学一(专升本)》模拟题及答案
相关主题
文本预览