matlab第5次实验
实验五(综合性实验)
系统的Simulink仿真分析
一、实验目的
1．熟悉Simulink 的操作方法；
2.掌握绘制系统模型的方法；
3.了解S函数的功能与设计方法。
二、实验设备
计算机、MATLAB软件
三、实验原理
在控制系统中，有初始状态为0的二阶微分方程y''+1.5y'+10y=2u'(t)+10u(t)，其中u(t)为单位阶跃函数，试建立系统模型并仿真。将上述微分方程改写为：y''=2u'(t)+10u(t) -1.5y'-10y。u经微分得到u'，y''经积分得到y'，y'在经积分得到y,而u'、u、y'经代数计算产生y''，根据此原理可建立系统模型并进行仿真，实验过程如下：
四、实验过程
4.1方法1：利用微分/积分器直接构造求解微分方程的模型
1、在MATLAB命令窗口中输入“simulink”进入仿真界面。
2、构建系统：（1）新建Simulink模型窗口（选择“File/New/Model”）,在Simulink Library Browser中将需要的模块拖动到新建的窗口中。
各模块如下：
Math Operations模块库中的Gain模块，它是增益。拖到模型窗口中后，双击
模块，在弹出的对话框中将‘Gain’分别改为‘10’、‘2’、‘1.5’，‘10’表示这4个增益系数。
Continuous模块库中的Integrator模块，它是积分模块；Derivative模块
，它是微分模块。
Math Operations模块库中的Add模块，它是加法模块，默认是两个输入相加，双击该模块，将‘List of Signs’框中的两个加号（++）后输入两个减
号（-），这样就改为了两个个加号，两个减号，用来表示4个信号的叠加。
输入信号在sources模块库，阶跃信号Step。
示波器则Sinks模块库，选择Scope。
3设置系统仿真参数，在模型编辑窗口的Simulation—>configuration parameter —> stop time栏把仿真的停止时间设置为10。
4.仿真操作。双击示波器图标，打开示波器窗口。
4.2方法2利用传递函数模块建模
对方程y''+1.5y'+10y=2u'(t)+10u(t)进行拉普拉斯变换，得到S2Y(S)=2SU(S)+10U(S)-1.5SY(S)-10Y(S)，求解传递函数得到：
在Continuous 模块库中有标准的传递函数（Transfer Fcn ）
模块可供调用。
双击传递函数，在以下界面中Numerator 栏填写[2 10]，Denominator 填写 [1 1.5 10]。得到。
构建如下图的仿真模型。
2设置系统仿真参数，在模型编辑窗口的Simulation —>configuration parameter —> stop time 栏把仿真的停止时间设置为10。
3.仿真操作。双击示波器图标，打开示波器窗口。
10
5.1102)()()(2+++==S S S s U s Y s G
四、实验要求
1.预习实验原理.
2.对实验内容编程,上机运行,记录运行结果(图或曲线或数据).
3.写出实验报告
MATLAB数学实验第二版答案(胡良剑)
数学实验答案 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
MATLAB实验报告(1-4)
信号与系统MATLAB第一次实验报告 一、实验目的 1.熟悉MATLAB软件并会简单的使用运算和简单二维图的绘制。 2.学会运用MATLAB表示常用连续时间信号的方法 3.观察并熟悉一些信号的波形和特性。 4.学会运用MATLAB进行连续信号时移、反折和尺度变换。 5.学会运用MATLAB进行连续时间微分、积分运算。 6.学会运用MATLAB进行连续信号相加、相乘运算。 7.学会运用MATLAB进行连续信号的奇偶分解。 二、实验任务 将实验书中的例题和解析看懂，并在MATLAB软件中练习例题，最终将作业完成。 三、实验内容 1.MATLAB软件基本运算入门。 1). MATLAB软件的数值计算： 算数运算 向量运算：1.向量元素要用”[ ]”括起来，元素之间可用空格、逗号分隔生成行向量，用分号分隔生成列向量。2.x=x0:step:xn.其中x0位初始值，step表示步长或者增量，xn为结束值。 矩阵运算：1.矩阵”[ ]”括起来；矩阵每一行的各个元素必须用”,”或者空格分开； 矩阵的不同行之间必须用分号”;”或者ENTER分开。2.矩阵的加法或者减法运算是将矩阵的对应元素分别进行加法或者减法的运算。3.常用的点运算包括”.*”、”./”、”.\”、”.^”等等。
举例：计算一个函数并绘制出在对应区间上对应的值。 2).MATLAB软件的符号运算：定义符号变量的语句格式为”syms 变量名” 2.MATLAB软件简单二维图形绘制 1).函数y=f(x)关于变量x的曲线绘制用语：>>plot(x,y) 2).输出多个图像表顺序：例如m和n表示在一个窗口中显示m行n列个图像，p 表示第p个区域，表达为subplot(mnp)或者subplot(m,n,p) 3).表示输出表格横轴纵轴表达范围：axis([xmax,xmin,ymax,ymin]) 4).标上横轴纵轴的字母：xlabel(‘x’),ylabel(‘y’) 5).命名图像就在subplot写在同一行或者在下一个subplot前：title(‘……’) 6).输出：grid on 举例1：
实验2-循环结构
实验2 循环结构 程序填空 1. 题目描述：马克思曾经出过这样一道趣味数学题：有30个人在一家小饭馆里用餐，其中有男人、女人和小孩。每个男人花了3先今，每个女人花了2先令，每个小孩花了1先令，一共花去50先令。问男人、女人以及小孩各有几个人。以下是完成此项工作的程序，请将未完成的部分填入，实现其功能，并在计算机上调试程序，以测试填入的部分是否正确。代码： #include
Matlab实验第一次实验答案
实验一Matlab使用方法和程序设计 一、实验目的 1、掌握Matlab软件使用的基本方法； 2、熟悉Matlab的数据表示、基本运算和程序控制语句 3、熟悉Matlab绘图命令及基本绘图控制 4、熟悉Matlab程序设计的基本方法 二、实验内容： 1、帮助命令 使用help命令，查找sqrt（开方）函数的使用方法； 解：sqrt Square root Syntax B = sqrt(X) Description B = sqrt(X) returns the square root of each element of the array X. For the elements of X that are negative or complex, sqrt(X) produces complex results. Remarks See sqrtm for the matrix square root. Examples sqrt((-2:2)') ans = 0 + 1.4142i 0 + 1.0000i
1.0000 1.4142 2、矩阵运算 （1）矩阵的乘法 已知A=[1 2;3 4]; B=[5 5;7 8]; 求A^2*B 解：A=[1 2;3 4 ]; B=[5 5;7 8 ]; A^2*B （2）矩阵除法 已知A=[1 2 3;4 5 6;7 8 9]; B=[1 0 0;0 2 0;0 0 3]; A\B,A/B 解：A=[1 2 3;4 5 6;7 8 9 ]; B=[1 0 0;0 2 0;0 0 3 ]; A\B,A/B （3）矩阵的转置及共轭转置
实验5_循环结构程序设计2题目及解答
4.3 循环结构程序设计2 【注意事项： （1）题中页码均为《C语言程序设计教程》中的页码 （2）题后为原程序文件名 （3）仅提交名为【*.c】的源程序文件，可一起提交，也可逐题提交 调试程序中常见问题及处理方法： （1）出现死循环时：关闭运行窗口即可。 】 1、按照下面的流程图书程序，用for语句实现。 【提示：参见p97—例4.23 ，并按照下面的左流程图书写007张三511.c，右流程图书写007张三512.c 程序运行结果： s=2550 】 2、用for语句编程，求300～500之间的水仙花数（如135：153=13+53+33）。（007张三52.c） 【程序运行结果： 300～500之间的全部水仙花数如下：??370??371??407 】 3、判断输入的数是否为素数。【提示：素数是只能被1和其自身和整除的数，】（007张三53.c） 【提示：参见p114—例4.36 程序运行时：程序运行结果： 请输入一个正整数：25?25不是素数 请输入一个正整数：12?125不是素数 请输入一个正整数：13?13是素数 】 4、判断输入的一个不超过32767的正整数是几位数。（007张三54.c）【提示：按照右侧的N-S图编程 运行时：运行结果： 请输入一个正整数（1～32767）: 135?135是3位数请输入一个正整数（1～32767）: 1356?1356是4位数请输入一个正整数（1～32767）: 28906?28906是5位数】
附加题： 1、求出100～200间的所有素数。（007张三5附加1.c ）。 【 运行结果： 101 103 107 109 113 127 131 137 139 149 151 157 163 167 173 179 181 191 193 197 199 】 2、根据以下公式求π值（直到最后一项小于10 –6 为止）（007张三5附加2.c ） 【 运行结果： pi=3.141590 】 3、百钱百鸡问题中国古代数学家张丘建在他的《算经》中提出了著名的“百钱买百鸡问题”：鸡翁一，值钱五，鸡母一，值钱三，鸡雏三，值钱一，百钱买百鸡，问翁、母、雏各几何？（007张三5附加3.c ） 【提示：在编程之前，需先根据题意列出方程， 如： 设鸡翁x 只，鸡母y 只，鸡雏z 只 则有下列方程： x+y+z=100 和 5x+3y+z/3=100 因为2个方程，3个未知数，故为不定方程，需用穷举法解题。 运行结果： x=0??? y=25??z=75?? x=4??? y=18??z=78?? x=8??? y=11??z=81?? x=12?? y=4???z=84?? 】 参考解答： 1、～2、略 3、 #include
matlab数学实验练习题
Matlab 数学实验 实验一 插值与拟合 实验内容： 预备知识：编制计算拉格朗日插值的M 文件。 1. 选择一些函数，在n 个节点上（n 不要太大，如5 ~ 11）用拉格朗日、分段线性、三次样条三种插值方法，计算m 个插值点的函数值（m 要适中，如50~100）。通过数值和图形输出，将三种插值结果与精确值进行比较。适当增加n ，再做比较，由此作初步分析。下列函数任选一种。 （1）、 ;20,sin π≤≤=x x y （2）、;11,)1(2/12≤≤--=x x y （3）、;22,c o s 10 ≤≤-=x x y （4）、22),exp(2≤≤--=x x y 2.用电压V=10伏的电池给电容器充电，电容器上t 时刻的电压为 ) (0)()(t e V V V t v ---=，其中0V 是电容器的初始电压，τ是充电常数。试由下面 一组t ，V 数据确定0V 和τ。 实验二 常微分方程数值解试验 实验目的： 1. 用MATLAB 软件求解微分方程，掌握Euler 方法和龙格-库塔方法； 2. 掌握用微分方程模型解决简化的实际问题。 实验内容：
实验三地图问题 1.下图是一个国家的地图，为了计算出它的国土面积，首先对地图作如下测量： 以由西向东方向为x轴，由南到北方向为y轴，选择方便的原点，并将从最西边界点到最东边界点在x轴上的区间适当地划分为若干段，在每个分点的y方向测出南边界点和北边界点的y坐标y1和y2，这样就得到了表中的数据（单位mm）。 根据地图的比例我们知道18mm相当于40km，试由测量数据计算该国土 的近似面积，并与它的精确值41288km2比较。
实验二 MATLAB程序设计 含实验报告
实验二 MATLAB 程序设计 一、 实验目的 1.掌握利用if 语句实现选择结构的方法。 2.掌握利用switch 语句实现多分支选择结构的方法。 3.掌握利用for 语句实现循环结构的方法。 4.掌握利用while 语句实现循环结构的方法。 5.掌握MATLAB 函数的编写及调试方法。 二、 实验的设备及条件 计算机一台（带有MATLAB7.0以上的软件环境）。 M 文件的编写： 启动MATLAB 后，点击File|New|M-File ，启动MATLAB 的程序编辑及调试器（Editor/Debugger ），编辑以下程序，点击File|Save 保存程序，注意文件名最好用英文字符。点击Debug|Run 运行程序，在命令窗口查看运行结果，程序如有错误则改正 三、 实验内容 1.编写求解方程02=++c bx ax 的根的函数（这个方程不一定为一元二次方程，因 c b a 、、的不同取值而定） ，这里应根据c b a 、、的不同取值分别处理，有输入参数提示，当0~,0,0===c b a 时应提示“为恒不等式!”。并输入几组典型值加以检验。 （提示：提示输入使用input 函数） 2.输入一个百分制成绩，要求输出成绩等级A+、A 、B 、C 、D 、E 。其中100分为A+，90分～99分为A ，80分～89分为B ，70分～79分为C ，60分～69分为D ，60分以下为E 。 要求：（1）用switch 语句实现。 （2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。 （提示：注意单元矩阵的用法） 3.数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与3相乘再加上1。重复此过程，最终得到的结果为1。如： 2?1 3?10?5?16?8?4?2?1 6?3?10?5?16?8?4?2?1 运行下面的程序，按程序提示输入n=1,2,3,5,7等数来验证这一结论。 请为关键的Matlab 语句填写上相关注释，说明其含义或功能。 4. 的值，调用该函数后，
(完整版)Matlab实验5选择结构程序结构
实验五、选择与循环结构 一、实验目的： 1、 掌握建立和执行M 文件的方法。 2、 掌握利用if 语句实现选择结构的方法。 3、 掌握利用switch 语句实现多分支选择结构的方法。 4、 掌握try 语句的使用。 5、 掌握利用for 语句实现循环结构的方法。 6、 掌握利用while 语句实现循环结构的方法。 7、 熟悉利用向量运算来代替循环的操作方法。 二、实验内容： 1、 列分段函数的值。 ?? ???--≠≠<≤+--≠<-+=其他且且,632,100,6530,6222x x x x x x x x x x x y 要求： (1) 用if 语句实现，分别输出x =-0.5,-3.0,1.0,2.0,2.5,3.0,5.0时的y 值。 提示：x 的值从键盘输入，可以是向量。 %homework_5_1_1.m x=input('请输入x 的值：x='); if (x<0 & x~=-3) y= x.*x + x - 6 elseif (x>=0 & x<10 & x~=2 & x~=3) y=x.*x-5.*x+6 else y=x.*x-x-6 end >> homework_5_1 请输入x 的值：x=[-0.5 -3.0 1.0 2.0 2.5 3.0 5.0] y = -5.2500 6.0000 -6.0000 -4.0000 -2.2500 0 14.0000 (2) 用逻辑表达式实现上述函数。 %homework_5_1_2.m x=input('请输入x 的值：x=') y=(x<0 & x~=-3).*(x.*x+x-6)... +(x>=0 & x<10 &x~=2 &x~=3).*(x.*x-5.*x+6)... +(x>=10 | x==-3 | x==3 | x==2).*(x.*x-x-6) >> homework_5_1_2 请输入x=[-0.5 -3.0 1.0 2.0 2.5 3.0 5.0] x = -0.5000 -3.0000 1.0000 2.0000 2.5000 3.0000 5.0000 y = -6.2500 6.0000 2.0000 -4.0000 -0.2500 0 6.0000
浅析Matlab数学实验报告
数学实验报告 姓名: 班级: 学号: 第一次实验任务 过程: a=1+3i; b=2-i; 结果: a+b =3.0000 + 2.0000i a-b =-1.0000 + 4.0000i a*b = 5.0000 + 5.0000i a/b = -0.2000 + 1.4000i 过程: x=-4.5*pi/180; y=7.6*pi/180; 结果: sin(abs(x)+y)/sqrt(cos(abs(x+y))) =0.2098 心得:对于matlab 中的角度计算应转为弧度。 (1)过程: x=0:0.01:2*pi; y1=sin(x); y2=cos(x); y3=exp(x); y4=log(x); plot(x,y1,x,y2,x,y3,x,y4) plot(x,y1,x,y2,x,y3,x,y4) 结果: (2)过程:>> subplot(2,2,1) >> plot(x,y1) >> subplot(2,2,2) >> plot(x,y2) ./,,,,2,311b a b a b a b a i b i a ?-+-=+=计算、设有两个复数 6,7,5.4)
cos()sin(2=-=++y x y x y x ,其中、计算的图形。 下分别绘制)同一页面四个坐标系)同一坐标系下(、在( x y e y x y x y x ln ,,cos ,sin 213==== >> subplot(2,2,3) >> plot(x,y3) >> subplot(2.2.4) >> subplot(2,2,4) >> plot(x,y4) 结果: 心得:在matlab中,用subplot能够实现在同一页面输出多个坐标系的图像,应注意将它与hold on进行区别,后者为在同一坐标系中划出多条曲线。 5、随机生成一个3x3矩阵A及3x2矩阵B,计算(1)AB,(2)对B中每个元素平方后得到的矩阵C,(3)sinB,(4)A的行列式,(5)判断A是否可逆,若可逆,计算A的逆矩阵,(6)解矩阵方程AX=B,(7)矩阵A中第二行元素加1,其余元素不变,得到矩阵D,计算D。 过程:A=fix(rand(3,3).*10) ; B=fix(rand(3,3).*10);
昆明理工大学MATLAB实验指导书(第二次实验)
************************ MATLAB上机指导书 ************************ 昆明理工大学机电学院 彭用新 2015年3月
实验三符号计算 一、操作部分：在命令窗口执行命令完成以下运算，记录运算结果。 1.findsym：帮助我们获取系统定义的自变量 f= sym('sin(a*x+b*y)'); findsym(f) 2.numden（获取分子分母）, sym2poly,（获取多项式时系数）poly2sym（根据多项式系 数获得符号表达式） [n,d]=numden(sym('x*x+y')+sym('y^2')) p=sym('2*x^3+3*x^2+4'); sym2poly(p) x=[2,3,0,4]; poly2sym(x) 3. collect ：合并同类项；expand：展开多项式；horner: 分解成嵌套形式；factor：因式 分解；simplify: 对表达式化简 syms x y; collect(x^2*y+y*x-x^2-2*x) collect((x+y)*(x^2+y^2+1), y) syms x y; expand((x-2)*(x-4)) syms x;horner(x^3-6*x^2+11*x-6) syms x;factor(x^3-6*x^2+11*x-6) syms x;simplify((x^2+5*x+6)/(x+2)) 4. finverse ：求得符号函数的反函数。 syms x y; finverse(1/tan(x)) f= x^2+y; finverse(f,y) finverse(f) https://www.wendangku.net/doc/689042714.html,pose 求符号函数的复合函数 syms x y; f = 1/(1 + x^2); g = sin(y); compose(f,g) 6. subs :表达式替换。 syms a b;subs(a+b,a,4)
Python实验循环结构
实验五 循环结构（续） 【实验目的】 1．掌握range()函数的使用 2. 掌握For 语句的使用。 3．掌握循环嵌套的使用。 4．掌握一些常用的算法（求素数、穷举法、求和与乘积、递推法）。 【实验内容】 000000000 1、 填空 （1） 在命令行中显示1 到10的数字 range(1,11) 。 （2） 在命令行中显示 1到10的所有偶数 range(2,11,2) 。 （3） 显示所有的大写字符。 for i in range(65,91) : print chr(i) (4) 显示1--100能被3整除的数，并且每5个数打印一行。 i=0 ‘记录能被3整除的个数 for j in range(1,101): if j%3==0: print j, i+=1 if i%5==0: print 2、输出斐波那契级数1、1、2、 3、5、8、13……，此级数项的规律是：前两项的值各为1，从第3项起，每一项是前2项的和。（要求一行输出6项）。 （1） 当某项的数值大于或等于30000时结束。 （2） 求前100项之和。 （1）1 1 2 3 5 8 13 21 34 55 89 144 233 377 610 987 1597 2584 4181 6765 10946 17711 28657 （2 3、编写程序，输入一个数，判断该数是否为完数。例如6是一个完数，因为6=1+2+3。 x=input ("Enter an integer:") i=1 while True : x=x-i i+=1 if x==i: print "A perfect number." break if x
matlab 数学实验 迭代 _ 蛛网图(免积分)
数学实验—实验报告（免积分） 一、实验项目：Matlab实验三—迭代 二、实验目的和要求 a.熟悉MATLAB软件的用户环境，掌握其一般目的命令和MATLAB数组操作与 运算函数； b.掌握MATLAB软件的绘图命令，能够熟练应用循环和选择结构实现各种循环 选择功能； c.借助MATLAB软件的绘图功能，对函数的特性进行探讨，广泛联想，大胆猜 想，发现进而证实其中的规律。 三、实验内容 问题一：将方程53 x x x +-+=改写成各种等价的形式进行迭代 5210 观察迭代是否收敛，并给出解释。 问题二：迭代以下函数，分析其收敛性。 4 f(x)=x-a 使用线性连接图、蛛网图或费根鲍姆图对参数a进行讨论和观察，会得到什么结论？ 问题一： （1）画图 x1=-6:0.01:6; x2=-3:0.01:3; x3=-1:0.01:1; x4=-0.8:0.01:-0.75; y1=x1.^5 +5*x1.^3-2*x1+1; y2=x2.^5 +5*x2.^3-2*x2+1; y3=x3.^5 +5*x3.^3-2*x3+1; y4=x4.^5 +5*x4.^3-2*x4+1; subplot(2,2,1),plot(x1,y1) ,title('图(1)') ,grid on, subplot(2,2,2),plot(x2,y2) ,title('图(2)'),grid on, subplot(2,2,3),plot(x3,y3) ,title('图(3)'),grid on, subplot(2,2,4),plot(x4,y4) ,title('图(4)') ,grid on,
MATLAB第二次上机实验报告
电子科技大学电子工程学院标准实验报告（实验）课程名称MATLAB与数值分析 学生姓名： 学号： 指导教师：
一、实验名称 实验二 线性方程组求解和函数的数值逼近 二、实验目的 通过上机实验，使学生对病态问题、线性方程组求解和函数的数值逼近方法有一个初步的理解。 实验涉及的核心知识点：病态方程求解、矩阵分解和方程组求解、Lagrange 插值。 实验重点与难点：算法设计和MATLAB 编程 三、实验内容 1. 对高阶多项式 ()()() ()()20 1 1220k p x x x x x k ==---=-∏ 编程求下面方程的解 ()190p x x ε+= 并绘图演示方程的解与扰动量ε的关系。 2. 对2 20n =，生成对应的Hilbert 矩阵，计算矩阵的条件数；通过先确定解获得常向量 b 的方法，确定方程组 ()n H x b = 最后，用矩阵分解方法求解方程组，并分析计算结果。 3. 对函数 ()2 1 125f x x = + []1,1x ∈- 的Chebyshev 点 ()()21cos 21k k x n π ?? -= ? ?+? ? ，1,2,,1k n =+ 编程进行Lagrange 插值，并分析插值结果。 四、实验数据及结果分析 1. 对高阶多项式
()()() ()()20 1 1220k p x x x x x k ==---=-∏ 编程求下面方程的解 ()190p x x ε+= 并绘图演示方程的解与扰动量ε的关系。 p=[1,-1]; for i=2:20 n=[1,-i]; p=conv(p,n); % 求多项式乘积 end m=zeros(1,21); % m 的最高次幂为20，有21项 hold on x=1:20; d=[-1,0,0.1,0.5,1]; for i=1:5 delt=d(i); m(2)=delt; y=(roots(p+m))'; % 求多项式的根 plot(x,y,'-o','color',[i/5,i/20,i/10]); end title('方程p(x)=0的解与扰动量delt 的关系') legend('delt=-1','delt=0','delt=0.1','delt=0.5','delt=1') 2468101214161820 010 20 30 40 50 60 方程p(x)=0的解与扰动量delt 的关系 delt=-1delt=0delt=0.1delt=0.5delt=1
实验五__循环结构程序设计
实验五循环结构程序设计 （4学时） 【实验前做的操作】 在本人所用计算机的逻辑盘（如D：盘或E：盘）上建立一个以本人姓名为名称的文件夹，然后在此父文件夹下建立子文件夹，名称为“实验五”，用来零时存放实验五需要上交的有关文件。 【实验目的】 1、掌握循环结构程序设计的基本思想。 2、理解循环结构程序段中语句的执行过程。 3、了解goto语句和if语句构成循环的设计方法。 4、掌握用for、while、do-while语句实现循环（一重循环和多重（重点是双重）循环）） 程序设计。 5、掌握for、while、do-while三种循环的区别和联系及它们之间的灵活转换。 6、在程序设计过程中实现一些常用算法。 7、进一步练习程序的跟踪调试技术，掌握运行到光标处的调试方法。 8、巩固利用VC++对C程序的查错方法，进一步提高修改程序错误的能力。 【实验内容】 一、程序运行到光标处的跟踪调试和变量值的监视练习（可不上 交、但必须操作、理解，掌握其方法） 1、程序运行到光标处的跟踪调试和变量值的监视示例(请仔细按以下各步骤进行操作)： 编程计算表达式“1 + 2 + 3 + ……+ 100”的值。 源程序（有错误的程序） #include
运行结果（改正后程序的运行结果） sum = 5050 （1）输入源程序，并以error5_1.c 文件名保存，编译程序，出现第一个错误信息： missing ';' before ')' 双击该错误信息，箭头指向“for ”这一行，错误信息指出在for 语句的括号里面应使用“;”，对for 语句括号里面进行仔细分析，发现把“;”写成了“,”。把“,”改为“;”后，重新编译，连接，都正确。 （2）开始调试，如果调试工具栏不可见，采用实验四或实验二介绍的方法，调出该工具栏（见实验四图4.2所示）。 （3）鼠标单击第6行，光标就在第六行前面闪烁，这就是当前的光标位置（如图5.1所示）。 图5.1 光标在程序中位置 （4）单击 （Run to Cursor(Ctrl+F10)），程序运行到光标的位置（如图5.3所示）。 此过程中，出现图5.2 “生成可执行文件”的对话框时，选“是(Y)”按钮。 图5.2 生成“error5_1.exe ”文件 在变量窗口中，第一次循环时i 的值为1，正确，而sum 的值是-858993460，不正确。仔细分析程序，发现sum 没有赋初值。在for 语句前面加一条语句sum = 0; ，重新编译、连 光标位置
南邮MATLAB数学实验答案(全)
第一次练习 教学要求：熟练掌握Matlab 软件的基本命令和操作，会作二维、三维几何图形，能够用Matlab 软件解决微积分、线性代数与解析几何中的计算问题。 补充命令 vpa(x,n) 显示x 的n 位有效数字，教材102页 fplot(‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形 在下面的题目中m 为你的学号的后3位（1-9班）或4位（10班以上） 1.1 计算30sin lim x mx mx x →-与3 sin lim x mx mx x →∞- syms x limit((902*x-sin(902*x))/x^3) ans = 366935404/3 limit((902*x-sin(902*x))/x^3,inf) ans = 0 1.2 cos 1000 x mx y e =，求''y syms x diff(exp(x)*cos(902*x/1000),2) ans = (46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算 22 11 00 x y e dxdy +?? dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1) ans = 2.1394 1.4 计算4 2 2 4x dx m x +? syms x int(x^4/(902^2+4*x^2)) ans = (91733851*atan(x/451))/4 - (203401*x)/4 + x^3/12 1.5 (10)cos ,x y e mx y =求 syms x diff(exp(x)*cos(902*x),10) ans = -356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x) 1.6 0x =的泰勒展式(最高次幂为4).
MATLAB数学实验报告
Matlab 数学实验报告
一、实验目的 通过以下四组实验,熟悉MATLAB的编程技巧,学会运用MATLAB的一些主要功能、命令，通过建立数学模型解决理论或实际问题。了解诸如分岔、混沌等概念、学会建立Malthu模型和Logistic 模型、懂得最小二乘法、线性规划等基本思想。 二、实验内容 2.1实验题目一 2.1.1实验问题 Feigenbaum曾对超越函数y=λsin(πｘ)（λ为非负实数）进行了分岔与混沌的研究，试进行迭代格式x k+1=λsin(πx k),做出相应的Feigenbaum图 2.1.2程序设计 clear;clf; axis([0,4,0,4]); hold on for r=0:0.3:3.9 x=[0.1]; for i=2:150 x(i)=r*sin(3.14*x(i-1)); end pause(0.5) for i=101:150
plot(r,x(i),'k.'); end text(r-0.1,max(x(101:150))+0.05,['\it{r}=',num2str(r)]) end 加密迭代后 clear;clf; axis([0,4,0,4]); hold on for r=0:0.005:3.9 x=[0.1];
for i=2:150 x(i)=r*sin(3.14*x(i-1)); end pause(0.1) for i=101:150 plot(r,x(i),'k.'); end end 运行后得到Feigenbaum图
2.2实验题目二 2.2.1实验问题 某农夫有一个半径10米的圆形牛栏，长满了草。他要将一头牛拴在牛栏边界的桩栏上，但只让牛吃到一半草，问拴牛鼻子的绳子应为多长？ 2.2.2问题分析 如图所示，E为圆ABD的圆心，AB为拴牛的绳子，圆ABD为草场，区域ABCD为牛能到达的区域。问题要求区域ABCD等于圆ABC 的一半，可以设BC等于x,只要求出∠a和∠b就能求出所求面积。先计算扇形ABCD的面积，2a÷π×πx2=2aπ2，再求AB的面积，用扇形ABE的面积减去三角形ABE的面积即可。
第二次数学实验报告Matlab 二维曲线绘图
《数学实验》报告实验名称 Matlab 二维曲线绘图 2011年 5月
一、【实验目的】 学习Matlab 绘图的运用，学会制作二维曲线，三维图形的绘画。 二、【实验任务】 P79 第3，5，9题。 1，在同一图形窗口画三个子图…… 2，绘制圆锥螺线的图像并加各种标注…… 3，画三维曲面z=5-x^2-y^2与平面z=3的交线。 三、【实验程序】 1. >> clear >> x=-pi:pi/50:4*pi; y1=x.*cos(x); y2=x.*tan(1./x).*sin(x.^3); y3=exp(1./x).*sin(x); subplot(3,1,1) plot(x,y1,'r*'),grid on title('y1=xcosx') xlabel('x轴'),ylabel('y轴') axis([-pi pi -pi pi]) gtext('y1=xcosx'),legend('y1=xcosx') subplot(3,1,2),plot(x,y2,'b'),grid on title('y=xtan(1/x)sin(x^3)') gtext('y=xtan(1/x)sin(x^3)') legend('y=xtan(1/x)sin(x^3)') axis([pi 4*pi -2 2]) subplot(3,1,3),plot(x,y3,'y'),grid on title('y=exp(1/x)sinx') xlabel('x轴'),ylabel('y轴') gtext('y=exp(1/x)sinx') legend('y=exp(1/x)sinx') axis([1 8 -3 3]) 2. >> clear >> t=0:pi/50:20*pi; x=t.*cos(pi/6.*t); y=t.*sin(pi/6.*t); z=2.*t; plot3(x,y,z) title('圆锥螺线') xlabel('x轴'),ylabel('y轴'),zlabel('z轴') >> t=0:pi/50:20*pi; x=t.*cos(pi/6.*t);
实验五 循环结构设计
实验五循环结构程序设计 【目的与要求】 熟练掌握while、do—while和for三种循环语句的应用。 【上机内容】 【示例】只要求运行前4个例题，并在/* */中说明该语句的作用或使用的方法等；典型例题仅供参考，不作要求。 【例1】求5！ main() { int n,t; n=1; t=1; while(t<=5) { n=n*t; t=t+1;/**/ } printf(“%d”,n); } 【例2】求和s=1!+2!+3! main() { int n,s=0,t=1; for(n=1;n<=3;n++) { t=t*n;/*计算阶乘*/ s=s+t;/* */ } printf(“%d”,s); } 【例3】求和s= 3+33+333 main() { int s=0,t=3,i; for(i=1;i<=3;i++) { s=s+t; t=10*t+3;/*生成3,33,333等数据*/ } printf(“%d”,s); } 【例4】打印500以内的“水仙花数”。“水仙花数”是一个三位数，其各位数立方和等于该数本身。main()
int i,j,k,n; for(n=100;n<=500;n++) { i=n/100;/* */ j=n/10-i*10; /* */ k=n%10; /* */ if(n= =i*i*i+j*j*j+k*k*k) /* */ printf(“ %d”,n); } } 【典型示例】 【例1】计算两个数的最大公约数；分析：首先，随机输入两个数m，n(默认m>n)；其次，算法：使k为m除以n的余数，如果m能被n整除，则k值为0，n为这两个数的最大公约数，否则，使k代替n，n代替m，重复以上过程，直到k值为0。 #include"stdio.h" main() { int m,n,k,result; printf("Enter two numbers:"); scanf("%d,%d",&m,&n); if(m>0&&n>0) /*限定两个正整数*/ { do { k= m % n; if(k= =0) result=n; else { m = n; n=k; } }while(k>0); /*循环取余求出最大公因子*/ printf("The greatest common divistor is:%d\n",result); } else printf("Nonpositive values not allowed\n"); } 【例2】打印下列图案 * *** ***** ******* main()
数学实验与数学软件(Mathmaticandmatlab)
数学软件与数学实验2013-2014学年度秋季学期期末试卷 专业：统计学 班级：11级2班 学号：20110723 姓名：晏静
一、按要求计算出下列表达式的值 (1)318, 3 162 53 ?? + ? ?? , 21 eπ+, 2.5 tg, 2 log15； (2)给出π的9位和e的10位近似值； (3)求658和4102的最大公约数及35和25的最小公倍数； (4)产生10个0与10之间随机数的一个表; (5)求虚数1453 i i i i +- -的实部,虚部,模,共轭,辐角。 (6)自己运用Table建立两个表，并进行表运算，如连接、并集、交、排序等操作。
二、因式分解 22212321332112322 1 22(1)()()()4;(2)21;x x x x x x x x x x x x x x x +++++---- 解： 三、解方程(组) 1234234124234-2+344-+-3(1)+31-73+3 x x x x x x x x x x x x x -=??=? ? +=??+=-? 65432(2)5232002000.x x x x x x -+--++= 四、求极限 () 20 (1)1sin ;(2);(3)56! ctg x n x n n n Lim x Lim n n →→∞ →∞++
(1) (2) (3) 五、求导数 32 22(1)()=ln(x+1+);(2)()=cos 2,; (3)=log (),Z . x f x x f f x e y x y Z xy x y y ???求的导数已知求求关于的二阶导 (1) (2) (3) 六、求下列定积分与不定积分： ()()()12201+sin ln 1+(1);(2);(3)sin (1+cos ) +1(1+)(2+-) x x dx dx x x x x x x ? ? ?2 2-(4)=0,=1,==.y D D x y y x I x e d σ??设是由直线围成的区域，计算的值 (1) (2)
matlab实验内容答案
实验报告说明： matlab 课程实验需撰写8个实验报告，每个实验报告内容写每次实验内容中标号呈黑体大号字显示的题目。 第一次实验内容： 实验一 MATLAB 运算基础 一、实验目的 1.熟悉启动和退出MA TLAB 的方法。 2.熟悉MA TLAB 命令窗口的组成。 3.掌握建立矩阵的方法。 4.掌握MA TLAB 各种表达式的书写规则以及常用函数的使用。 二、实验内容 1.先求下列表达式的值，然后显示MA TLAB 工作空间的使用情况并保存全部变量。 （1）2 2sin 8511z e ?= + （2 ）12ln(2 z x =+ ，其中2120.45 5i +? ? =? ?-?? （3）0.30.33sin(0.3), 3.0, 2.9, 2.8,,2.8,2.9,3.02 a a e e z a a --= +=--- 提示：利用冒号表达式生成a 向量，求各点的函数值时用点乘运算。 （4）2 2 2 01 41 1221 23 t t z t t t t t ?≤
12344347873657A -????=??????，131203327B -???? =????-?? 求下列表达式的值： （1）A+6=B 和A-B+I(其中I 为单位矩阵)。 （2）A*B 和A.*B 。 （3）A^3和A^.3 。 （4）A/B 和B\A 。 （5）[A ，B]和[A([1,3],;);B^2] 。 3.设有矩阵A 和B 12345678910111213141516171819202122232425A ????????=????????， 30 161769 23497041311B ?? ?? -?? ??=-?? ????? ? （1） 求它们的乘积C 。 （2） 将矩阵C 的右下角3×2子矩阵赋给D （3） 查看MA TLAB 工作空间使用情况。 4.完成下列操作： （1）求[100，999]之间能被21整除的数的个数。 提示：先利用冒号表达式，再利用find 和length 函数。 （2）建立一个字符串向量，删除其中的大写字母。 提示：利用find 函数和空矩阵。 第二次实验内容： 实验三 选择结构程序设计 一、实验目的 1. 掌握建立和执行M 文件的方法。 2. 掌握利用if 语句实现选择结构的方法。 3. 掌握利用switch 语句实现多分支选择结构的方法。 4. 掌握try 语句的使用。 二 、实验内容
- 《数学实验与Matlab》程序
- MATLAB数学实验
- Matlab数学实验知识点与函数集
- Matlab数学实验一2015答案版
- Matlab数学实验报告4
- 北京科技大学《matlab和数学实验》绘图.ppt
- 数学实验与Matlab
- Matlab数学实验一2015(答案版)
- MATLAB数学实验答案(全)
- matlab数学实验练习题
- 数学实验matlab练习题
- matlab数学实验答案 (胡良剑)版本
- 北科大Matlab_数学实验报告1~6次(全)
- 南邮MATLAB数学实验答案(全)
- matlab 数学实验 实验报告 欧拉公式 ROSSLER微分方程
- MATLAB数学实验报告
- matlab数学实验体会
- 数学实验和Matlab介绍.ppt
- Matlab数学实验一2015(标准答案版)
- matlab数学实验