北航MATLAB第一次作业
1.变量名如何声明，变量名必须遵守什么规则、是否区分大小写。
（1）变量声明：不需要声明
（2）变量名须遵守的规则：○1变量名必须以字母开头，且只能由字母、数字或下划线3类符号组成，不能含有空格和标点符号（如（），。%）等；○2变量名区分字母的大小写；
○3变量名必须是不含空格的单个词；○4关键字（如if、while）不能作为变量名。
（3）需要区分大小写。
2.试说明分号、逗号、冒号的用法。
（1）分号：在语句最后加分号可使命令或语句的运行结果不显示在屏幕上，矩阵的行与行之间用分号；
（2）逗号：一条命令行可以输入若干条命令，各命令之间以逗号或者分号分隔
（3）冒号：利用冒号生成向量。用于生成一维数值数组，表示一维数组的全部元素或多维数组的某一维的全部元素。
3.说明函数ones（）、zeros（）、eye（）、linspace（）说明它的用法。
（1）ones（）函数：全部元素都为1的常数矩阵；
（2）zeros（）函数：全部元素都为0的矩阵；
（3）eye（）函数：单位矩阵；
（4）linspace（）函数：如a=linspace(n1,n2,n3),表示在线性空间上，行矢量的值从n1到n2，数据个数为n。
作业二：
1.使用logspace（）创建1~4π的有10个元素的行向量。
>>logspace(0,log10(4*pi),10)
ans =
Columns 1 through 6
1.0000 1.3248 1.7550
2.3249
3.0799
4.0801
Columns 7 through 10
5.4051 7.1605 9.4858 12.5664
2.列举出两种创建向量A=[13 15 17 19 21]的方法，及命令语句。
第一种：
>> A=[13 15 17 19 21]
A =
13 15 17 19 21
>> A=[13,15,17,19,21]
A =
13 15 17 19 21 第三种：
>> A=linspace(13,21,5)
A =
13 15 17 19 21 第四种：
>> A=13:2:21
A =
13 15 17 19 21
matlab 作业
实验一 1、熟悉MATLAB的窗口结构（命令窗口、历史命令窗口、工作区窗口、当前目录 窗口） 2、掌握命令窗口中基本命令的使用 3、在命令窗口中，给定圆的半径r，求得圆的周长c和面积s，并查看工作区窗 口的变化 4、将r,c,s变量保存到磁盘文件abc.mat中，并删除内存变量r,c,s，查看工作 区窗口的变化 5、将abc.mat文件中变量装入内存，查看工作区窗口的变化 6、将历史命令窗口中的命令再装入命令窗口中使用 7、改变当前目录，查看当前目录窗口的变化 8、掌握命令窗口中 cd,quit,help,date,dir,ls,what,who,clocl,fix(clock),format,save,loa d,clc,clear等命令的使用 9、注意各种MATLAB版本的差别 实验二 1、在命令窗口中，输入长方形的长和宽，求长方形的周长和面积 2、输入三角形的三条边（要满足构成三角形的条件），求三角形的周长和面积 3、掌握MATLAB中各标准函数的使用（sin,cos,sind,fix,mod,…） 4、用fprintf输出各种类型的数据（如fprintf('a=%d\n',123) a=123 >> fprintf('b=%f\n',123.456) b=123.456000 >> fprintf('c=%c\n','A') c=A……) 实验三 1、在编辑窗口中：输入学生成绩，输出该成绩的等级。等级规定如下：[90， 100]为A等，[80，90)为B等，[70，80)为C等，[60，70)为D等，[0，60)为E等。要求用if和 switch两种方法实现。 2、商场购物，100件以下，不优惠，100~199件95折，200~399件90折，400~799 件85折，800~1499件80折，1500件以上，75折。输入所购货物的单价、件数，求实际付款数目。要求用if和 switch两种方法实现（在编辑窗口中实现）。 实验四 1、求两个正整数的最大公约数和最小公倍数（在编辑窗口中实现，命令窗口中 调用）。 2、求100~300内所有素数（在编辑窗口中实现，命令窗口中调用）。
matlab期末大作业
电气学科大类 Modern Control Systems Analysis and Design Using Matlab and Simulink Title: Automobile Velocity Control Name: 巫宇智 Student ID: U200811997 Class：电气0811
电气0811 巫宇智 Catalogue Preface (3) The Design Introduction (4) Relative Knowledge (5) Design and Analyze (6) Compare and Conclusion (19) After design (20) Appendix (22) Reference (22)
Automobile Velocity Control 1.Preface: With the high pace of human civilization development, the car has been a common tools for people. However, some problems also arise in such tendency. Among many problems, the velocity control seems to a significant challenge. In a automated highway system, using the velocity control system to maintain the speed of the car can effectively reduce the potential danger of driving a car and also will bring much convenience to drivers. This article aims at the discussion about velocity control system and the compensator to ameliorate the preference of the plant, thus meets the complicated demands from people. The discussion is based on the simulation of MATLAB. Key word: PI controller, root locus
北邮模式识别课堂作业答案(参考)
第一次课堂作业 1.人在识别事物时是否可以避免错识 2.如果错识不可避免，那么你是否怀疑你所看到的、听到的、嗅 到的到底是真是的，还是虚假的 3.如果不是，那么你依靠的是什么呢用学术语言该如何表示。 4.我们是以统计学为基础分析模式识别问题，采用的是错误概率 评价分类器性能。如果不采用统计学，你是否能想到还有什么合理地分类 器性能评价指标来替代错误率 1.知觉的特性为选择性、整体性、理解性、恒常性。错觉是错误的知觉，是在特定条件下产生的对客观事物歪曲的知觉。认知是一个过程，需要大脑的参与.人的认知并不神秘，也符合一定的规律，也会产生错误 2.不是 3.辨别事物的最基本方法是计算 . 从不同事物所具有的不同属性为出发点认识事物. 一种是对事物的属性进行度量，属于定量的表示方法(向量表示法 )。另一种则是对事务所包含的成分进行分析，称为定性的描述(结构性描述方法)。 4.风险 第二次课堂作业 作为学生，你需要判断今天的课是否点名。结合该问题(或者其它你熟悉的识别问题，如”天气预报”)，说明: 先验概率、后验概率和类条件概率 按照最小错误率如何决策 按照最小风险如何决策 ωi为老师点名的事件,x为判断老师点名的概率 1.先验概率: 指根据以往经验和分析得到的该老师点名的概率,即为先验概率 P(ωi ) 后验概率: 在收到某个消息之后，接收端所了解到的该消息发送的概率称为后验概率。 在上过课之后,了解到的老师点名的概率为后验概率P(ωi|x) 类条件概率:在老师点名这个事件发生的条件下,学生判断老师点名的概率p(x| ωi ) 2. 如果P(ω1|X)>P(ω2|X)，则X归为ω1类别 如果P(ω1|X)≤P(ω2|X)，则X归为ω2类别 3.1)计算出后验概率 已知P(ωi)和P(X|ωi)，i=1,…，c，获得观测到的特征向量X 根据贝叶斯公式计算 j=1,…，x
matlab上机作业
第四次 上机作业 1、 从键盘输入一个4位整数，按照如下规则加密后输出。加密规则：每位数字 都加上7，然后用和除以10的余数取代该数字；再把第一位与第三位交换，第二位与第四位交换。 Clear X=ones(1,4); X （1）=input(’输入第一位：‘)； X （2）=input(’输入第二位：‘)； X （3）=input(’输入第三位：‘)； X （4）=input(’输入第四位：‘)； X=rem(7+x,10); Y=1000.*x(3)+100.*x(4)+10.*x(1)+x(2) 2、 分别用if 和switch 语句实现以下计算，其中a 、b 、c 的值从键盘输入。 ??? ? ??? <≤+<≤+<≤++=5 .55.3, ln 5.35.1, sin 5.15.0,2x x c b x x b a x c bx ax y c a=input(‘请输入a ：’)； b=input(‘请输入b ：’)； c=input(‘请输入c ：’)； If(x>=0.5&&x<=1.5) y=a.*x^2+b.*x+c Elseif(x>=1.5&&x<=3.5) y=a.*(sin(b))^c+x
Elseif(x>=3.5&&x<=5.5) y=log(abs(b+c./x)) end a=input(‘请输入a：’)； b=input(‘请输入b：’)； c=input(‘请输入c：’)； Switch x case(x>=0.5&&x<=1.5) y=a.*x^2+b.*x+c case(x>=1.5&&x<=3.5) y=a.*(sin(b))^c+x case(x>=3.5&&x<=5.5) y=log(abs(b+c./x)) end 3、产生20个两位随机整数，输出其中小于平均值的偶数。Clear al ;close all ;clc; X=fix(rand(1,20)*89)+10; Disp([‘20个随机数是：’,num2str(x)]); X1=mean(x); Disp([‘平均值为：’,num2str(x1)]); N=find(rem(x,2)==0&x 系统仿真导结课作业 一、概述 建模与仿真技术已成功地应用于航空航天、生产制造、交通运输、信息、生物、医学、材料、能源、教育、军事、社会、经济等众多领域;并成功地应用于产品研制的全生命周期,包括需求分析、方案论证、概念设计、初步设计、详细设计、生产制造、试验试飞、运行、维护、训练等各个阶段。仿真科学与技术正是从其广泛的应用中获得了日益强大的生命力，而仿真技术的发展反过来使得其得到愈来愈广泛的应用。 广义而言,仿真是采用建模的方法和物理的方法对真实环境客观事物进行 抽象、映射、描述和复现。基于系统原理、理论、定律、系统数据等应用计算机技术、软件技术和信息技术建立仿真环境(虚拟环境) ,在仿真环境中对客观事物进行研究。客观事物包括真实环境中的实体/系统、自然环境(地形、大气、海洋、空间)、和人的行为(操作、决策、推理)。仿真环境包括模型、数据、软件、物理效应设备、计算机等。 计算机仿真的三要素是系统、模型、计算机，三个基本活动是模型设计、模型执行、模型分析。计算机仿真的三要素和三个基本活动的关系关系如图1 所示。 系统 模型设计模型分析 模型执行 模型计算机 图1 从模型设计到模型分析经历的过程，即对实物进行仿真可概括为以下几个方面： 1) 问题的描述; 2) 建立概念模型; 3) 建立仿真模型; 4) 收集数据; 5) 编写程序; 6) 在计算机进行模型试验; 7) 模型和数据的验证; 8) 仿真结果显示; 9) 仿真结果分析和评估。 仿真是建立模型在计算机上运行,但这属于数学仿真,随着技术的发展,许多应用 领域建立仿真系统时除了模型外还要求将实物和人员包含在仿真回路中。 由此可见，基于仿真设计与传统设计的方法和流程两者有很大区别(图2) ,基于仿真设计可以在计算机上建立虚拟样机,对产品的外形、结构、强度、动力 一、 单项选择题(本大题共5小题，共25.0分) 1. 传播信息的电缆、电磁波等称为： (5.0分 ) A ：媒介 选择： B B ：媒材 C ：媒体 D ：以上三个均不对 2. 在演示文稿制作时，如果在每一页幻灯片上添加版权信息，则可以在下列哪一种视图模式中进行操作设置： (5.0分 ) A ：讲义母版 选择： B B ：幻灯片母版 C ：备注母版 D ： 备注页视图 3. 在演示文稿中，如果仅需要查看幻灯片内的文字信息，而不需要观看其缩略图内容，则可以选择下列哪一个选项卡： (5.0分 ) A ：大纲选项卡 选择： A B ：幻灯片选项卡 C ：属性选项卡 D ：以上三个均不对 4. 在计算机中用构成一幅图的位置、维数、大小和形状的指令来进行描述的是： (5.0分 ) A ：图形 选择： B B ：图像 C ：图形与图像 D ：以上三个均不对 5. 适用于表现比较细节、层次以及丰富色彩的是 ： (5.0分 ) A ：图形 选择： C B ：图像 C ：图形与图像 D ：以上三个均不对 二、 多项选择题(本大题共2小题，共10.0分) 1. 多媒体技术所涉及到以下哪些形式: (5.0分 ) 选择：A,B,C A ：媒介 B ：媒材 C ：媒体 D ：以上三个均不对 2. 评价多媒体课件质量的指标有 : (5.0分 ) A ：科学性 选择： A,B,C,D B ：技术性 C ：艺术性 D ：教育性 三、 判断题(本大题共3小题，共15.0分) 1. 多媒体课件在开发时其艺术性必须服务于教育性。 (5.0分) √ × 选择：√ 2. 多媒体课件在开发时要求必须使用多种媒体元素。 (5.0分) √ × 选择：× 3. 多媒体课件脚本要严格按照模板制作，不能有个性化调整。 (5.0分) √ × 选择：× 四、 开题报告 论文题目：基于matlab的时域有限差分法的电磁仿真研究（10分） 学院：电气工程及其自动化学院学号：1103000105姓名：__杨志刚___ 一、论文选题的目的和意义（300字以内；15分） 时域有限差分法，因具有多种优点被运用到电磁场理论研究的各个方面，而且其使用成效和应用领域还在迅速扩大和提高，在现代电磁场理论研究中具有很大的重要性和很强的可操作性。但是同时这种方法也存在一定的缺陷，主要表现在对无边界问题需要吸收边界条件处理，有色散误差，消耗内存大等方面。本课题在利用时域有限差分法对一些实际的算例进行实验仿真和验证，同时对这种方法在解决实际问题的缺陷进行一定程度的研究和分析。 Matlab作为一种工程仿真工具得到了广泛应用。用于时域有限差分法，可以简化编程，使研究者的研究重心放在FDTD法本身上，而不必在编程上花费过多的时间。 二、国内外关于该论题的研究现状和发展趋势（500字以内；15分） 时域有限差分方法作为一种典型的全波时域分析方法，因其原理直观、编程简便、实用性强在目前的计算电磁学领域内被人们广泛深入地研究，并取得巨大应用成功的方法。时域数值技术的一个突出优点是可以给出关于问题空间的丰富的时域信息，而且经过简单的时频变换，即可得到宽带范围的频域信息，相对频域方法显著地节约了计算量。最近几十年，是电磁场数值计算时域技术蓬勃发展的时期，各具优势和特色的新颖时域算法层出不穷。 但是到目前为止国内关于时域有限差分法中的PML 算法文献较少，其中绝大多数文献集中在综述和应用方面。而在国际的学报和杂志上对于这方面的文献非常多。时域有限差分法经过了三十年多年的高速发展之后，仍然还是计算电磁学制高点的研究热潮，而且其应用的范围和成效还在迅速的扩大和提高。本课题正是利用时域有限差分法的基础理论，利用matlab对一些实际的电磁场问题进行仿真研究。 三、论文的主攻方向、主要内容、研究方法及技术路线（1000字左右；40分） 通过对时域有限差分法理解基础之上，利用matlab仿真软件按照这种方法编程，实现对三种情况下的电磁场情况的仿真研究。 MATLAB仿真 期末大作业 姓名：班级：学号：指导教师： 2012春期末大作业 题目：设单位负反馈控制系统前向通道传递函数由)()(21s G s G 和串联，其中： ) 1(1)()(21++==s A s G s K s G A 表示自己学号最后一位数（可以是零），K 为开环增益。要求： （1）设K=1时，建立控制系统模型，并绘制阶跃响应曲线（用红色虚线，并标注坐标和标题）；求取时域性能指标，包括上升时间、超调量、调节时间、峰值时间； （2）在第（1）问中，如果是在命令窗口绘制阶跃响应曲线，用in1或者from workspace 模块将命令窗口的阶跃响应数据导入Simulink 模型窗口，用示波器显示阶跃响应曲线；如果是在Simulink 模型窗口绘制阶跃响应曲线，用out1或者to workspace 模块将Simulink 模型窗口的阶跃响应数据导入命令窗口并绘制阶跃响应曲线。 （3）用编程法或者rltool 法设计串联超前校正网络，要求系统在单位斜坡输入信号作用时，速度误差系数小于等于0.1rad ，开环系统截止频率s rad c /4.4''≥ω，相角裕度大于等于45度，幅值裕度大于等于10dB 。 仿真结果及分析： （1）、（2）、将Simulink模型窗口的阶跃响应数据导入命令窗口并绘制阶跃响应曲线 通过在Matlab中输入命令： >> plot(tout,yout,'r*-') >> title('阶跃响应曲线') 即可得出系统阶跃响应曲线，如下： 求取该控制系统的常用性能指标：超调量、上升时间、调节时间、峰值时间的程序如下： G=zpk([],[0,-1],5)。 S=feedback(G,1)。 数理经济学第一次作业题 董志勇 1. You can write your answer in either English or in Chinese or in both. 2. You have to do this paper independently; giving identical answers will be assumed cheating, which will lead to severe punishment. 3. Answer to the point; be precise and avoid unnecessary length. 1、单个垄断企业分别在两个市场上实现价格歧视，求证其曲线从下面与各曲线的水平之和相交。 2、设，为两种投入要素，共同生产一种产品，该产品的价格是，生产函数是，的价格是，的价格是，，是使得利润最大化的要素投入量，求证：和不会同时小于0。 3、生产函数是，求出要素需求曲线，对根据该生产函数进行生产而且追求利润最大化的企业进行比较静态分析，证明对该企业而言，交叉效应项的符号为负。 4、一个垄断者在两个市场上销售同一种产品，市场1的收益为，市场2的收益为，两个市场的总成本为，假设在市场1中，每销售一个单位的产品要征税t。 (a)证明t增加将会使得市场1销售得产品量降低。 (b)有关t增加对市场2中的产品的影响，最大化假设意味着什么？ (c)即使假设二阶条件总是成立，市场1中的税收的增加也可能导致总产出，请说明在什么情况下会是这样的。 (d)假设根据政府法律，市场2的产出固定在前面所述的利润最大化水平上，与不加管制的情况相比，在受到政府管制的情况下，市场1中的产出对税收增加得反应在绝对量上较小，请加以证明，并对此进行直观的经济解释。 5、垄断者在两个市场中销售其产品，收入函数分别为和。总产出的函数。在两个市场中，每销售一个单位的产品，征收同样的税收t。 (a)求解，和。如果有解的话，利润最大化隐含了哪些偏导的符号？ (b)假设产出固定不变，求解，能够确定该式的符号？ 第一次上机实习 1．编写一个程序ex1_1.m，对输入的方阵计算其伴随矩阵。 a)要求第一行帮助行注明这个函数的用途，能在matlab命令窗口可以用help或lookfor 查到该帮助； b)对输入的变量进行检测，如果输入的不是方阵或者矩阵元素不是数，则报错,退出程 序。 function A=ex1_1(X) % 建立函数A=ex1_1(X)，求矩阵X的伴随矩阵 [a,b]=size(X); if a~=b||~isnumeric(X) || ~ismatrix(X) error('输入的不是方阵或者矩阵元素不是数!'); end for i=1:a for j=1:b A(i,j)=(-1)^(i+j)*X(i,j)*det(X([1:i-1,i+1:end],[1:j-1,j+1:end])); end end A=A'; 2．已知表达式y=6*x^5+4*x^3+2*x^2-7x+10，x的范围是[0，100]，使用三阶拟合和五阶的方法得出多项式的表达式，并编程在图中绘制出原曲线、三阶拟合和五段拟合的曲线。 clc clear all x=linspace(0,100,2000); p1=[6 0 4 2 -7 0]; y=polyval(p1,x); p3=polyfit(x,y,3); p5=polyfit(x,y,5); y1=polyval(p3,x); y2=polyval(p5,x); plot(x,y,'r-o') hold on grid on plot(x,y1,'b-.',x,y2,'y--*'); grid on 三条曲线 放大图（原曲线与五阶拟合完全重合）： 3．读如图像5.tif（450*450），做fourier变换， a)将变换之后所得到的矩阵只保留前10行10列，后10行10列的数据，其他位置都 设置为0，然后再作反fourier变换，观察结果； 原图： 选 题 说 明 本人选做第2、4、5、9、11、12、13、14、16、19、24 题。 作业内容题目2：问题描述：在[0 ， 2π]范围内绘制二维曲线图y=cos(5x)*sin(x) （1）问题分析 这是一个二维绘图问题，先写出x的取值范围，再用plot函数画出y的图像。 （2）软件说明及源代码 >> x = 0:pi/100:2.*pi; y=cos(5*x).*sin(2*x); >> plot(x,y) （3）实验结果 题目4：问题描述：创建符号函数并求解，要求写出步骤和运行结果 (1)创建符号函数f=ax2+bx+c (2)求f=0的解 （1）问题分析 这是符号计算问题，首先要确定符号变量，然后创建符号函数，最后利用subs函数求解特值。 （2）软件说明及源代码 >> syms a b c x f; f=a*x^2+b*x+c; subs(f,0) （3）实验结果 ans = c 题目5：问题描述：求积分 （1）问题分析 这是符号计算的积分求解问题，首先需要确定符号变量，然后利用int函数计算积分。 （2）软件说明及源代码 >> syms x y; y=sqrt(1-2*sin(2*x)); >> int(y,x,0,pi/2) （3）实验结果 ans = ellipticE(-pi/4, 4)*1i - ellipticE(pi/4, 4)*1i - ellipticE(-pi/6, 4)*2i + ellipticE(pi/6, 4)*2i 题目9：问题描述：按水平和竖直方向分别合并下述两个矩阵： （1）问题分析 这是考查矩阵的基本操作，首先定义矩阵，然后合并矩阵。 （2）软件说明及源代码 >> A=[1,0,0;1,1,0;0,0,1]; B=[2,3,4;5,6,7;8,9,10]; >> a=[A,B],b=[A;B] （3）实验结果 a = 1 0 0 2 3 4 1 1 0 5 6 7 0 0 1 8 9 10 b = 1 0 0 1 1 0 0 0 1 2 3 4 5 6 7 8 9 10 题目11：问题描述：计算z＝yx2+3y2x+2y3的和： （1）问题分析 这是符号计算问题，首先确定符号变量，然后构造函数，最后利用diff函数进行求导。 （2）软件说明及源代码 >> syms x y z; >> z=y*x^2+3*y^2*x+2*y^3; >> diff(z,y,1),diff(diff(z,y,1),x,1) （3）实验结果 ans = x^2 + 6*x*y + 6*y^2 1 一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全部水仙花数。 for m=100:999 m1=fix(m/100); %求m的百位数字 m2=rem(fix(m/10),10); %求m的十位数字 m3=rem(m,10); %求m的个位数字 if m==m1*m1*m1+m2*m2*m2+m3*m3*m3 disp(m) end end 2.从键盘输入若干个数，当输入0时结束输入，求这些数的平均值和它们之和。 sum=0; n=0; val=input('Enter a number (end in 0):'); while (val~=0) sum=sum+val; n=n+1; val=input('Enter a number (end in 0):'); end if (n > 0) sum mean=sum/n end 3. 若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6是完数。求[1,500]之间的全部完数。 for m=1:500 s=0; for k=1:m/2 if rem(m,k)==0 s=s+k; end end if m==s disp(m); end end 4. 从键盘上输入数字星期，在屏幕上显示对应英文星期的单词。 function week n=input('input the number:'); if isempty(n) errror('please input !!') end if n>7|n<1 error('n between 1 and 7') end switch n case 1 disp('Monday') case 2 disp('Tuesday') case 3 disp('Wednesday') case 4 disp('Thursday') case 5 disp('Friday') case 6 disp('Saturday') case 7 disp('Sunday') end 5. 某公司销售电脑打印机的价格方案如下： （）如果顾客只买一台打印机，则一台的基本价格为$150。 （）如果顾客购买两台以上打印机，则第二台价格为$120。 （）第三台以后，每台$110。 写一段程序分别计算出购买1--10台打印机所需的钱数。打印机台数可以在程序开始处指定，或通过input命令读入。运行程序，计算出购买10台打印机的总价格。 写出程序，生成分别购买1--10台打印机所需价格的图表（使用fprintf命令输出图表，不允许手算）。 x=input('请输入购买的打印机台数：'); for m=1:x if m<=1 y(m)=150*m; elseif m<=2 y(m)=150+120*(m-1); else y(m)=150+120+110*(m-2); y(1,m)=y(m); end end y(x) plot(1:m,y,'r*--') 语文园地二 1、精湛财富属于朴实器重喉结例外赞赏闷闷不乐狼吞虎咽面呈菜色 2、维持生计得到宽容回忆往事动作熟练神情严肃生活拮据 3、(1)行果 (2)言而有信 (3)勿施于人学而不厌，诲人不倦 4、(1)深信(2)自信(3)守信(4)诚信 (1)寡言少语(2)豪言壮语(3)花言巧语(4)胡言乱语 5、(1)严肃道德难题，把车留下来还是还给库伯 （2）小明这次语文考得很差，试卷发下来后，他战战兢兢地递给了妈妈。妈妈的脸顿时由晴转阴，大声嚷道：“你上课都听了些什么呀？这么简单的题目也会错？” 6、(1)不小心，文中指小男孩不小心掉进了大西洋。肯定，一定文中指船长肯定孩子掉进了大海。(2)肯定/确定用力/尽力/奋力的确/确实 (3)他实在游不动了，准备放弃，加上没人发现他掉进大海里，他觉得自己活不下去了认为船长一定会返航去救他的小男孩相信船长一定会去救他，说明船长被孩子所信任 (4)一个人能被他人相信也是一种幸福。他人在绝望时想起你，相信你会给予拯救更是一种无可替代的幸福，我们要真诚地去对待别人。 九、自然之道 1、沙滩侦察蠢事棒球帽愚不可及无遮无拦气喘吁吁 2、啄逐逐啄彻切彻切 4、(1)这个词突出了向导非常不情愿，生怕惹出麻烦，对同伴们善良的呼喊非常理解，但也感到非常无奈 (2)这个词突出了向导悲伤地叹息，对自 己之前的所作所为感到非常羞愧，表达了“自己”自责、悔恨的心理。5、(1)看着说飞走答道抱起 (2)R 艹 5 ① T 页 7 ③ (3)眼前的一幕是指，当幼龟踌躇不前时，一只嘲鸫突然飞来，它用尖嘴啄幼龟的头，企图把它拉到沙滩上去。 (4) ③ 十一、蝙蝠和雷达 1、清朗蝙蝠飞蛾避开蚊子揭开障碍荧光屏 2、bǔ líng qī zhuàng cāng ruì 3、启发灵敏(灵活) 仿照清晰 4、横七竖八地拉了许多绳子系着许多铃铛 5.(1)√×× (2)①波浪②天上的云/春天的风小羊/妈妈的手飘来飘去/吹在脸上真舒服。 （3）蝙蝠用嘴发射超声波，超声波向四周传播，遇到障碍物就反射回来，反射回来的那部分被蝙蝠用耳朵接收，经过大脑处理，就知道哪里有障碍物了。 6、嘴巴超声波耳朵青蛙－电子蛙眼蜻蜓－飞机蝴蝶－卫星螳螂－镰刀顺风耳－电话 2015.1.9 上机作业题报告 JONMMX 2000 1.Chapter 1 1.1题目 设S N =∑1j 2?1 N j=2 ，其精确值为 )1 1 123(21+--N N 。 （1）编制按从大到小的顺序1 1 131121222-+ ??+-+-=N S N ，计算S N 的通用程序。 （2）编制按从小到大的顺序1 21 1)1(111222-+ ??+--+-= N N S N ，计算S N 的通用程序。 （3）按两种顺序分别计算64210,10,10S S S ,并指出有效位数。（编制程序时用单精度） （4）通过本次上机题，你明白了什么？ 1.2程序 1.3运行结果 1.4结果分析 按从大到小的顺序，有效位数分别为：6，4，3。 按从小到大的顺序，有效位数分别为：5，6，6。 可以看出，不同的算法造成的误差限是不同的，好的算法可以让结果更加精确。当采用从大到小的顺序累加的算法时，误差限随着N 的增大而增大，可见在累加的过程中，误差在放大，造成结果的误差较大。因此，采取从小到大的顺序累加得到的结果更加精确。 2.Chapter 2 2.1题目 （1）给定初值0x 及容许误差ε，编制牛顿法解方程f(x)=0的通用程序。 （2）给定方程03 )(3 =-=x x x f ,易知其有三个根3,0,3321= *=*-=*x x x ○1由牛顿方法的局部收敛性可知存在,0>δ当),(0δδ+-∈x 时，Newton 迭代序列收敛于根x2*。试确定尽可能大的δ。 ○2试取若干初始值，观察当),1(),1,(),,(),,1(),1,(0+∞+-----∞∈δδδδx 时Newton 序列的收敛性以及收敛于哪一个根。 （3）通过本上机题，你明白了什么？ 2.2程序 利用Matlab 实现Romberg 数值积分算法 一、内容摘要 针对于某些多项式积分，利用Newton —Leibniz 积分公式求解时有困难，可以采用数值积分的方法，求解指定精度的近似解，本文利用Matlab 中的.m 文件编写了复化梯形公式与Romberg 的数值积分算法的程序，求解多项式的数值积分，比较两者的收敛速度。 二、数值积分公式 1.复化梯形公式求解数值积分的基础是将区间一等分时的Newton —Cotes 求积公式: I =(x)[f(a)f(b)]2 b a b a f dx -≈ +? 其几何意义是，利用区间端点的函数值、与端点构成的梯形面积来近似(x)f 在区间[a,b]上的积分值，截断误差为： 3" (b a)()12 f η-- (a,b)η∈ 具有一次的代数精度，很明显，这样的近似求解精度很难满足计算的要求，因而，可以采用将积分区间不停地对分，当区间足够小的时候，利用梯形公式求解每一个小区间的积分近似值，然后将所有的区间加起来，作为被求函数的积分，可以根据计算精度的要求，划分对分的区间个数，得到复化梯形公式： I =1 1 (b a)(b a) (x)dx [f(a)f(b)2(a )]2n b a k k f f n n -=--≈+++∑? 其截断误差为： 2" (b a)h ()12 R f η--= (a,b)η∈ 2.Romberg 数值积分算法 使用复化的梯形公式计算的数值积分，其收敛速度比减慢，为此，采用Romberg 数值积分。其思想主要是，根据I 的近似值2n T 加上I 与2n T 的近似误差，作为新的I 的近视，反复迭代，求出满足计算精度的近似解。 用2n T 近似I 所产生的误差可用下式进行估算： 12221 ()3 n n n I T T T -?=-=- 新的I 的近似值： 122 n n j T T -=?+ j =（0 1 2 ….） Romberg 数值积分算法计算顺序 i=0 (1) 002T i=1 (2) 102T (3) 012T i=2 (4) 202T (5) 112T (6) 022T i=3 (7) 302T (8) 212T (9) 122T (10) 032T i=4 (11) 402T (12) 312T (13) 222T (14) 132T … … … … 其中，第一列是二阶收敛的，第二列是四阶收敛的，第三列是六阶收敛的，第四列是八阶收敛的，即Romberg 序列。 杭州电子科技大学信息工程学院《科学计算与MATLAB》期末大作业 给出程序、图、作业分析，程序需加注释。 1. 试编写名为fun.m 的MATLAB 函数，用以计算下述的值： ?? ? ??-<->=t t n t t t n t f 的)4/sin()(si 对所有)4/sin(其他情况)sin(的)4/sin()(si 对所有)4/sin()(ππππ 绘制t 关于函数f(t)的图形，其中t 的取值范围为ππ66≤≤-t ，间距为10/π。 function y=fun()%定义函数 % t=-6*pi:pi/10:6*pi; %定义变量范围 y = (sin(pi/4)).*(sin(t)>sin(pi/4))+(sin(-pi/4)).*(sin(t) 2.解以下线性方程组 ??? ??=+=++=--3 530 42231 321321x x x x x x x x A=[2 -1 -1;1 1 4;3 0 5];%输入矩阵 B=[2;0;3]; %输入矩阵 X = A\B %计算结果 3.已知矩阵? ? ??? ???? ???=44434241 3433323124232221 14131211A 求： (1)A(2:3,2:3) (2)A(:,1:2) (3)A(2:3,[1,3]) (4)[A,[ones(2,2)；eye(2)]] A=[11 12 13 14;21 22 23 24;31 32 33 34;41 42 43 44];%输入矩阵A(2:3,2:3) %输出矩阵 A(:,1:2) %输出矩阵 A(2:3,[1,3]) %输出矩阵 [A,[ones(2,2);eye(2)]] %输出矩阵 人教版小学四年级下册语文课堂作业本答案 一、古诗词三首 1、敬亭山洞庭湖青螺 2.×√××√ 3.(1)完了，没有了悠闲孤独众孤因怀才不遇而感到的孤独、寂寞的心情 (2)月光下的洞庭湖洞庭湖的君山翠白青 (3)胜过仿佛，像蓝草，一种植物，叶蓝绿色，可提取青蓝色染料5、江南春唐杜牧千里莺啼绿映红，水村山郭酒旗风。南朝四百八十寺，多少楼台烟雨中。 二、桂林山水 1、无瑕攀登屏障骆驼危峰兀立波澜壮阔奇峰罗列 2、sì juàn dào m?ng 3、波澜壮阔水平如镜峰峦雄伟红叶似火衣服（图画）山峰（山峦） 4、(1)秋天，落叶仿佛像一个会跳舞的精灵，从树上翩翩飞舞而下。 (2)老房子的围墙真危险呐，手一碰，墙壁上的石灰好像要掉下来似的。 5、(1)M 辶 6 ④ D 斤 7 ② (2)静清绿奇秀险 (3) 绿树红花竹筏小舟云雾山水 6、田野上的野花真多啊，多得像天上的星星数也数不清；田野上的野花真美啊，美得仿佛那是一条五彩的锦缎；田野上的野花真香啊，香得引来群群蜜蜂争相采蜜！（我游览过水平如镜的西湖，玩赏过飞泻而下的瀑布，却从没看到过秦淮河的水。秦淮河的水真凉啊，凉得让你想马上和河水来 个亲密的接触；秦淮河的水真清啊，清得就像是一块彩色的大玻璃；秦淮河的水真亮啊，当河边的万家灯火映在水中，像是水底铺满了熠熠生辉的大钻石，使河水变得无比明亮、美丽；秦淮河的水真野啊，刚才还是块大玻璃，现在就波涛起伏，使得河上的画舫也一个劲儿地摇摆，再加上朱雀桥洞内随风而舞的狂草，更体现了秦淮河的粗犷之美；秦淮河的水真神啊，要不，它怎么会吸引来那么多的小雨点儿？想当年大自然在创造世界时，一定挥舞着最大的画笔，而放弃了那些精细的刻刀，毫无雕琢地在南京勾勒出了这样一条粗犷、神奇、美丽的秦淮河吧？） 三、记金华的双龙洞 1、浙江臀部稍微额角擦伤蜿蜒依据石钟乳 2、一簇臀部稍微蜿蜒民族宫殿树梢大碗 3、(1)②① (2)①② 4、（1）小妹妹总是时而笑，时而哭，时而忧，时而喜，真是捉摸不透。（六 月的天气，时而风，时而雨，时而凉，时而热，真是变化多端。）（2）即使这件事情再困难，我们也要尽自己最大的力量去做好它。 5、(1)根据名称欣赏 (2)，，、、， (3)在洞顶的那两条龙，一条黄龙，一条青龙。 (4)蜿蜒在洞顶的双龙变化多端，颜色各异的石钟乳和石笋。 6、怎样美的笔盒呢？让人一看就爱不释手，是这样美的笔盒。（怎样简陋的小屋呢？裂缝的石灰墙，风化的水泥地，一扇破旧的窗，一张油漆脱落的书桌，根本找不到一样值钱的东西，就是这样一间简陋的小屋。） 语文园地一 1、洞庭玩赏扩散泰山拥挤孔隙仰卧油桐水平如镜拔地而起突兀 第一次 上机作业 1．写出完成下列操作的命令： 1) 将矩阵A 的第2～5行中的第1，3，5列元素赋给矩阵B B=A(:,1:3:5) 2) 删除矩阵A 的第7号元素 A(7)=[] 3) 将矩阵A 的每个元素值加30 A=A+30 4) 求矩阵A 的大小和维数 Whos size(A) 5) 将含有12个元素的向量x 转换成3×4矩阵 Y=reshape(x,3,4) 2．下列命令执行后，L1、L2、L3、L4的值分别是多少？ A=1:9; B=10－A; L1=A= =B; L2=A< =5; L3=A>3&A<7; L4=find(A>3&A<7); L1= 0 0 0 0 1 0 0 L2=1 1 1 1 1 0 0 L3=0 0 0 1 1 1 0 L4=4 5 6 3．已知： A=????? ???????---14.35454.9632053256545410778.01023 完成下列操作： (1)取出A的前3行构成矩阵B，前两列构成矩阵C，右下角3×2子矩阵构成矩 阵D，B与C的乘积构成矩阵E B=A(1:3,:) C=A(:,1:2) D=A(2:4,3:4) E=B*C (2)分别求E 第二次 上机作业 1、 求下列矩阵的主对角线元素、上三角阵、下三角阵、秩、范数、条件数和迹。 （1）????? ???????--=901511250324153211A （2）??????-=2149.824343.0B 1. A=[1,-1,2,3;5,1,-4,2;3,0,5,2;11,15,0,9] D=diag(A) C=triu(A) B=tril(A) E=rank(A) F=trace(A) a1=norm(A,1) a2=norm(A,inf) a3=norm(A,inf) c1=cond(A) c1=cond(A,1) c2=cond(A,2) c3=cond(A,inf) 2. B=[0.43,43,2;-8.9,4,21] D=diag(B) C=triu(B) B=tril(B) E=rank(B) F=trace(B) a1=norm(B,1) a2=norm(B,inf) a3=norm(B,inf) c1=cond(B) c1=cond(B,1) c2=cond(B,2) c3=cond(B,inf) 2、 求矩阵A 的特征值和相应的特征向量。 ???? ??????=225.05.025.0115.011A A=[1,1,0.5;1,1,0.25;0.5,0.25,2] [V ,D]=eig(A) 3、 下面是一个线性方程组： ???? ??????=????????????????????52.067.095.06/15/14/15/14/13/14/13/12/1321x x x （1） 求方程的解。 （2） 将方程右边向量元素3b 改为0.53，再求解，并比较3b 的变化和解的相 对变化。 （3） 计算系数矩阵A 的条件数并分析结论。 A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6] B=[0.95,0.67,0.52] X=inv(A)*b c1=cond(A,1) c2=cond(A,2) c3=cond(A,inf) 4、 利用Matlab 提供的randn 函数生成符合正态分布的10×5随机矩阵A,进行 如下操作： (1)A 各列元素的均值和标准方差 (2)A 的最大元素和最小元素 (3)求A 每行元素的和以及全部元素之和 (4)分别对A 的每列元素按升序、每行按降序排列 X=randn(10,5) M=mean(X) D=std(X) m=max(X) n=min(X) P=sum(X,2) sum(p) 4.10 上机操作步骤 1在MatLab 的命令窗口输入如下命令序列: clf subplot(1,2,1) hold on grid on n=1:1000; m=1./n.*cos(n*pi/2); plot(n,m,'k.') 观察数列的散点图22,当n 趋于无穷大时,数列趋于 0 subplot(1,2,2) hold on grid on n=500:10000; m=1./n.*cos(n*pi/2); plot(n,m,'k.') fplot('0.001',[500,10000]) fplot('-0.001',[500,10000]) axis([500,10000,-0.005,0.005]) 观察图23,当001.0=ε时,可以取N= 1000 ,当n>N 时有επε<< -2 co s n 1n . 图22 图23 2 在MatLab 的命令窗口输入如下命令序列: clf subplot(1,2,1) hold on grid on fplot('x.*x',[1,3]) 观察函数图24, 当2x →时,2x y =的极限是 4 subplot(1,2,2) hold on grid on fplot('x.*x',[1.9,2.1]) fplot('4.001',[ 1.9,2.1]) fplot('3.999',[ 1.9,2.1]) axis([1.9997,2.0005,3.9989,4.0011]) % 调整显示图形的范围是该实验的重点 观察图25,当001.0=ε时, δ取 0.003 δ<-<2 0x 时,001.04<-y ? 图24 图25 3 在MatLab 的命令窗口输入: syms x limit((2.^x-log(2.^x)-1)./(1-cos(x)),x,0) 运行结果为 ans = log(2)^2 理论上用洛必达法则计算该极限: x x x cos 112ln 2lim 0x ---→= 1 4 在MatLab 的命令窗口输入如下命令序列: （1）syms x y=sqrt(x+2)*(3-x)^4/(x+1)^5 diff(y,x) ％求一阶导数 运行结果 ＝y'1/2/(x+2)^(1/2)*(3-x)^4/(x+1)^5-4*(x+2)^(1/2)*(3-x)^3/(x+1)^5-5*(x+2)^(1/2)*(3-x )^4/(x+1)^6 x=1; eval(y) ％求导数在x ＝1处的值 运行结果 1'=x y = 0.8660系统仿真结课作业
西南大学含弘在线课堂PowerPoint多媒体课程第一次作业答案
科技写作结课作业(时域有限差分法的Matlab仿真开题报告)
春MATLAB仿真期末大作业
第一次作业
matlab 上机作业1,中国石油大学
MATLAB大作业
Matlab 上机题及答案
人教版小学四年级下册语文课堂作业本答案
东南大学-数值分析上机题作业-MATLAB版
利用Matlab实现Romberg数值积分算法----系统建模与仿真结课作业
《科学计算与MATLAB》期末大作业
人教版小学四年级下册语文课堂作业本答案
matlab上机作业
matlab上机作业
MATLAB结课作业