作业要求：
1）编写程序并上机实现，提交作业文档，包括打印稿（不含源程序）和电子稿（包含源程序），以班为单位交，作业提交截止时间6月24日。
2）作业文档内容：问题描述、问题求解算法（方案）、matlab程序、结果分析、本课程学习体会、列出主要的参考文献。打印稿不要求matlab程序，但电子稿要包含matlab程序。
3）作业文档字数不限，但要求写实，写出自己的理解、收获和体会，有话则长，无话则短。不要抄袭复制，可以参考网上、文献资料的内容，但要理解，要变成自己的语言，按自己的思路组织内容。
4）从给出的问题中至少选择一题(多做不限，但必须独立完成，严禁抄袭)。
5）大作业占过程考核的20%，从完成情况、工作量、作业文档方面评分。
第一类：绘制图形。（b级）
问题一：斐波那契（fibonacci）螺旋线，也称**螺旋线（golden spiral），是根据斐波那契数列画出来的螺旋曲线，自然界中存在许多斐波那契螺旋线的图案，是自然界最完美的经典**比例。斐波那契螺旋线，以斐波那契数为边的正方形拼成的长方形，然后在正方形里面画一个90度的扇形，连起来的弧线就是斐波那契螺旋线，如图所示。
问题二：绘制谢尔宾斯基三角形（sierpinski********）是一种分形，由波兰数学家谢尔宾斯基在2024年提出，它是一种典型的自相似集。其生成过程为：
取一个实心的三角形（通常使用等边三角形），沿三边中点的连线，将它分成四个小三角形，然后去掉中间的那一个小三角形。接下来对其余三个小三角形重复上述操作，如图所示。
问题三：其他分形曲线或图形。分形曲线还有很多，教材介绍了科赫曲线，其他还有皮亚诺曲线、分形树、康托（g.
cantor）三分集、julia集、曼德布罗集合（mandelbrot set），等等。这方面的资料很多（如请分析构图原理并用matlab实现。
问题四：模拟掷骰子游戏：掷1000次骰子，统计骰子各个点出现的次数，将结果以下表的形式显示，并绘制出直方图。
问题五：利用matlab软件绘制一朵鲜花，实现一定的**效果。
提示：二维/三维绘图，对花瓣、花蕊、叶片、花杆等的形状和颜色进行详细设置。
第二类：插值与拟合。（b级）
问题一：有人对汽车进行了一次实验，具体过程是，在行驶过程中先加速，然后再保持匀速行驶一段时间，接着再加速，然后再保持匀速，如此交替。注意，整个实验过程中从未减速。
在一组时间点上测得汽车的速度如表所示。
1）分别使用最近点插值、线性插值、三次埃尔米特插值和三次样条插值进行计算[0,110]时间段50个时间点的速度。
2）绘制插值图形并标注样本点。
问题二：估算矩形平板各个位置的温度。已知平板长为5m，宽为3m，平板上3×5栅格点上的温度值为44，25，20，24，30；42，21，20，23，38；25，23，19，27，40。
1）分别使用最近点插值、线性插值和三次样条插值进行计算。
2）用杆图标注样本点。
3）绘制平板温度分布图。
问题三：自行车道的设计。对9条道路上的自行车道宽度以及自行车与过往机动车之间的平均距离进行测量，数据如表所示。
1）对数据进行线性拟合。
2）绘制拟合曲线和样本点。
3）如果自行车与过往机动车之间安全距离的最小距离是1.8m，试计算相应的自行车道宽度的最小值。
问题四：在水资源工程学中，水库的大小与为了蓄水而拦截的河道中的水流速度密切相关。对于某些河流来说，这种长时间的历史水流记录很难获得。
然而通常容易得到过去若干年间关于降水量的气象资料。鉴于此，推导出流速与降水量之间的关系式往往特别有用。只要获得那些年份的降水量数据，就可以利用这个关系式计算出水流速度。
下表是在被水库拦截的某河道中测得的数据。
1）对数据进行线性拟合。
2）绘制拟合曲线和样本点。
3）如果某年的降水量是120cm，利用拟合直线估算当年的水流速度。
4）若流域面积为1100km2，估计在其他过程中，如蒸发、深层地下水渗透和消耗用途，损失的降水量占总体降水量的比例。
问题五：假设有已知实测数据如下表所示：
假设已知该数据可能满足的原型函数为，试求出满足数据的最小二乘解a,b,c,d的值。
提示：曲线拟合并绘图分析。
第三类：定积分问题。（b级）
问题一：地球密度随着离中心（r=0）距离的变化而变化，不同半径处的密度如表所示，试估算地球质量。
问题二：河道平均流量q（m3/s）可使用速度和深度的乘积的积分来计算（河道横截面不规则），公式如下。
其中v(x)是离岸x（m）距离处的水速（m/s）,h(x)是离岸x距离处的水深（m）。根据收集到的河道离岸不同距离处的水速v和水深h（如表所示），估计流量。
第四类：线性方程组求解。（b级）
问题一：多项式插值指的是采用唯一的n-1次多项式对n个数据点进行拟合。该多项式的一般形式为：
p(x)=p1xn-1+p2xn-2+…+pn-1x+pn
确定这些系数的一种直接方法是，建立n个线性代数方程，然后求解。已知一个四次多项式通过5个点，如表所示。
1）建立线性方程组，并求解得到多项式的系数。
2）计算该线性方程组系数矩阵的条件数，并进行解释。
3）绘制多项式曲线并求其零点。
问题二：如图所示，5个反应器通过导管连接在一起。每根导管中化学物的传输率等于流速（q，单位是m3/s）乘以化学物浓度（c，单位是mg/m3）。
若系统达到稳定状态，流入和流出每个反应器的质量相等。例如，对于第一个反应器来说，质量守恒可表示为：
q01c01+q31c3=q15c1+q12c1
1）使用lu分解计算平衡方程系数矩阵的逆矩阵。
2）求各反应器中化学物的稳态浓度。
问题三：静定桁架受力分析。
1）如图所示，求力和反作用力。
2）求受力平衡方程系数矩阵的逆矩阵，对于逆矩阵第二行中的零，作何解释。
3）将节点1的力改为方向向上，计算这种改变对h2和v2的影响。
4）将节点1的力撤销，而在节点1和2处施加1500n的水平外力，求节点3处垂直反作用力（v3）。
第五类：一元方程求解。（b级）
问题一：在热力学中，下列多项式将干燥空气的零压力比热cp（单位为kj/（kgk））与温度（单位为k）关联起来了：
cp=0.99403+1.671×10-4t+9.7215×10-8t2-9.5838×10-11t3+1.9520×10-14t4
1）绘制在t=0~1200k范围内，cp随温度变化的曲线。
2）求对应于1.1kj/（kgk）比热的温度。
问题二：在化学工程中，将水蒸汽（h2o）加热到足够高的温度，使得大部分水发生分解或分离而形成氧气（o2）和氢气（h2）：
h2oh2+o2
如果假定其中只存在这一种化学反应，那么已经发生分解的h2o所占比列x可以表示为：
其中k为该反应的平衡系数，pt为混合物的总压强。如果pt=4，且k=0.05，那么求满足该式子的x值。
第六类：最优化问题。（b级）
问题一：最大利润问题。某公司生产两种产品的产量分别为x，ykg，其相应的成本满足以下函数：
c(x,y)=x2+2xy+2y2+2000
已知产品x的**为200元/kg，产品y的**为300元/kg，并假定两种产品全部售完，试求使公司获得最大利润产品产量以及公司的最大利益润。
问题二：作用在螺旋桨上的总阻力可以通过下式估计：
摩擦力升力。
其中，d=阻力， =飞行高度与海平面之间的大气密度比(ratioofairdensity)，w=重量，v=速度。
如图所示，当速度增加时，对阻力的两个部分受到的影响是不同的。摩擦阻力随速度的增加而增加，但由升力引起的阻力却随速度的增加而下降。二者的结合导致一个最小的阻力。
1）如果=0.6、w=16000，求最小阻力及阻力最小时的速度值。
2）进行敏感性分析以确定当w为12000~20000的过程中，最优值是如何变化的，取=0.6。
螺旋桨上阻力与速度的关系图。
问题三：如图所示，一个梯子通过支撑角分别与两个面接触，梯子的最大可能长度可以通过计算下面函数取值最小时的值而确定。
对于==2m的情况，绘制l随变化的图形，其中的取值范围为45°~135°。
通过一个墙角连接两个墙面的梯子。
问题四：对于一架稳定水平航行的喷气机，推力与阻力平衡，升力与重力平衡(如图所示)。在这种情况下，当阻力与速度的比例最小时，会出现最佳巡航速度。阻力可以用下式计算：
其中是零升力时的阻力系数，是升力系数，ar是展弦比。在稳定水平飞行的。
情况下，升力系数可以用下式计算：
其中w是喷气机重量（n），是空气密度（kg/）,是速度(m/s),a是机翼平面面积()，然后阻力可以用下式计算：
在稳定水平飞行中，喷气机受到的四个主要力。
使用这些公式，确定在海平面上10千米飞行的670kn喷气机的最佳稳定巡航速度。在计算中应用以下参数：a=150，ar=6.5， =0.018， kg/。
matlab大作业
matlab语言及应用。周一晚作业廖育州07010559 第一组程序。绘制选用潍柴发动机wd615.50时汽车行驶性能曲线图，ft ua图clear，clc 清除内存变量，清屏 n为发动机的转速。n 1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,...
MATLAB大作业
实验七图像的形态 算。一 实验目的。1.掌握图像的奇异值分解与合成。2.能够利用模板操作实现图像的腐蚀运算3.能够利用模板操作实现图像的膨胀运算4.掌握图像的开 闭运算。二 实验内容。1.读下面的程序，分析该程序要完成什么样的功能。a 1 2 1 2 3 2 1 2 1 u d v svd a nx...
MATLAB大作业
1.设计一个程序，要求用户用键盘输入10次，输入内容可以是数字 英文 汉字 或者其他文字或者符号。如果输入的是数字，则把存入一个数字数组当中 如果输入的是英文的字符，则存入一个英文字符数组当中 如果是其他字符，则存入一个字符串数组中。当10次内容都输入完毕后，先输出英文数组中的内容，然后输出字符串数...