数学建模与数学实验 作业
数学实验与数学建模 实验二
数学实验与数学建模实验二2.圆钢原材料每根长5.5米，现需要A,B,C三种圆钢材料，长度分别为3.1m,2.1m,1.2m,数量分别为100根，200根，400根，试安排下料方式，使所需圆钢原材料的总数最少。
约束条件为:目标函数：则列lingo程序如下：min=x1+x2+x3+x4+x5;x1+x2>=100;x1+x3+2*x4>=200;2*x2+2*x3+x4+4*x5>=400;@gin(x1);@gin(x2);@gin(x3);@gin(x4);@gin(x5);End运行结果如下：结合lingo数据得出结论：方案一和三没有采用，方案二和四用去100跟原材料，方案五用去25跟原材料，一共用去225根原材料，即为最优。
3．住宅小区服务中心选址：某地新建一个生活住宅区，共有20栋住宅楼，小区内所有道路都是东西或南北走向，开发商拟在小区内修建一个服务中心，地址选在离所有楼房的总路程最小的地方。
为保证建筑物之间有足够的空间，服务中心的位置与其他楼房位置之间的距离不能少于30米（已经考虑了所有建筑的占地面积），请你确定服务中心的位置。
设初始点X0=[20,20],设（ai,bi）(i=1,......,20)为第i栋住宅楼的坐标：a=[29.74 4.9 69.32 65.0 98.3 55.27 40.0 19.8 62.5 73.3 37.58 0.98 41.98 75.37 79.38 92.0 84.47 36.77 62.08 73.13],b=[19.39 90.48 56,92 63.18 23.44 54.88 93.16 33.5 65.5 39.19 62.73 69.9 39.72 41.37 65.52 43.5 34.6 75.2 12.32 86.7]。
约束条件：目标函数：解：列出lingo式子如下：model:sets:zl/1..20/:x,y,e;endsetsdata:x=29.74,4.9,69.32,65.0,98.3,55.27,40.0,19.8,62.5,73.3,37.58,0.98,41.98,75.37,79.38,92.0,84.47,36.77,62.08,7 3.13;y=19.39,90.48,56.92,63.18,23.44,54.88,93.16,33.5,65.5,39.19,62.73,69.9,39.72,41.37,65.52,43.5,34.6,75.2,12.32,86.7;e=20,20;enddatamin=@sum(zl(i):(((x(i)-px)^2)^(1/2)+((y(i)-py)^2)^(1/2)));@for(zl(i):(x(i)-px)^2+(y(i)-py)^2>=900);end运行结果如下：得出结论为：当服务中心的坐标为[1.281228,9.897984]时，离所有的楼房的总路程最小。
数学建模实验报告
数学建模实验报告一、实验目的1.通过具体的题目实例, 使学生理解数学建模的基本思想和方法, 掌握数学建模分析和解决的基本过程。
2、培养学生主动探索、努力进取的的学风, 增强学生的应用意识和创新能力, 为今后从事科研工作打下初步的基础。
二、实验题目（一）题目一1.题目: 电梯问题有r个人在一楼进入电梯, 楼上有n层。
设每个乘客在任何一层楼出电梯的概率相同, 试建立一个概率模型, 求直到电梯中的乘客下完时, 电梯需停次数的数学期望。
2.问题分析（1）由于每位乘客在任何一层楼出电梯的概率相同, 且各种可能的情况众多且复杂, 难于推导。
所以选择采用计算机模拟的方法, 求得近似结果。
（2）通过增加试验次数, 使近似解越来越接近真实情况。
3.模型建立建立一个n*r的二维随机矩阵, 该矩阵每列元素中只有一个为1, 其余都为0, 这代表每个乘客在对应的楼层下电梯（因为每个乘客只会在某一层下, 故没列只有一个1）。
而每行中1的个数代表在该楼层下的乘客的人数。
再建立一个有n个元素的一位数组, 数组中只有0和1,其中1代表该层有人下, 0代表该层没人下。
例如:给定n=8;r=6（楼8层, 乘了6个人）,则建立的二维随机矩阵及与之相关的应建立的一维数组为：m =0 0 1 0 0 01 0 0 0 0 00 0 0 0 0 00 1 0 0 0 00 0 0 0 0 00 0 0 0 0 10 0 0 0 1 00 0 0 1 0 0c = 1 1 0 1 0 1 1 14.解决方法（MATLAB程序代码）:n=10;r=10;d=1000;a=0;for l=1:dm=full(sparse(randint(1,r,[1,n]),1:r,1,n,r));c=zeros(n,1);for i=1:nfor j=1:rif m(i,j)==1c(j)=1;break;endcontinue;endends=0;for x=1:nif c(x)==1s=s+1;endcontinue;enda=a+s;enda/d5.实验结果ans = 6.5150 那么, 当楼高11层, 乘坐10人时, 电梯需停次数的数学期望为6.5150。
数学建模与数学实验习题答案
数学建模与数学实验习题答案数学建模与数学实验习题答案数学建模和数学实验习题是数学学习中的重要组成部分，通过这些习题，我们可以更好地理解和应用数学知识。
本文将介绍数学建模和数学实验习题的一些答案和解题方法，帮助读者更好地掌握数学学习。
一、数学建模数学建模是将数学方法和技巧应用于实际问题的过程。
在数学建模中，我们需要将实际问题抽象为数学模型，并通过数学方法进行求解和分析。
下面是一个简单的数学建模问题和其解题过程。
问题：某工厂生产产品A和产品B，每天的产量分别为x和y。
产品A的生产成本为10x+20y，产品B的生产成本为15x+10y。
如果工厂每天的总成本不超过5000元，且产品A的产量必须大于产品B的产量，求工厂一天最多能生产多少个产品。
解题过程：首先，我们需要建立数学模型来描述这个问题。
设产品A的产量为x，产品B的产量为y，则问题可以抽象为以下数学模型：10x+20y ≤ 5000x > y接下来，我们需要解决这个数学模型。
首先，我们可以通过图像法来解决这个问题。
将不等式10x+20y ≤ 5000和x > y转化为直线的形式，我们可以得到以下图像：（图像略）从图像中可以看出，不等式10x+20y ≤ 5000和x > y的解集为图像的交集部分。
通过观察图像，我们可以发现交集部分的最大值为x=250，y=125。
因此，工厂一天最多能生产250个产品A和125个产品B。
除了图像法，我们还可以通过代数法来解决这个问题。
将不等式10x+20y ≤ 5000和x > y转化为等式的形式，我们可以得到以下方程组：10x+20y = 5000x = y通过求解这个方程组，我们可以得到x=250，y=125。
因此，工厂一天最多能生产250个产品A和125个产品B。
二、数学实验习题数学实验习题是通过实际操作和实验来学习数学知识和技巧的一种方式。
下面是一个关于概率的数学实验习题和其答案。
习题：一枚硬币抛掷10次，求出现正面的次数为偶数的概率。
数学建模 -实验报告1
������������⁄������������ = ������������(1 − (������ + ������)) − ������1������∗������,
（4 − 3）
������������∗⁄������������ = −������1������∗������ + ������2������
二、 问题分析
建立肿瘤细胞增长模型时，我们可以从自由增长模型开始分析，引进 Logistic 阻滞增长模型，构成肿瘤细胞增长初步框架。再者肿瘤细胞不同于普 通细胞，其生长受到人体自身免疫系统的制约。于是综合考虑正常细胞转化，癌 细胞增殖，癌细胞死亡，癌细胞被效应细胞消除等情况，建立动力学方程。并对 模型进行适当简化求解。在放射治疗方案的设计中，我们可以引入放射生物学中 广泛接受的 LQ 模型对问题进行分析，由于放疗对人体伤害相当大，因此我们采 取分次逐次放疗的方式进行治疗。我们具体分两种情形进行讨论，一是在总剂量 一定的条件下，不同的分次剂量组合对生物效应的影响；二是在产生相同生物效 应的情况下，分析最优的分次剂量组合。
易算出癌细胞转入活动期已有 300 多天，故如何在早期发现癌症是攻克癌症的关键之一 （2）手术治疗常不能割去所有癌细胞，故有时需进行放射疗法。射线强度太小无法杀
死癌细胞，太强病人身体又吃不消且会使病人免疫功能下降。一次照射不可能杀死全部癌细 胞，请设计一个可行的治疗方案（医生认为当体内癌细胞数小于 100000 个时即可凭借体内 免疫系统杀灭）。
进一步简化，根据(4-4)，(4-5)式可知，效应细胞������∗和复合物������有出有进．假 设出入保持平衡，则有
������ + ������∗ = C (C 为常数)
《数学建模与数学实验》(第三版)6.5习题作业2
1.根据物理定律K K K R I V =,R I P 2=,建立如下模型：（1）：目标函数为：∑==412k k k R IP 约束条件⎪⎪⎩⎪⎪⎨⎧===++=≤≤8,6,41023213214I I I I I I I I R I k k k 1)直接计算求解183214=++=I I I I()K K k K K K K R I I R I P ∑∑====41412min min=K k K V I∑=41min现在K V 一定，要想求P 的最小值，只需K I 最小即可。
又因为K I 已知，代入数据即可求解。
即218282624min 44332211⨯+⨯+⨯+⨯=+++=V I V I V I V I P2）有K I 已知及K V 的取值范围，可得K R 的取值范围。
min =I1^2*R1+I2^2*R2+I3^2*R3+I4^2*R4;I1=4;I2=6;I3=8;I4=18;R1>=1/2;R2>=1/3;R3>=1/4;R4>=1/9;R1<=5/2;R2<=5/3;R3<=5/4;R4<=5/9;EndGlobal optimal solution found.Objective value: 72.00000Total solver iterations: 0Variable Value Reduced Cost I1 4.000000 0.000000 R1 0.5000000 0.000000 I2 6.000000 0.000000 R2 0.3333333 0.000000 I3 8.000000 0.000000 R3 0.2500000 0.000000 I4 18.00000 0.000000 R4 0.1111111 0.000000Row Slack or Surplus Dual Price 1 72.00000 -1.000000 2 0.000000 -4.000122 3 0.000000 -4.000081 4 0.000000 -4.000061 5 0.000000 -4.000027 6 0.000000 -16.00000 7 0.000000 -36.00000 8 0.000000 -64.00000 9 0.000000 -324.0000 10 2.000000 0.000000 11 1.333333 0.000000 12 1.000000 0.000000 13 0.4444444 0.000000（2）：目标函数：∑==412k k k R I P 约束条件为：⎪⎪⎩⎪⎪⎨⎧≤≤===≤≤++=628,6,4263213214k kk kI V V V V R V I I II1)183214=++=I I I I()K K k K KK K R I I R I P ∑∑====41412min min=K k K V I ∑=41min)min(44332211V I V I V I V I P +++=要使P 最小，取4V =0，则)min(332211V I V I V I P ++=现在K V 一定，要想求P 的最小值，只需K I 最小即可。
数学建模计算实验
学时：4学时 实验目的：掌握用Lindo求解线性规划问题的方法，能够阅读Lindo结果 报告。
实验内容：
解：
实例2：求解书本上P130的习题1。列出线性规划模型，然后用
Lindo求解，根据结果报告得出解决方案。
投资规划问题
某银行经理计划用一笔资金进行有价证券的投资，可供购进的证券
以及其信用等级、到期年限、收益如下表所示。按照规定，市政证券的
解：设投资证券A,B,C,D的金额分别为
（百万元），按照
规定限制1000万元的资金约束，则线性规划模型为：
0.043 +0.054*0.5 +0.050*0.5 +0.044*0.5 +0.045
实验三：用Lingo求解非线性规划问题
学时：2学时 实验目的：掌握用Lingo求解非线性规划问题的方法。 实验内容：
考虑如下的在线DVD租赁问题。顾客缴纳一定数量的月费成为会 员，订购DVD租赁服务。会员对哪些DVD有兴趣，只要在线提交订 单，网站就会通过快递的方式尽可能满足要求。会员提交的订单包括多 张DVD，这些DVD是基于其偏爱程度排序的。网站会根据手头现有的 DVD数量和会员的订单进行分发。每个会员每个月租赁次数不得超过2 次，每次获得3张DVD。会员看完3张DVD之后，只需要将DVDa放进网 站提供的信封里寄回（邮费由网站承担），就可以继续下次租赁。请考 虑以下问题：
数学实验与数学建模作业(1)
3.1.假设船的模型为一个标准的矩形，即船头与船尾宽度一致，以便于分析问题。
3.2.假设船的最大长度为船外侧与河道的边沿刚好接触时的情况，方便求解。
四、符号பைடு நூலகம்明
L：船的长度。
X：船外侧与横轴的夹角为x。
五、模型的建立与求解
5.1.建立船长度L与x的关系：
L=10/Sinx+12/Cosx-5*Tanx-5/Tanx
Mathematica建立船渡直角弯模型
一、问题重述
有一艘宽度为5m的船，欲驶过某河道的直角湾，河道的宽度如下图所示，试问：要驶过直角湾,船的长度不能超过多少米?(精确到0.01m)
图一
二、问题分析
如图所示，由于船需要拐过的弯为一个直角，所以船的最大长度有一限制，即当船头船尾所成直线恰好与直角的内角相切时，设驳船长度为L，要使驳船能驶过直角湾，假定驳船外侧与河道的边沿刚好接触，则河道内侧的角点到驳船内侧的距离不能大于5m，否则无法通过.因而问题归结为求L的最大值。所以只需要建立船长L与船外侧与横轴的夹角x的关系，并利用mathematica求出L的最大值即可。
所以船欲驶过该直角湾，船的最大长度为21.04m。
5.2.画出函数的图形：
源代码：
f[x_]:=10/Sin[x]+12/Cos[x]-5*Tan[x]-5/Tan[x]
Plot[f[x],{x,0,Pi/2},PlotRange{-10,100}]
结果：
结果分析：由图可知，函数有唯一的极值点，大约在0.75附近。
5.3.画出直观的图像：
源代码：
Plot[ ,{x,0,Pi/2},PlotRange->{-60,60}]
结果：
数学建模与数学实验课后习题答案
P594•学校共1002名学生，237人住在A 宿舍，333人住在B 宿舍，432 人住在C 宿舍。
学生要组织一个10人的委员会，使用Q 值法分配各 宿舍的委员数。
解:设P 表示人数，N 表示要分配的总席位数。
i 表示各个宿舍（分别取 A,B,C ）， p i 表 示i 宿舍现有住宿人数， n i 表示i 宿舍分配到的委员席位。
首先，我们先按比例分配委员席位。
23710 A 宿舍为:n A ==2.365 1002 333"0 B 宿舍为：n B =3.323 1002 432X0 C 宿舍为：n C =4.3111002现已分完9人，剩1人用Q 值法分配。
经比较可得，最后一席位应分给 A 宿舍。
所以，总的席位分配应为： A 宿舍3个席位，B 宿舍3个席位，C 宿舍4个席位。
QA23722 3= 9361.5 Q B33323 4 = 9240.7 Q C4322 4 5=9331.2商人们怎样安全过河傻麴删舫紬削＜ I 11山名畝臥蹄峨颂禮训鋤嫌邂 韻靖甘讹岸讎鞍輯毗匍趾曲展 縣確牡GH 錚俩軸飙奸比臥鋪謎 smm 彌鯉械即第紘麵觎岸締熾 x^M 曲颁M 删牘HX …佛讪卜过樹蘇 卜允棘髒合 岡仇卅毘冋如;冋冋1卯;砰=口 於广歎煙船上觸人敦% V O J U;xMmm朗“…他1曲策D 咿川| thPl,2卜允隸策集合 刼為和啊母紳轉 多步贱 就匚叫=1入“山使曲并按 腿翻律由汩3』和騒側），模型求解 -穷举法〜编程上机 ■图解法S={(x ?jOI x=o, j-0,1,2,3;X =3? J =0,1,2,3； X =»*=1,2}J规格化方法，易于推广考虑4名商人各带一随从的情况状态$=（xy¥）~ 16个格点 允许状态〜U ）个。
点 , 允许决策〜移动1或2格； k 奇）左下移；&偶，右上移. 右，…，必I 给出安全渡河方案评注和思考［廿rfn片,rfl12 3xmm賤縣臓由上题可求：4个商人，4个随从安全过河的方案。
《数学建模与数学实验》上机实验报告
成都信息工程大学《数学建模与数学实验》上机实验报告专业信息与计算科学班级姓名学号实验日期成绩等级教师评阅日期[问题描述]下表给出了某一海域以码为单位的直角坐标Oxy 上一点（x，y）（水面一点）以英尺为单位的水深z，水深数据是在低潮时测得的，船的吃水深为5英尺，问在矩形区域（75，200）x （-50，150）里那些地方船要避免进入。
[模型]设水面一点的坐标为（x,y,z），用基点和插值函数在矩形区域（75，200）*（-50，150）内做二维插值、三次插值，然后在作出等高线图。
[求解方法]使用matlab求解：M文件：water.mx=[129 140 103.5 88 185.5 195 105.5 157.5 107.5 77 81 162 162 117.5];y=[7.5 141.5 23 147 22.5 137.5 85.5 -6.5 -81 3 56.5 -66.584 -33.5];z=[-4 -8 -6 -8 -6 -8 -8 -9 -9 -8 -8 -9 -4 -9];cx = 75:0.5:200;cy = -50:0.5:150;[cx,cy]=meshgrid(cx,cy);作出曲面图：代码如下：>> water>> cz=griddata(x,y,z,cx,cy,'cubic');>> meshz(cx,cy,cz)>> xlabel('X'),ylabel('Y'),zlabel('Z')>>作出等高线图：代码如下：>> water>> cz=griddata(x,y,z,cx,cy,'cubic');>> figure(2)>> contour(cx,cy,cz,[-5,-5],'r')>> hold on>> plot(x,y,'*')>> xlabel('X'),ylabel('Y')[结果]插值结果等值图：[结果分析及结论]根据等值图可看出：红色区域为危险区域，所以船只要避免进入。
《数学建模与实验》习题库a
Hyundai 车价$12,400 预付$500 月利率 6.5%直到 48 个月
你每个月为买车最多能付 475 美元。利用动力系统模型来决定你应该买哪家公司的车。
1.2 习题
1.从引进到 Tasmania 岛的新环境里的养群数量的增长得到下面的数据。
年
1814 1824 1834 1844 1854 1864
1810
7,240,000
1820
9,638,000
1830
12,866,000
1840
17,069,000
1850
23,192,000
1860
31,443,000
1870
38,558,000
1880
50,156,000
1890
62,948,000
1900
75,995,000
1910
91,972,000Leabharlann 数量 125275
830
1200 1750 1650
根据数据画图形，能看出某种趋势吗?画出 1814 年后数量变化对年份的图形。构建一个
能合理地近似描述你所观察到的变化的离散动力系统。
2.下列数据表示从 1790 年到 2000 年的美国人口数据
年份
人口
1790
3,929,000
1800
5,308,000
《数学建模与实验》习题库 a
感谢信息与计算科学 02 级的五位同学, 作为毕业设计英文翻译任务完成了此习题库的 构建工作, 他(她)们的工作分别为:
刘 静: 第 1, 4 章; 朱佳琦: 第 2, 3, 6 章; 李新颖: 第 5, 7 章; 朱晓强: 第 8, 9, 10 章; 甘永生: 第 11, 12 章.
数学建模与数学实验第五版课后答案4
数学建模与数学实验第五版课后答案4.41、27.下列各函数中，奇函数的是（）[单选题] *A. y=x^(-4)B. y=x^(-3)(正确答案)C .y=x^4D. y=x^(2/3)2、4．点（-3，-5）关于x 轴的对称点的坐标为（）[单选题] *A（-3，5）(正确答案)B（-3，-5）C（3，5）D（3，-5）3、1．(必修1P5B1改编)若集合P={x∈N|x≤2 022}，a=45，则( ) [单选题] * A．a∈PB．{a}∈PC．{a}?PD．a?P(正确答案)4、2．在+3，﹣4，﹣8，﹣，0，90中，分数共有（）[单选题] *A．1个B．2个C．3个(正确答案)D．4个5、13．在海上，一座灯塔位于一艘船的北偏东40°方向，那么这艘船位于灯塔（）[单选题] *A．南偏西50°方向B．南偏西40°方向(正确答案)C．北偏东50°方向D．北偏东40°方向6、5．已知集合A={x|x=3k＋1，k∈Z}，则下列表示不正确的是( ) [单选题] *A．－2∈AB．2 022?AC．3k2＋1?A(正确答案)D．－35∈A7、45．下列运算正确的是（）[单选题] *A．（5﹣m）（5+m）＝m2﹣25B．（1﹣3m）（1+3m）＝1﹣3m2C．（﹣4﹣3n）（﹣4+3n）＝﹣9n2+16(正确答案)D．（2ab﹣n）（2ab+n）＝4ab2﹣n28、10．如图是丁丁画的一张脸的示意图，如果用表示左眼，用表示右眼，那么嘴的位置可以表示成（）．[单选题] *A．（1，0）B（-1，0）(正确答案)C（-1，1）D（1，-1）9、6．方程x2=3x的根是（）[单选题] *A、x = 3B、x = 0C、x1 =-3, x2 =0D、x1 =3, x2 = 0(正确答案)10、13．在数轴上，下列四个数中离原点最近的数是（）[单选题] *A．﹣4(正确答案)B．3C．﹣2D．611、13．设x∈R，则“x3（x的立方）>8”是“|x|>2”的( ) [单选题] * A．充分而不必要条件(正确答案)B．必要而不充分条件C．充要条件D．既不充分也不必要条件12、35．若代数式x2﹣16x+k2是完全平方式，则k等于（）[单选题] * A．6B．64C．±64D．±8(正确答案)13、f（x）=-2x+5在x=1处的函数值为（）[单选题] *A、-3B、-4C、5D、3(正确答案)14、13.不等式x+3＞5的解集为（）[单选题] *A. x＞1B. x＞2(正确答案)C. x＞3D. x＞415、掷三枚硬币可出现种不同的结果（）[单选题] *A、6B、7C、8(正确答案)D、2716、6、已知点A的坐标是，如果且，那么点A在（）[单选题] *x轴上y轴上x轴上，但不能包括原点(正确答案)y轴上，但不能包括原点17、49．若（x+2）（x﹣3）＝7，（x+2）2+（x﹣3）2的值为（）[单选题] * A．11B．15C．39(正确答案)D．5318、1.计算| - 5 + 3|的结果是[单选题] *A. - 2B.2(正确答案)C. - 8D.819、4.一个数是25，另一个数比25的相反数大- 7，则这两个数的和为[单选题] *A.7B. - 7(正确答案)C.57D. - 5720、30.圆的方程+=4,则圆心到直线x-y-4=0的距离是（）[单选题] *A.√2(正确答案)B.√2/2C.2√2D.221、用角度制表示为（）[单选题] *30°(正确答案)60°120°-30°22、19．对于实数a、b、c,“a>b”是“ac2(c平方)>bc2(c平方) ; ”的（）[单选题] * A．充分不必要条件B．必要不充分条件(正确答案)C．充要条件D．既不充分也不必要条件23、10．下列四个数中，属于负数的是（）．[单选题] *A-3(正确答案)B 3C πD 024、6.有15张大小、形状及背面完全相同的卡片，卡片正面分别画有正三角形、正方形、圆，从这15张卡片中任意抽取一张正面的图形既是轴对称图形，又是中心对称图形的概率是1/3?，则正面画有正三角形的卡片张数为（）[单选题] *Ａ.3Ｂ.5Ｃ.10(正确答案)Ｄ.1525、若m·23＝2?，则m等于[单选题] *A. 2B. 4C. 6D. 8(正确答案)26、下列各式计算正确的是( ) [单选题] *A. (x3)3＝x?B. a?·a?＝a2?C. [(－x)3]3＝(－x)?(正确答案)D. －(a2)?＝a1?27、13．下列说法中，正确的为（）．[单选题] * A．一个数不是正数就是负数B. 0是最小的数C正数都比0大(正确答案)D. -a是负数28、函数式?的化简结果是（）[单选题] *A.sinα-cosαB.±（sinα-cosα）(正确答案)C.sinα·cosαD.cosα-sinα29、6．下列说法正确的是（）．[单选题] * A．不属于任何象限的点不在坐标轴上就在原点B．横坐标为负数的点在第二、三象限C．横坐标和纵坐标互换后就表示另一个点D．纵坐标为负数的点一定在x轴下方(正确答案)30、下列函数是奇函数的是（）[单选题] *A、f（x）=3x(正确答案)B、f（x）=4xC、f（x）= +2x-1D、f（x）=。
数学建模与数学实验课程设计题目与参考答案
数学建模与数学实验课程设计题目1、一元线性回归问题在某产品表明腐蚀刻线，下表是试验活得的腐蚀时间（x）与腐蚀深度（y）间的一组数据。
试研究两变量（x，y）之间的关系。
其中：(秒)()。
要求：1）画出散点图，并观察y与x的关系；=+，求出a与b的值；2）求y关于x的线性回归方程：y a bx3）对模型和回归系数进行检验；4）预测x=120时的y的置信水平为0.95的预测区间。
5）编程实现上述求解过程。
注：参考书目：1、《概率论与数理统计》，浙江大学编，高等教育出版社。
2、《数学实验》，萧树铁主编，高等教育出版社。
2、 多元线性回归问题根据下述某猪场25头育肥猪4个胴体性状的数据资料，试进行瘦肉量y 对眼肌面积（x1）画出散点图y 与x1，y 与x2，y 与x3并观察y 与x1，x2， x3的关系；2）求y 关于x1，x2， x3的线性回归方程：0112233y a a x a x a x =+++-----（1），求出0123,,,a a a a 的值；3）对上述回归模型和回归系数进行检验；4）再分别求y 关于单个变量x1，x2, x3的线性回归方程：10111y a a x =+----（2），20222y a a x =+-----（3）,30333y a a x =+--- --（4）求出ij a 的值；分别求y 关于两个变量x1，x2, x3的线性回归方程：10111122y a a x a x =++----(2’)，20211222y a a x a x =++---（3’）,30311322y a a x a x =++ --- --（4’）求出系数ij a 的值；并说明这六个回归方程对原来问题求解的优劣。
5）编程实现上述求解过程。
注：参考书目：1、《概率论与数理统计》，浙江大学编，高等教育出版社。
2、《数学实验》，萧树铁主编，高等教育出版社。
3、优化理论中的线性规划问题---生产安排。
《数学建模与数学实验》(第三版)6.5习题作业
1.电路问题一电路由三个电阻123R R R 、、并联，再与电阻4R 串联而成，记k R 上电流为k I ，电压为k V ，在下列情况下确定k R 使电路总功率最小(1,2,3,4)k =： 1）1234,6,8,2k I I I ===≤V ≤10; 2）1234,6,8,2k V V V I ===≤≤6;1）解：根据建立2;P I R U IR ==数学模型为：W=min 421k k k I R =∑123412346..82(1,...,4)kI I s t I I I I Ik I ⎧⎪=⎪=⎪⎪=⎨⎪=++⎪⎪=⎪⎩k k 10≤R ≤I用Lingo 求解：min =I1^2*R1+I1^2*R1+I2^2*R2 结果：+I3^2*R3+I4^2*R4;I1=4;I2=6;I3=8; I4=18; R1>1/2; R2>1/3; R3>1/4; R4>1/9; end2）解：根据建立2;P I R U IR ==数学模型为：W=min 421k k k I R =∑ 4123112233R =4/I ;..R =6/I ;R =8/I ;2I I I I s t I =++⎧⎪⎪⎪⎨⎪⎪⎪⎩k ≤≤6(k =1,...,4);用Lingo 求解：min =I1^2*R1+I2^2*R2+I3^2*R3 结果：+I4^2*R4;I4=I1+I2+I3;I1<6; I2<6;I3<6;I4<6; 《数学建模与数学实验》（第三版）6.5习题作业专业 班级 姓名 学号12340.50000.33330.25000.1111R R R R =⎧⎪=⎪⎨=⎪⎪=⎩ , 1234 4.00006.00008.000018.0000I I I I =⎧⎪=⎪⎨=⎪⎪=⎩ 80P = 112233440.5835976E+08 0.6854038E-07 0.1586609E+08 0.3781429E-06 1.3333 6.000000 0.4752196E+27 6.000000R I R I R I R I ==⎧⎧⎪⎪==⎪⎪⎨⎨==⎪⎪⎪⎪==⎩⎩0.1710790E+29P =R1=4/I1; R2=6/I2; R3=8/I3; end3.（设计最优化问题）要设计和发射一个带有X 射线望远镜和其他科学仪器的气球。
数学实验与数学建模MATLAB实验报告78
数学实验与数学建模MATLAB实验报告78数学实验与数学建模实验报告学院：信息科学与⼯程学院专业班级：姓名：学号：习题七1.求下列微分⽅程的通解（1）x y x y dx dy -+=（2）yxx y y +=cos ' （3）（xcosy+sin2y ）y`=1 (4)x ey y y x2cos 3=-'+''(5) x y e y y x 2cos 3'''=-+解：(1)dsolve('Dy=(y+x)/(y-x)','x')(2)dsolve('Dy=cos(y/x)+x/y','x')(3)dsolve('(x*cos(y)+sin(2*y))*Dy=1','x')(4)dsolve('D2y+3*Dy-y=exp(x)*cos(2*x)','x')(5)dsolve('D2y+3*Dy-y=exp(x*cos(2*x))','x')ans=exp(1/2*(-3+13^(1/2))*x)*C2+exp(-1/2*(3+13^(1/2))*x)*C1-1/13*13^(1/2)*(-Int(exp (1/2*x*(3-13^(1/2)+2*cos(2*x))),x)*exp(x*13^(1/2))+Int(exp(1/2*x*(3+13^(1/2)+2*cos(2*x))),x))*exp(-1/2*(3+13^(1/2))*x) 2．求下列初值问题的解(1)==-++-+=10)2(212222y x y y x x dx dy xy xy (2)????===++==V dt dx x x a t t x dt dx n dt x d 000222,02解：(1) dsolve('x^2+2*x*y-y^2+(y^2+2*x*y-x^2)*Dy=0','y(1)=1','x')(2) dsolve('D2x+2*n*Dx+a^2*x=0','x(0)=x0','Dx(0)=V0','t')ans =1/2*(n*x0+(n^2-a^2)^(1/2)*x0+V0)/(n^2-a^2)^(1/2)*exp((-n+(n^2-a^2)^(1/2))*t )-1/2*(n*x0-(n^2-a^2)^(1/2)*x0+V0)/(n^2-a^2)^(1/2)*exp((-n-(n^2-a^2)^(1/2))*t)3．求微分⽅程组=--=++t te y x dtdy e y x dtdx 235的通解.解：[x,y]=dsolve('Dx+5*x+y=exp(t)','Dy-x-3*y=exp(2*t)','t')x =-4*exp((-1+15^(1/2))*t)*C2+exp((-1+15^(1/2))*t)*C2*15^(1/2)-4*exp(-(1+15^(1/2))*t)*C1-exp(-(1+15^(1/2))*t)*C1*15^(1/2)+2/11*exp(t)+1/6*exp(2*t) y =exp((-1+15^(1/2))*t)*C2+exp(-(1+15^(1/2))*t)*C1-1/11*exp(t)-7/6*exp(2*t) 4．求下列初值问题的解(1)⽅程组+=+=11x dtdy y dt dx满⾜=-=0)0(2)0(y x 的特解。
《数学建模与实验》实验报告
1、输入数据： >> x=[20 25 30 35 40 45 50 55 60 65]'; 实 X=[ones(10,1) x]; Y=[13.2 15.1 16.4 17.1 17.9 18.7 19.6 21.2 22.5 24.3]'; [b,bint,r,rint,stats]=regress(Y,X) 验 b,bint,stats 得出结果： >> b = 步 9.1212 0.2230 bint = 骤 8.0211 0.1985 Stats = 0.9821 439.8311 0.0000 10.2214 0.2476
2 ˆ t 得回归模型： s 9.1329 65.8896 489.2946t
《数学建模与实验》实验报告
实验名称 班级 实验目的 数学 09-1 实验六 MATLAB 回归分析 姓名 学号 25 号
1.熟悉 MATLAB 基本命令； 2.掌握回归分析的方法。 1、考察温度 x 对产量 y 的影响，测得下列 10 组数据：
温度（℃） 20 产量（kg） 13.2 25 15.1 30 16.4 35 17.1 40 17.9 45 18.7 50 19.6 55 21.2 60 22.5 65 24.3
从残差图可以看出，所有数据的残差离零点均较近，且残差 的置信区间均包含零点，这说明回归模型 y=9.1212+0.2230x 能较好的符合原始数据. 预测及作图： >> z=b(1)+b(2)*x plot(x,Y,'k+',x,z,'r') 当 x=42℃时，产量估值 z=18.4885 . 预测区间：[16.3581，20.6206]（置信度 95%） 2、作二次多项式回归： >>xi =0:2:20; yi =[0.6 2.0 4.4 7.5 11.8 17.1 23.3 31.2 39.6 49.7 61.7]; [p,S]=polyfit(xi, yi,2) 得到结果： p= 0.1403 S= R: [3x3 double] df: 8 normr: 1.1097
《数学建模与数学实验》课程作业
《数学建模与数学实验》实验报告学院班级姓名学号二零一二年六月《数学建模与数学实验》课程作业一、简要说明MATLAB有那几个主要的界面？说明其作用是什么？1．与Windows的窗口界面类似，有File、Edit、Option、Windows、HelpFile菜单项：实现有关文件的操作。
Edit菜单项：用于命令窗口的编辑操作。
View菜单项：用于设置MATLAB集成环境的显示方式。
Web菜单项：用于设置MATLAB的Web操作。
Window菜单项：主窗口菜单栏上的Window菜单，只包含一个子菜单Close all，用于关闭所有打开的编辑器窗口，包括M-file、Figure、Model和GUI窗口。
Help菜单项：为MATLAB的学习提供在线和系统自带的帮助信息。
2.窗口（１）命令窗口。
用于输入命令并显示除图形以外的所有执行结果。
（２）工作空间窗口。
用于存储各种变量和结果的内存空间，显示工作空间中所有变量的名称、大小、字节数和变量类型说明，可对变量进行观察、编辑、保存和删除。
（３）当前目录窗口和搜索路径。
可以显示或改变当前目录，还可以显示当前目录下的文件并提供搜索功能。
（４）命令历史记录窗口。
自动保留自安装起所有用过的命令的历史记录，并且还标明了使用时间，从而方便用户查询。
（５）启动平台窗口。
帮助用户方便地打开和调用MATLAB的各种程序、函数和帮助文件。
二、简要说明你对数学建模的看法。
应用数学知识解决实际问题，并了解到相关数学软件的使用三、输入下面的矩阵A、B并完成相应的运算；5200210000830052A ⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎣⎦1000120021301214B ⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎣⎦1．求出矩阵A 的逆矩阵、矩阵A 的秩、矩阵A 所对应的行列式的值、A^2; 解：命令>> A=[5 2 0 0;2 1 0 0;0 0 8 3;0 0 5 2] A =5 2 0 0 2 1 0 0 0 0 8 3 0 0 5 2 矩阵A 的逆矩阵 >> inv(A) ans =1.0000 -2.0000 0 0 -2.0000 5.0000 0 0 0 0 2.0000 -3.0000 0 0 -5.0000 8.0000 矩阵A 的秩 >> rank(A) ans =4矩阵A 所对应的行列式的值 >> det(A)ans =1 A^2 >> A^2 ans =29 12 0 0 12 5 0 0 0 0 79 30 0 0 50 192.求出矩阵A 的伴随矩阵、矩阵A 的特征值及特征向量、矩阵A 对应的上三角矩阵和下三角矩阵及将矩阵、将矩阵A 化为最简的阶梯型矩阵；解：矩阵A 的伴随矩阵>> det(A)*inv(A)ans =1.0000 -2.0000 0 0-2.0000 5.0000 0 00 0 2.0000 -3.00000 0 -5.0000 8.0000 矩阵A的特征值及特征向量>> [V,D]=eig(A,'nobalance')V =1.0000 -0.4142 0 00.4142 1.0000 0 00 0 1.0000 -0.37980 0 0.6330 1.0000D =5.8284 0 0 00 0.1716 0 00 0 9.8990 00 0 0 0.1010 矩阵A对应的上三角矩阵>> triu(A)ans =5 2 0 00 1 0 00 0 8 30 0 0 2矩阵A对应的下三角矩阵>> tril(A)ans =5 0 0 02 1 0 00 0 8 00 0 5 2最简阶梯型矩阵>> rref(A)ans =1 0 0 00 1 0 00 0 1 00 0 0 12.完成下列矩阵的运算A*B、A/B、A\B、A.*B、A./B；解: >>A=[5 2 0 0;2 1 0 0;0 0 8 3;0 0 5 2]>>B=[1 0 0 0;1 2 0 0;2 1 3 0;1 2 1 4]>> A*Bans =7 4 0 03 2 0 019 14 27 1212 9 17 8>> A/Bans =4.0000 1.0000 0 01.5000 0.5000 0 0-3.6250 -1.9583 2.4167 0.7500-2.2500 -1.2500 1.5000 0.5000 >> A\Bans =-1.0000 -4.0000 0 03.0000 10.0000 0 01.0000 -4.0000 3.0000 -12.0000-2.0000 11.0000 -7.0000 32.0000 >> A.*Bans =5 0 0 02 2 0 00 0 24 00 0 5 8>> A./BWarning: Divide by zero.ans =5.0000 Inf NaN NaN2.0000 0.5000 NaN NaN0 0 2.6667 Inf0 0 5.0000 0.5000四、解下面的线性方程组；（1）123412423412342583692254760 x x x xx x xx x xx x x x⎧+-+=⎪--=⎪⎨-+=-⎪⎪+-+=⎩>> A=[2 1 -5 1;1 -3 0 -6;0 2 -1 2;1 4 -7 6]>> b=[8 9 -5 0]>> rank(A)ans =4>> rank([A,b'])ans =4运行结果：r(A)=r(A|b)=n，则线性方程组存在唯一解>> A\b'ans =3.0000-4.0000-1.00001.0000（2）123123123231 2252 353 x x xx x xx x x⎧++=⎪++=⎨⎪++=⎩命令>> A=[1 2 3;2 2 5;3 5 1]>> b=[1 2 3]>> rank(A)ans =3>> rank(A,b')ans =2运行结果：r(A)≠r(A|b)，则线性方程组无解; 最小二乘意义上的近似解>> A\b'ans1五、解决下列高等数学中的问题；1.求出下列极限的值； （1）设1/1()1xf x e-=+，求当1,0,0,x x x x +- 时函数的极限；命令>> syms x>> f=1/(1+exp(-1/x)) f =1/(1+exp(-1/x)) >> limit(f,x,1) ans =1/(1+exp(-1))>> limit(f,x,0,'right') ans = 1>> limit(f,x,0,'left') ans = 0>> limit(f,x,inf) ans = 1/22.求出下列函数的导数值； （1）求出函数22cos x x y ee--=的一阶导数；命令>> syms x>> y=exp(-x^2)*cos(exp(-x^2)) y =exp(-x^2)*cos(exp(-x^2)) >> diff(y,x) ans =-2*x*exp(-x^2)*cos(exp(-x^2))+2*exp(-x^2)^2*sin(exp(-x^2))*x （2）求出函数(23)x y x e =+的2阶及4阶导数； 命令>> syms x>> y=(2*x+3)*exp(x (2*x+3)*exp(x) >> diff(y,2) ans =4*exp(x)+(2*x+3)*exp(x) >> diff(y,4)ans =8*exp(x)+(2*x+3)*exp(x)（3）求出函数22()2ln[()]x yz e x y+=++的2422,,,z z z zx y x y x y抖抖抖抖抖偏导数导数；命令>> syms x y z>> z=log(exp(2*(x+y^2))+(x^2+y))z =log(exp(2*x+2*y^2)+x^2+y)>> diff(z,x)ans =(2*exp(2*x+2*y^2)+2*x)/(exp(2*x+2*y^2)+x^2+y)>> diff(z,y)ans =(4*y*exp(2*x+2*y^2)+1)/(exp(2*x+2*y^2)+x^2+y)>> diff(diff(z,x),y)ans =8*y*exp(2*x+2*y^2)/(exp(2*x+2*y^2)+x^2+y)-(2*exp(2*x+2*y^2)+2*x)/(exp(2*x+2* y^2)+x^2+y)^2*(4*y*exp(2*x+2*y^2)+1)>> diff(diff(z,x,2),y,2)ans =16*exp(2*x+2*y^2)/(exp(2*x+2*y^2)+x^2+y)+64*y^2*exp(2*x+2*y^2)/(exp(2*x+2*y ^2)+x^2+y)-32*y*exp(2*x+2*y^2)/(exp(2*x+2*y^2)+x^2+y)^2*(4*y*exp(2*x+2*y^2)+1)+ 2*(4*exp(2*x+2*y^2)+2)/(exp(2*x+2*y^2)+x^2+y)^3*(4*y*exp(2*x+2*y^2)+1)^2-(4*exp( 2*x+2*y^2)+2)/(exp(2*x+2*y^2)+x^2+y)^2*(4*exp(2*x+2*y^2)+16*y^2*exp(2*x+2*y^2))-128*y^2*exp(2*x+2*y^2)^2/(exp(2*x+2*y^2)+x^2+y)^2+64*(2*exp(2*x+2*y^2)+2*x)/(exp (2*x+2*y^2)+x^2+y)^3*y*exp(2*x+2*y^2)*(4*y*exp(2*x+2*y^2)+1)-16*(2*exp(2*x+2*y^ 2)+2*x)/(exp(2*x+2*y^2)+x^2+y)^2*exp(2*x+2*y^2)-64*(2*exp(2*x+2*y^2)+2*x)/(exp(2* x+2*y^2)+x^2+y)^2*y^2*exp(2*x+2*y^2)-6*(2*exp(2*x+2*y^2)+2*x)^2/(exp(2*x+2*y^2) +x^2+y)^4*(4*y*exp(2*x+2*y^2)+1)^2+2*(2*exp(2*x+2*y^2)+2*x)^2/(exp(2*x+2*y^2)+x ^2+y)^3*(4*exp(2*x+2*y^2)+16*y^2*exp(2*x+2*y^2))3．求出下列积分的值；（1）ln tansin cosxdxx x ò命令>> syms x>> log(tan(x))/(sin(x)*cos(x)) ans =log(tan(x))/sin(x)/cos(x)>> int(f,x) ans =-dilog(1-i*exp(i*x))+dilog(exp(i*x)+1)+log(exp(i*x)+1)*log(1/2-1/2*exp(i*x)+1/2*i*(e xp(i*x)+1))+log(exp(i*x)-1)*log(i*(1-exp(i*x)^2)/(exp(i*x)^2+1))-log(exp(i*x))*log(i*(1-ex p(i*x)^2)/(exp(i*x)^2+1))-1/2*log(exp(i*x)+1)^2-log(2)*log(1/2*exp(i*x)-1/2)+log(exp(i*x)+1)*log(1/2-1/2*exp(i*x)-1/2*i*(exp(i*x)+1))+dilog(1/2-1/2*exp(i*x)-1/2*i*(exp(i*x)+1))+d ilog(1/2-1/2*exp(i*x)+1/2*i*(exp(i*x)+1))-log(exp(i*x))*log(1-i*exp(i*x))-log(exp(i*x))*log (1+i*exp(i*x))+log(exp(i*x))*log(exp(i*x)+1)-dilog(exp(i*x))-dilog(1+i*exp(i*x))+log(exp(i *x)+1)*log(i*(1-exp(i*x)^2)/(exp(i*x)^2+1))+log(exp(i*x)-1)*log(1/2+1/2*exp(i*x)-1/2*i*(e xp(i*x)-1))+log(exp(i*x)-1)*log(1/2+1/2*exp(i*x)+1/2*i*(exp(i*x)-1))-log(exp(i*x)-1)*log(1/2*exp(i*x)+1/2)+dilog(1/2+1/2*exp(i*x)+1/2*i*(exp(i*x)-1))+dilog(1/2+1/2*exp(i*x)-1/2*i *(exp(i*x)-1))-1/2*log(exp(i*x)-1)^2（2）83xdxò命令>> syms x>> f=x/(1+x)^(1/2) f =x/(1+x)^(1/2) >> int(f,x,3,8) ans = 32/3（3）计算二重积分22121x xx dydx y蝌>> syms x y 命令>> f=x^2/y^2 f =x^2/y^2>> int(int(f,y,1/x,x),x,1,2)Warning: Explicit integral could not be found. > In sym.int at 58 ans = 9/4六、绘制下列函数的图形（1）1sin(),[0.1,0.1]y x x=?命令>> x=-0.1:0.001:0.1 >> y=sin(1./x)Warning: Divide by zero. >>plot(x,y)-0.1-0.08-0.06-0.04-0.020.020.040.060.080.1-1-0.8-0.6-0.4-0.200.20.40.60.81（2）sin(cos ),[0,3]y x x x x p =++01234567891024681012命令>> x=0:pi/100:3*pi >> y=x+sin(x+cos(x)) >>plot(x,y)（3）20y x xy e +-= 命令>> syms x y>> f=x^2+x*y-exp(y) >> ezplot(f)xyx 2+x y-exp(y) = 0-6-4-20246-6-4-2246（4） 22ln(1)z x y =+- >> x=-1:0.1:1 >> y=-1:0.1:1>> [x,y]=meshgrid(x,y) >> z=log(x^2+y^2-1) >> mesh(x,y,z)-11七、谈谈你对数学建模和数学实验选修课程的看法和改进意见。
数学建模实习作业
碎纸片的拼接复原摘要图像碎片自动拼接复原是需要借助计算机把大量碎片重新拼接复原成初始图像的完整模型，这一研究在考古、刑侦犯罪、古生物学、医学图像分析、遥感图像处理以及壁画保存复原等方面具有广泛、实际的应用。
本文主要解决碎纸机破碎文档的自动拼接复原问题。
我们利用图像数字化技术，借助Matlab软件将图像转化为矩阵。
通过建立数学模型，运用矩阵论、自定义相似度方法、遗传算法等方法，对数据进行处理，实现对图像碎片自动拼接，从而将所给碎片拼接复原为完整图像。
我们首先把碎片图形进行二值化处理，根据所给纵切黑白碎片边缘的像素关系(相邻两张碎片，一张碎片矩阵右边的像素与另一张碎片左边的像素相同 )，我们采和自定义相似度算法，利用附件求出碎片间的相似度，然后根据所需要满足的条件即相似度最大原则，建立了纵切碎片拼接模型一及其算法，运用Matlab编程实现该模型，并得到碎片复原结果（见附录1）。
关键词：碎片拼接矩阵论图形二值化相似度模型一、问题重述1.1背景：破碎文件的拼接和复原对于司法物证复原、历史文献再现和军事情报获取等方面都有极其重要的作用。
于是碎纸片的拼接复原技术便成为图像处理与模式识别领域中的一个崭新典型的应用。
图像配准是图像拼接复原的基础，而且图像配准算法的计算量一般非常大，因此图像拼接复原技术的发展很大程度上取决于图像配准技术的创新。
本文将通过图像提取技术获取一组碎纸片的形状、颜色、文字等信息，然后利用计算机进行相应的处理从而实现对这些碎纸片的自动拼接复原。
1.2重述：该题研究的是如何对碎纸片进行拼接复原。
传统上，拼接复原工作需由人工完成，准确率较高，但是效率低。
随着计算机技术的发展，当碎纸片数量巨大的时候，人们试图开发碎纸片的自动拼接技术，以提高拼接复原的效率。
对于给定的来自同一页印刷文字文件的碎纸机破碎纸片（仅纵切），建立碎纸片拼接复原模型和算法，并针对附件给出的中文文件的碎片数据进行拼接复原。
如果复原过程需要人工干预，写出干预方式及干预的时间节点。
数学实验与数学建模实验报告
数学实验与数学建模实验报告学院：湘雅医学院专业班级：不告诉你姓名：郝甍学号：完成时间：2012年6月22日承诺书本人承诺所呈交的数学实验与数学建模作业都是本人通过学习自行进行编程独立完成，所有结果都通过上机验证，无转载或抄袭他人，也未经他人转载或抄袭。
若承诺不实，本人愿意承担一切责任。
承诺人：2012年 6 月20 日注意事项如下：1、2012年6月22日(第十八周星期五)之前，将电子文档发送到邮箱：xuanyunqin@（word文档命名：姓名＋学号＋数学实验作业）2、2012年6月22日(第十八周星期五)，将实验报告电子打印稿交到物理楼数学实验室办公室，过时不再受理。
谢谢同学们合作!!!数学实验学习体会（每个人必须要写1500字以上，占总成绩的20%）通过几周的学习，我对MATLAB数学实验与建模有了更生的认识。
我感受到MATLAB强大的运算能力和实用性。
但最深刻的感受就是：要不断地用它。
MATLAB是个好工具，但如果不用他来解决问题，只知道一点语法，那是连皮毛都没有学到的。
还有就是程序设计，对于程序的运行效率非常有帮助。
有的时候，编出来的程序能够运行，但是耗时太长，程序虽然没有错，但是不是和实际应用。
这就需要对程序的结构和算法问题进行改进，要时刻思考多动脑，找到十一的解题途径。
还有就是学习MATLAB要多动手，找一个习题实际操作一下或者找一个实际的程序来动手编一下，才能更好地对MATLAB有所了解，进一步巩固知识。
要在编程的过程中学习，程序需要什么只是再去补充，变成是一点一点积累的，需做一些随手笔记，我就是在这个时候有所懈怠才发懵的。
当然，除了要去用它以外，辅导书可以很大程度上提高我们的知识与技能，通过模仿别人编写的程序，可以大大加快我们掌握它的进度，并且学到一些课堂中所没有的知识.实验一图形的画法1. 做出下列函数的图像：（1）)2sin()(22--=xxxxy，22≤≤-x（分别用plot、fplot）（2）22/9/251x y+=（用参数方程）(3) 在同一图形窗口中，画出四幅不同图形（用subplot命令）：1cos()y x=，2sin(/2)y x pi=-，23cos()y x x pi=-，sin()4xy e=（]2,0[π∈x）（1）>> x=-2:0.001:2;>> y=x.^2.*sin(x.^2-x-2);>> plot(x,y)>> fplot('x.^2.*sin(x.^2-x-2)',[-2,2])（2）>> t=0:0.001:2*pi;>> x=9*cos(t);>> y=25*sin(t);>> plot(x,y)（3）>> x=0:0.01:2*pi;>> figure(1);>> subplot(2,2,1);>> y1=cos(x);>> plot(y1);>> subplot(2,2,2);>> y2=sin(x-pi/2);>> plot(y2);>> subplot(2,2,3);>> y3=(x.^2).*cos(x-pi);>> plot(y3);>> subplot(2,2,4);>> y4=exp(sin(x));>> plot(y4)2作出极坐标方程为)cos1(2tr-=的曲线的图形.>> t=linspace(-2*pi,2*pi,1000);>> r=2*(1-cos(t));>> plot(r) 图1.1.1 图1.1.2图1.1.33 作出极坐标方程为10/t e r =的对数螺线的图形.>> t=0:0.01:2*pi; >> polar(t,exp(t/10))4 绘制螺旋线⎪⎩⎪⎨⎧===t z t y t x ,sin 4,cos 4在区间［0，π4］上的图形.在上实验中，显示坐标轴名称。