3．计算
∑=10
1
!k k ，要求：
（1）用循环结构并使用上体的M 函数文件 ①s=0; for k=1:10 f=f1(k); s=s+f; end s
结果：s =
4500244
②function f=f4(n) f=0;
for i=1:n
f=f+f1(i); end
（2）用help 命令查询sum 和prod 的用法，并用这两个指令完成计算。
n=input('enter a number:'); sum=0; for i=1:n
sum=sum+prod([1:i]); end sum
4．用尽可能多的方法编写Matlab 函数mysum1.m 来计算和式f(n)=
∑=+n
i i i 1
)*/11ln(.其中，n 为函数的
输入变量，函数输出为f （n ）的值，分析指出他们执行效率的高低，给出理由。 ①n=input('enter a number:'); s=0;
for i=1:n
t=log(1+1/(i*i)); s=s+t;
end s
②function f=mysum1(n) f=0; k=0;
while k f=f+log(1+1/(k*k)); end ③function f=mysum11(n) if n==1 f=log(2); elseif n>1 f=mysum11(n-1)+log(1+1/(n*n)); else error('wrong input parameter!'); end ④function f=mysum12(n) k=1:n; s=log(1+1./(k.*k)); f=sum(s); 5．、用尽可能多的方法编写Matlab 函数mysum1.m 来计算如下和式： f(n)=1/(1*4)+1/(4*7)+…+1/((3n-2)*（3n+1）) 其中，n 为函数的输入变量，函数输出为f(n)的值。 ①function f=mysum2(n) f=0; for k=1:n t=1/((3*k-2)*(3*k+1)); f=f+t; end ②function f=mysum21(n) f=0; k=0; 数学软件与数学实验作业 while k k=k+1; t=(3*k-2)*(3*k+1); f=f+1/t; end ③function f=mysum22(n) if n==1 f=1/4; elseif n>1 t=(3*n-2)*(3*n+1); f=mysum22(n-1)+1/t; else error('wrong input parameter!'); end ④function f=mysum23(n) k=1:n; t=(3.*k-2).*(3.*k+1); f=sum(1./t); 6．写一个程序用于得到阶乘不超过10^100的最小整数。 for k=1:5000000000000 f=1; for i=1:k f=f*i; end if(f<=10^100) continue; end break; end k 结果：k = 70 7．写一个函数rs=f(s),对输入放入字符串变量s,删除其中的小写字母，然后将原来的大写字母变为小写字母，得到rs的返回 function rs=fy(s) k=length(s); p=[]; for i=1:k if s(i)>='a'&s(i)<='z' p=[p,i]; elseif s(i)>='A'&s(i)<='Z' s(i)=char(s(i)+32); end end s(p)=[]; rs=s; 8. 设计函数min_element.m,其功能是在一个二维矩阵中找出其最小元素，函数定义如下： [minEle,row,column]=min_element(mat rix) 其中，matrix是一个二维矩阵，minEle是返回的最小元素的值，row是其所在行，coulmn是其所在列。要求： （1）在调用函数时，要写出调用的语句； （2）不适用Matlab函数min编写；(3)使用min函数和find函数编写。function [min,row,column]=min_element(x) n=size(x); min=x(1); for i=1:n(1) for j=1:n(2) if x(i,j) min=x(i,j); [row,column]=find(x<=min); end end end 9、若一个自然数是素数，且它的个位数 字位置经过任意的对换之后仍为素数，则 该数为绝对素数，例如，113是绝对素数，试求所有三位数的绝对素数。 function [x1,x2,x3,x4,x5,x6]=f(x) a=rem(x,10); b=rem(fix(x/10),10); c=fix(x/100); x1=x; x2=c*100+a*10+b; x3=a*100+c*10+b; x4=a*100+b*10+c; x5=b*100+a*10+c; x6=b*100+c*10+a; for i=100:999 [x1,x2,x3,x4,x5,x6]=f(i); if isprime(x1)~=0&isprime(x2)~=0&is prime(x3)~=0&isprime(x4)~=0&ispr ime(x5)~=0&isprime(x6)~=0 i end end 运行结果 i = 113 i = 131 i = 199 i = 311 i = 337 i = 373 i = 733 i = 919 i = 991 北京交通大学海滨学院考试试题 课程名称：数学实验2010-2011第一学期出题教师：数学组适用专业: 09机械, 物流, 土木, 自动化 班级：学号：姓名： 选做题目序号： 1.一对刚出生的幼兔经过一个月可以长成成兔, 成兔再经过一个月后可以 繁殖出一对幼兔. 如果不计算兔子的死亡数, 请用Matlab程序给出在未来24个月中每个月的兔子对数。 解: 由题意每月的成兔与幼兔的数量如下表所示： 1 2 3 4 5 6 ··· 成兔0 1 1 2 3 5··· 幼兔 1 0 1 1 2 3··· 运用Matlab程序： x=zeros(1,24); x(1)=1;x(2)=1； for i=2:24 x(i+1)=x(i)+x(i-1); end x 结果为x = 1 1 2 3 5 8 13 21 3 4 5 5 89 144 233 377 610 987 1597 2584 4181 6765 1094 6 7711 2865 7 46368 2.定积分的过程可以分为分割、求和、取极限三部分, 以1 x e dx 为例, 利用 已学过的Matlab 命令, 通过作图演示计算积分的过程, 并与使用命令int() 直接积分的结果进行比较. 解：根据求积分的过程，我们先对区间[0,1]进行n 等分， 然后针对函数x e 取和，取和的形式为10 1 i n x i e e dx n ξ=≈ ∑ ? ，其中1[ ,]i i i n n ξ-?。这里取i ξ为区间的右端点，则当10n =时，1 x e dx ?可用10 101 1.805610 i i e ==∑ 来近似计算， 当10n =0时，100 100 1 01 =1.7269100 i x i e e dx =≈ ∑?，当10n =000时，10000 10000 1 1 =1.718410000 i x i e e dx =≈ ∑ ?. 示意图如下图，Matlab 命令如下： x=linspace (0,1,21); y=exp(x); y1=y(1:20); s1=sum(y1)/20 y2=y(2:21); s2=sum(y2)/20 plot(x,y); hold on for i=1:20 fill([x(i),x(i+1),x(i+1),x(i),x(i)],[0,0,y(i),y(i),0],'b') end syms k;symsum(exp(k/10)/10,k,1,10)；%n=10 symsum(exp(k/100)/100,k,1,100)；%n=100 symsum(exp(k/10000)/10000,k,1,10000)；%n=10000 数学软件与数学实验作业 一.《数学软件》练习题（任选12题，其中19－24题至少选2题）： 3.对下列各式进行因式分解. (1). syms x y >> factor(x^5-x^3) (2). syms x y >> factor(x^4-y^4) (3). syms x >> factor(16-x^4) (4). syms x >> factor(x^3-6*x^2+11*x-6) (5). syms x y >> factor((x+y)^2-10*(x+y)+25) (6). syms x y >> factor(x^2/4+x*y+y^2) (7). syms x y a b >> factor(3*a*x+4*b*y+4*a*y+3*b*x) (8). syms x >> factor(x^4+4*x^3-19*x^2-46*x+120) 5.解下列方程或方程组. (1).solve('(y-3)^2-(y+3)^3=9*y*(1-2*y)') (2). solve('3*x^2+5*(2*x+1)') (3). solve('a*b*x^2+(a^4+b^4)*x+a^3*b^3','x') (4). solve('x^2-(2*m+1)*x+m^2+m','x') (5). [x,y]=solve('4*x^2-9*y^2=15','2*x-3*y=15') 6.计算极限. (1). syms x f=(exp(x)-exp(-x))/sin(x); limit(f,x,0) (2) syms x >> f=(x/(x-1)-1/log(x)); >> limit(f,x,1) (3). syms x >> f=(1-cos(x))/x^2; >> limit(f,x,0) 一、实验内容 P206第六题 function f=wuyan2(c) y=[3.9 5.3 7.2 9.6 12.9 17.1 23.2 31.41 38.6 50.2 62.9 76.0 92.0 106.5 123.2 131.7 150.7 179.3 204.0 226.5 251.4 281.4] t=[0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200 210] f=y-c(1)/(1+c(1)/3.9-1)*exp^(-c(2)*t) c0=[1 1] c=lsqnonlin('wuyan2',c0) P206第七题 function f=wuyan1(c) q=[0.4518 0.4862 0.5295 0.5934 0.7171 0.8964 1.0202 1.1963 1.4928 1.6909 1.8548 2.1618 2.6638 3.4634 4.6759 5.8478 6.7885 7.4463 7.8345 8.2068 8.9468 9.7315 10.5172 11.7390 13.6876 ]; k=[0.0911 0.0961 0.1230 0.1430 0.1860 0.2543 0.3121 0.3792 0.4754 0.4410 0.4517 0.5595 0.8080 1.3072 1.7042 2.0019 2.2914 2.4941 2.8406 2.9855 3.2918 3.7214 4.3500 5.5567 7.0477]; l=[4.2361 4.3725 4.5295 4.6436 4.8179 4.9873 5.1282 5.2783 5.4334 5.5329 6.4749 6.5491 6.6152 6.6808 6.7455 6.8065 6.8950 6.9820 7.0637 7.1394 7.2085 7.3025 7.3470 7.4432 7.5200]; f=q-c(1)*k.^c(2).*l.^c(3) c0=[1 1 1] c=lsqnonlin('wuyan1',c0) c = 0.4091 0.6401 1.1446 a=0.4091 α=0.6401 β=1.1446 P239第五题 c=[-20 -30]; A=[1 2;5 4]; b=[20 70]; v1=[0 0]; [x,f,ef,out,lag]=linprog(c,A,b,[],[],v1) z=-f x = 10.0000 5.0000 实验报告 课程名称：数学软件姓名： 学院： 专业： 年级： 学号： 指导教师： 职称： 年月日 实验项目列表 附件三： 实验报告（二） 系：专业：年级：姓名学号：实验课程： 实验室号：_ 实验设备号：实验时间： 指导教师签字：成绩： 1. 实验项目名称：符号计算基础与符号微积分 2. 实验目的和要求 1.掌握定义符号对象的方法 2.掌握符号表达式的运算法则以及符号矩阵运算 3.掌握求符号函数极限及其导数的方法 4.掌握求符号函数定积分和不定积分的方法 3. 实验使用的主要仪器设备和软件 方正商祺N260微机；MATLAB7. 0或以上版本 4. 实验的基本理论和方法 （1）符号函数;sym(x)；syms a b …… （2）平方根：sqrt(x) （3）分解因式：factor（s） （4）符号表达式化简：simplify（s） （5）逆矩阵：inv(x) （6）下三角矩阵：tril（x） （7）矩阵行列式的值:det(x) （8）符号函数求极限：limit （f ，x ，a ）；limit （f ，x ，a ，‘right ’） （9）符号函数求导：diff （f ，v ，n ） （10）符号函数求不定积分：int （f ，v ） （11）符号函数求定积分：int （f ，v ，a ，b ） 5. 实验内容与步骤 （描述实验中应该做什么事情，如何做等，实验过程中记录发生的现象、中间结果、最终得到的结果，并进行分析说明） （包括：题目，写过程、答案） 题目： 1. 已知x=6，y=5，利用符号表达式求 y x x z -++= 31。 提示：定义符号常数)'5(')'6('sym y sym x ==，。 >> x=sym('6'); >> y=sym('5'); >> z=(x+1)/(sqrt(3+x)-sqrt(y)) z = 7/(3-5^(1/2)) 2. 分解因式：44y x - >> syms x y; >> A=x^4-y^4; >> factor(A) ans = (x-y)*(x+y)*(x^2+y^2) 3. 化简表达式 (1)2121sin cos cos sin ββββ- (2) 123842+++x x x (1) >> syms x y; >> f1=sin(x)*cos(y)-cos(x)*sin(y); 实验2 方程模型及其求解算法 一、实验目的及意义 [1] 复习求解方程及方程组的基本原理和方法； [2] 掌握迭代算法； [3] 熟悉MATLAB软件编程环境；掌握MATLAB编程语句(特别是循环、条件、控制等语句)； [4] 通过范例展现求解实际问题的初步建模过程； 通过该实验的学习，复习和归纳方程求解或方程组求解的各种数值解法（简单迭代法、二分法、牛顿法、割线法等），初步了解数学建模过程。这对于学生深入理解数学概念，掌握数学的思维方法，熟悉处理大量的工程计算问题的方法具有十分重要的意义。 二、实验内容 1．方程求解和方程组的各种数值解法练习 2．直接使用MATLAB命令对方程和方程组进行求解练习 3．针对实际问题，试建立数学模型，并求解。 三、实验步骤 1．开启软件平台——MATLAB，开启MATLAB编辑窗口； 2．根据各种数值解法步骤编写M文件 3．保存文件并运行； 4．观察运行结果(数值或图形)； 5．根据观察到的结果写出实验报告，并浅谈学习心得体会。 四、实验要求与任务 基础实验 1．用图形放大法求解方程x sin(x) = 1. 并观察该方程有多少个根。 画出图形程序： x=-10:0.01:10; y=x.*sin(x)-1; y1=zeros(size(x)); plot(x,y,x,y1) MATLAB运行结果： -10-8-6-4-20246810 -8-6 -4 -2 2 4 6 8 扩大区间画图程序： x=-50:0.01:50; y=x.*sin(x)-1; y1=zeros(size(x)); plot(x,y,x,y1) MATLAB 运行结果： -50-40-30-20-1001020304050 由上图可知，该方程有偶数个无数的根。 数学与软件科学学院实验报告 一、实验目的 (1) 掌握C语言环境下结构体和共用体类型变量的定义和使用方法； (2) 掌握结构体类型数组的概念和使用； (3) 掌握指向结构体变量的指针变量、尤其是链表概念； 二、实验内容 1.首先熟悉结构体类型变量的基本声明方法、结构体类型变量的内存分配原则、初始化和引用结构体变量及其成员变量的基本方法；然后掌握结构体变量的输入、输出方法。(参见教材例7.1，请给该例加上输入功能) #include for(i=0;i<3;i++) printf("%5s:%d\n",leader[i].name,leader[i].count); } 2.基于结构体数组的应用实验。 (1) 有n个学生，每个学生的数据包括学好(num)、姓名(name[20])、性别(sex)、年龄(age)，以及三门课程的成绩(score[3])。要求：在main()函数中输入这些学生的这些数据，然后设计一个函数count()来计算每个学生的总分和平均分，最后， 打印出所有数据信息(包含原来输入的学生原始数据信息和求解出来的新信息)。#include 《数学建模与数学实验》上机报告（第 1 次） 一、上机训练目的、题目或内容（简述综述）等 题目一：数学软件（MathType5.2、MATLAB 、Maple、Mathematica4.0、LINGO8.0）安装调试；基本命令使用（变量赋值、定义函数、过程控制、绘图命令、拟合、线性规划、非线性规划）；高等数学实验（绘图，极限，求导，积分，解微分方程）；线性代数实验（矩阵基本运算，线性方程组求解，解超定方程组，优化命令）。调试运行给定的两个程序： 题目二： 1、以两种方式打开MATLAB 工作窗口，进入MATLAB 6.0 的工作环境，并尝试用不同的方式退出。（这个在报告里面说明方法就可以） 2、尝试、熟悉MATLAB 6.0 的各栏菜单以及各个工具栏的功能。（自己掌握，报告里面就不写了） 3、绘制函数y=cos(5x+2)/sin(3x+1) 的图像，并求解当x=2 时的函数值。 4、练习并熟练掌握MATLAB 的帮助命令，学会利用MATLAB 的帮助信息。 5、求矩阵Ａ=的行列式、逆的特征根；Ｂ=，解方程ＢＸ= 6、两个矩阵Ａ=Ｂ=将矩阵改为3行3列的矩阵，作加、减、乘和除（左 除，右除）运算，同事运用数组运算法则进行运算，比较二者计算结果有何异同。 二、数学模型或求解分析或算法描述程序命令图形等 题目一： 1) c=[6,3,4]; A=[0,1,0]; b=[50]; Aeq=[1,1,1]; beq=[120]; vlb=[30,0,20]; vub=[]; [x,fval]=linprog(c,A,b,Aeq,beq,vlb,vub) 2) function f=fun3(x); f=-x(1)-2*x(2)+(1/2)*x(1)^2+(1/2)*x(2)^2 x0=[1;1]; A=[2 3 ;1 4]; b=[6;5]; Aeq=[];beq=[]; VLB=[0;0]; VUB=[]; [x,fval]=fmincon('fun3',x0,A,b,Aeq,beq,VLB,VUB) 题目二： 3. x=2; y=cos(5*x+2)./sin(3*x+1) x=[-10:0.01:10]; y=cos(5*x+2)./sin(3*x+1); plot(x,y) 数学软件实验报告 学院名称：理学院专业年级： 姓名：学号： 课程：数学软件实验报告日期：2014年12月6日 实验七SIMULINK建模与工具箱的使用 一．实验目的 MATLAB 具有丰富的可用于各种专业方向的工具箱，这些工具箱已经形成了MATLAB 的系列产品。特别是动态仿真建模工具箱，更是成为许多工具箱的基础。本次实验的目的就是要使大家了解MA TLAB工具箱使用的基本方法，以及如何查询工具箱，主要掌握系统优化工具箱的使用和系统动态仿真建模工具箱的使用。 二．实验要求 MATLAB系统的工具箱十分的丰富，并且随着版本的不断升级，其工具箱还在不断地增加。通过本次实验，要求了解MA TLAB系统工具箱的分类与查询，会使用系统优化工具箱解决一些实际问题。能建立系统仿真方框图，并进行系统仿真模拟。 三．实验内容 最优化工具箱 非线性最小化函数 fgoalattain 多目标达到优化 constr 有约束最小化 fminbnd 有边界最小化 fminunc使用梯度法的无约束最小化 fminsearch 使用简单法的无约束最小化 fzero 非线性方程求解（数量情况） fsolve 非线性方程求解 lsqnonlin 非线性最小二乘 fminimax 最小的最大解 fseminf 半无穷区间最小化 2.矩阵问题的最小化 linprog 线性规划 quadprog 二次规划 lsqnonneg 非负线性最小二乘 lsqlin 约束线性最小二乘 第十章 10.1线性优化 >> f=[-5 4 2]; >> a=[6 -1 1;1 2 4]; >> b=[8 10]; >> 1b=[-1 0 0]; >> ib=[-1 0 0]; >> ub=[3 2]; >> [x,fval,exitflag,output,lambda]=linprog(f,a,b,[],[],ib,ub) Optimization terminated. x = 1.3333 0.0000 0.0000 fval = -6.6667 exitflag = 1 output = iterations: 7 algorithm: 'large-scale: interior point' cgiterations: 0 message: 'Optimization terminated.' constrviolation: 0 lambda = ineqlin: [2x1 double] 东华大学M A T L A B数学实验第二版答案(胡良 剑) -CAL-FENGHAI-(2020YEAR-YICAI)_JINGBIAN 数学实验答案 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 = 实验一：了解数学软件MATLAB 实验目的与要求： 了解MATLAB的基本特点，掌握MATLAB界面上的主要窗口，熟练掌握MATLAB的帮助系统。 实验内容： 1启动按钮 打开MATLAB主界面以后，单击“Start”按钮，显示一个菜单，利用“START”菜单机器子菜单中的选项，打开MATLAB的有关工具。 2命令窗口 命令窗口（Command Window）是用于输入数据，运行MA TLAB函数和脚本并显示结果的主要工具之一。命令窗口没有打开时，从“Desktop”菜单中选择“Command Window”选项可以打开它。 在命令窗口中键入命令并执行：a=[123;456;789] 在上述语句末尾加分号“；”，结果是什么？请与不加分号的情况作比较。 功能。 命令历史窗口（command history）显示命令窗口中最近输入的所有语句。先关闭历史窗口，再分别用“Desktop”菜单打开它和用command history命令打开它。 （1）将命令历史窗口中的语句复制到命令窗口中； （2）直接双击命令历史窗口中的语句。 4工作空间窗口 清空工作空间的命令是：clear 清空命令窗口的命令是：clc 在命令窗口中键入： t=0:pi/4:2*pi y=sin(t) 在命令窗口中键入：who,看运行结果； 在命令窗口中键入：whos,看运行结果； 在命令窗口中键入：whos y,看运行结果。 退出MATLAB时，工作空间中的内容随之清除。可以将当前空间中的部分或全部变量保存到一个MA T文件中，它是一个二进制文件，扩展名为.mat。然后可以在以后使用它时载入它，请练习这一操作过程。 5帮助浏览器 使用帮助浏览器可以搜索和查询所有Math Works产品的文档和演示。帮助浏览器是集成到MA TLAB桌面的一个HTML查看器。请打开帮助浏览器，熟悉它。 分别用Help函数和doc函数获取format函数的帮助，进而说明format函数的功能 数学软件与实验第一次上机作业 上机时间：2013-4-10 地点：E204 班级:071111 学号:07111014 姓名:曹红兴xdhjtang@https://www.doczj.com/doc/e99602891.html, 学号、姓名、MATLAB、第一次作业 1.计算三角形三边分别为a，b，c中c边对应内角的角度 >> a = 3; b = 3; c = 3; >> acos((a^2+b^2-c^2)/2/a/b) *180/pi ans = 60.0000 >> a = 3; b = 4; c = 5; >> acos((a^2+b^2-c^2)/2/a/b) *180/pi ans = 90 >> a = 3; b = 4; c = 20; >> acos((a^2+b^2-c^2)/2/a/b) *180/pi ans = 1.8000e+002 -1.9715e+002i 2.试分别生成5 阶的单位阵、8 阶均匀分布的随机矩阵及其下三角 矩阵，要求矩阵元素为介于10~99之间整数 >> C=eye(5,5) C = 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 >> N=randsrc(8,8,[10:99]) N = 59 21 72 34 19 76 25 52 66 28 54 28 22 15 45 68 12 23 58 60 24 87 22 12 65 27 50 67 65 94 12 85 42 13 21 47 61 98 94 60 14 67 54 28 14 87 37 86 54 35 86 95 93 80 36 41 27 58 88 17 75 56 39 50 >> Z=tril(N) Z = 83 0 0 0 0 0 0 0 91 96 0 0 0 0 0 0 21 24 81 0 0 0 0 0 92 97 96 45 0 0 0 0 66 96 69 68 72 0 0 0 18 53 13 25 38 54 0 0 35 82 86 73 95 50 20 0 59 22 94 12 13 68 54 72 3.生产列向量x=[1, 3, 10, 15, 20, 25, 30, 35, 40] >> x=[1;3;10;15;20;25;30;35;40] x = 1 3 10 15 20 25 30 35 40 重庆科技学院 数学实验与数学软件课程设计 课程名称：菜单与对话框设计 开课学期：_2014-2015-1 学院:__ 数理学院 开课实验室：_数学实验与建模实验室_ 学生姓名: 谭云文 专业班级: 应数13-2班 __ 学号:___ 20134432214 _ 实验十二 我们本次实验做的是菜单与对话框设计，所谓菜单与对话框的设计包括在图形用户界面中。而图形用户界面是由窗口、菜单、对话框等各种图形元素组成的用户界面。因为在这种用户界面中，用户的操作既生动形象，又方便灵活，这是它的一大特点。 在MATLAB中，基本的图形用户界面对象包含3类：用户界面控件对象、下 拉式菜单对象和快捷菜单对象，可以设计出界面友好、操作方便的图形用户界面。 其中MATLAB用户菜单对象是图形窗口的子对象，所以菜单设计总在某一个图形 窗口中进行。MATLAB的图形窗口有自己的菜单栏。为了建立用户自己的菜单系 统，可以先将图形窗口的MenuBar属性设置为none，以取消图形窗口默认的菜 单，然后再建立用户自己的菜单。对话框是用户与计算机进行信息交流的临时窗 口，在现代软件中有着广泛的应用。在软件设计时，借助于对话框可以更好地满 足用户操作需要，使用户操作更加方面灵活。为了更便捷地进行用户界面设计， MATLAB提供了图形用户界面开发环境，这使得界面设计在可视化状态进行，设计过程中变得简单直观，实现了“所见即所得”。 例1 一、实验目的 1. 掌握plot菜单的方法。 2. 掌握建立控件对象的方法。 3. 掌握对话框设计的方法。 二、实验内容 设计图1所示的菜单。 菜单条上仅有Plot菜单，其中有Sine Wave、Cosine Wave和Exit共3个命令。若选择了其中的Sine Wave命令，则将绘制出正弦曲线；若选择了其中的Cosine Wave命令，则将绘制出余弦曲线；如果选择了Exit命令，则将关闭窗口。 程序如下： screen=get(0,'ScreenSize'); W=screen(3);H=screen(4); figure('Color',[1,1,1],'position',[0.2*H,0.2*H,0.5*W,0.3*H],... 'Name','图形演示系统','NumberTitle','off','Menubar','none'); %plot hplot=uimenu(gcf,'Label','&Plot'); uimenu(hplot,'Label','Sine Wave','Call',... ['t=-pi:pi/20:pi;','plot(t,sin(t));',... 'set(hgon,''Enable'',''on'');',... 'set(hgoff,''Enable'',''on'');',... 'set(hbon,''Enable'',''on'');',... 'set(hboff,''Enable'',''on'');']); uimenu(hplot,'Label','Cosine Wave','Call',... ['t=-pi:pi/20:pi;','plot(t,cos(t));',... 'set(hgon,''Enable'',''on'');',... 'set(hgoff,''Enable'',''on'');',... 'set(hbon,''Enable'',''on'');',... 'set(hboff,''Enable'',''on'');']); uimenu(hplot,'Label','&Exit','Call','close(gcf)'); 三、运行结果 1.点击SineWave函数将出现我们所需要的图像，如图： 2点击CosineWave函数将出现我们所需要的图像，如图： 二、判断题（“对”在题号前（）中打√×）（10分） （√）1、误差是指测量值与真值之差，即误差=测量值－真值，如此定义的误差反映的是测量值偏离真值的大小和方向，既有大小又有正负符号。 （×）2、残差（偏差）是指测量值与其算术平均值之差，它与误差定义一样。（√）3、精密度是指重复测量所得结果相互接近程度，反映的是随机误差大小的程度。 （√）4、测量不确定度是评价测量质量的一个重要指标，是指测量误差可能出现的范围。 （×）7、分光计设计了两个角游标是为了消除视差。 （×）9、调节气垫导轨水平时发现在滑块运动方向上不水平，应该先调节单脚螺钉再调节双脚螺钉。 （×）10、用一级千分尺测量某一长度（Δ仪=0.004mm），单次测量结果为N=8.000mm，用不确定度评定测量结果为N=（8.000±0.004）mm。 三、简答题（共15分） 1．示波器实验中，（1）CH1（x）输入信号频率为50Hz，CH2（y）输入信号频率为100Hz；（2）CH1（x）输入信号频率为150Hz，CH2（y）输入信号频率为50Hz；画出这两种情况下，示波器上显示的李萨如图形。（8分） 差法处理数据的优点是什么？（7分） 答：自变量应满足等间距变化的要求，且满足分组要求。（4分） 优点：充分利用数据；消除部分定值系统误差 四、计算题（20分，每题10分） 1、用1/50游标卡尺，测得某金属板的长和宽数据如下表所示，求金属板的面 解：（1）金属块长度平均值：)(02.10mm L = 长度不确定度： )(01.03/02.0mm u L == 金属块长度为：mm L 01.002.10±= %10.0=B （2分） （2）金属块宽度平均值：)(05.4mm d = 宽度不确定度： )(01.03/02.0mm u d == 金属块宽度是：mm d 01.005.4±= %20.0=B （2分） （3）面积最佳估计值：258.40mm d L S =?= 不确定度：2222222 221.0mm L d d s L s d L d L S =+=??? ????+??? ????=σσσσσ 相对百分误差：B ＝%100?S s σ=0.25％ （4分） （4）结果表达：21.06.40mm S ±= B ＝0.25％ （2分） 注：注意有效数字位数，有误者酌情扣 5、测量中的千分尺的零点误差属于已定系统误差；米尺刻度不均匀的误差属于未 第一次练习 教学要求：熟练掌握Matlab 软件的基本命令和操作，会作二维、三维几何图形，能够用Matlab 软件解决微积分、线性代数与解析几何中的计算问题。 补充命令 vpa(x,n) 显示x 的n 位有效数字，教材102页 fplot(‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形 在下面的题目中m 为你的学号的后3位（1-9班）或4位（10班以上） 计算30sin lim x mx mx x →-与3 sin lim x mx mx x →∞- 程序： syms x limit((1001*x-sin(1001*x))/x^3,x,0) 结果： 程序： syms x limit((1001*x-sin(1001*x))/x^3,x,inf) 结果： 0 cos 1000 x mx y e =，求''y 程序： syms x diff(exp(x)*cos(1001*x/1000),2) 结果： -2001/1000000*exp(x)*cos(1001/1000*x)-1001/500*exp(x)*sin(1001/1000*x) 计算 2 2 11 00 x y e dxdy +?? 程序： dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1) 结果： 计算4 2 2 4x dx m x +? 程序： syms x int(x^4/(1000^2+4*x^2)) 结果： (10)cos , x y e mx y =求 程序： syms x diff(exp(x)*cos(1000*x),10) 结果： 给出 0x =的泰勒展式(最高次幂为4). 程序： syms x taylor(sqrt(1001/1000+x),5) 结果： Fibonacci 数列{}n x 的定义是121,1x x ==， 12,(3,4,)n n n x x x n --=+=L 用循环语句编程给出该数列的前20项（要求将结果用向量的形式给出）。 程序： x=[1,1]; for n=3:20 x(n)=x(n-1)+x(n-2); end x 结果： Columns 1 through 10 1 1 2 3 5 8 13 21 3 4 5 5 Columns 11 through 20 89 144 233 377 610 987 1597 2584 4181 6765 《数学实验》第一次上机实验 1． 设有分块矩阵?? ? ???= ????22322333S O R E A ，其中E,R,O,S 分别为单位阵、随机阵、零阵和对角阵，试通过数值计算验证?? ????+= 22 S 0RS R E A 。 程序及结果： E=eye(3); %创建单位矩阵E% R=rand(3,2); %创建随机矩阵R% O=zeros(2,3); %创建0矩阵% S=diag(1:2); %创建对角矩阵% A=[E,R;O,S]; %创建A 矩阵% B=[E,(R+R*S);zeros(2,3),S^2] %计算等号右边的值% A^2 %计算等号左边的值% 运行结果： B = 1.00 0 0 1.63 2.74 0 1.00 0 1.81 1.90 0 0 1.00 0.25 0.29 0 0 0 1.00 0 0 0 0 0 4.00 ans = 1.00 0 0 1.63 2.74 0 1.00 0 1.81 1.90 0 0 1.00 0.25 0.29 0 0 0 1.00 0 0 0 0 0 4.00 2．某零售店有9种商品的单件进价（元）、售价（元）及一周的销量如表1.1，问哪种商品的利润最大，哪种商品的利润最小；按收入由小到大，列出所有商品及其收入；求这一周该10种商品的总收入和总利润。 表1.1 1)程序： a=[7.15 8.25 3.20 10.30 6.68 12.03 16.85 17.51 9.30]; b=[11.10 15.00 6.00 16.25 9.90 18.25 20.80 24.15 15.50]; c=[568 1205 753 580 395 2104 1538 810 694]; 数学软件与数学实验_15 林府标Tel : 668960 Email ： 三维爱的方程式 3 2222323 99 10 480 x y z x z y z ?? ++---= ? ?? 图像 一、单项选择题 1. 下列变量名中( )是合法的. A. P72_5_01 B. x*y C.x/y D. end 2.下列那些变量是合法的？（） A. fxjgxy-11 B. P_1_3 C. 函数g D. 7.x 3.下列那些变量是合法的？（） A. 999abcd B. 11 C. P.m D. A2 4.下列标量中比0大的最大数是（） A. exp(1) B. eps C. realmax D. realmin 5. 在MATLAB命令窗口输入 >> syms a b c d e f g y ; S=a+2*b+3*c+4*d+5*e*f*g*y; findsym(S,1) 则该命令运行结果是（） A. a+2*b+3*c+4*d+5*e*f*g*y B. 5*e*f*g*y C. a b c d e f g y D. y 6. 在MATLAB命令窗口输入 >> syms a b c d e f g y ; S=a+2*b+3*c+4*d+5*e*f*g*y; findsym(S) 则该命令运行结果是（） A. a+2*b+3*c+4*d+5*e*f*g*y B. 5*e*f*g*y C. a b c d e f g y D. y 7. 在循环结构中跳出循环，但继续执行循环语句的下一语句的命令是( )． A.return B.break C. continue D. keyboad 8. 在MATLAB中下列数值的表示不正确的是( ). A. log(2) B. log3(3) C. log10(10) D. log2(2) 9. 在MATLAB中下列数值的表示不正确的是( ). A. 1.7977e+308 B. 2.2251e-308 C. +251 D. e^2 10. 在MATLAB命令窗口输入X=logspace(1,2,99999);a=X(99999)/X(1) 大学物理实验试题（一） 一、单项选择题（每小题3分，共10小题） （1）.在光栅测量波长的实验中，所用的实验方法是[ ] （A）模拟法（B）干涉法（C）稳态法（D）补偿法 （2）.用箱式惠斯登电桥测电阻时，若被测电阻值约为4700欧姆，则倍率选[ ] （A）0.01 （B） 0.1 （C）10 （D） 1 （3）.用某尺子对一物体的长度进行15次重复测量，计算得A类不确定度为0.01mm，B类不确定度是0.6mm，如果用该尺子测量类似长度，应选择的合理测量次数为 （A）1次（B）6次（C）15次（D） 30次 （4）.用惠斯登电桥测电阻时，如果出现下列情况，试选择出仍能正常测 量的情况[ ] （A）有一个桥臂电阻恒为零（B）有一个桥臂电阻恒为无穷大 （C）检流计支路不通（断线）（D）电源与检流计位置互换 （5）.研究二极管伏安特性曲线时，正确的接线方法是[ ] （A）测量正向伏安特性曲线时用外接法；测量反向伏安特性曲线时用内接法（B）测量正向伏安特性曲线时用内接法；测量反向伏安特性曲线时用外接法（C）测量正向伏安特性曲线时用内接法；测量反向伏安特性曲线时用内接法（D）测量正向伏安特性曲线时用外接法；测量反向伏安特性曲线时用外接法（6）.在测量钢丝的杨氏模量实验中，预加1Kg砝码的目的是[ ] （A）消除摩擦力（B）使系统稳定 （C）拉直钢丝（D）增大钢丝伸长量 （7）.调节气垫导轨水平时发现在滑块运动方向上不水平，应该[ ] （A）只调节单脚螺钉（B）先调节单脚螺钉再调节双脚螺钉（C）只调节双脚螺钉（D）先调节双脚螺钉再调节单脚螺钉（8）.示波管的主要组成部分包括[ ] （A）磁聚集系统、偏转系统、显示屏（B）电子枪、偏转系统、显示屏（C）电聚集系统、偏转系统、显示屏（D）控制极、偏转系统、显示屏（9）.分光计设计了两个角游标是为了消除[ ] 数学软件与数学实验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) 《数学软件与数学实验》考试题型示例 一、单项选择题 1．利用赋值语句和表达式可完成某些复杂计算，例如在MA TLAB命令窗口中键入命令，Vname=sum(2.^[0:63])/(4.0e+10)，可计算出对应的数据，在这一语句中如果省略了变量名Vname及等号，MA TLAB将用缺省变量名（）显示计算结果 A）eps；B）ans；C）NaN；D）pi 2．下面有关MATLAB变量名和函数名的说法中，错误的说法是( ) (A) 变量名的第一个字符必须是一个英文字母 (B) 变量名可由英文字母、数字和下划线混合组成 (C) 变量名不得包含空格和标点，但可以有下连字符 (D) 变量名和函数名对于英文的大小使用没有区别 3．某城市电视塔地理位置：北纬30度35.343分，东经104度2.441分，在MA TLAB中用变量B=[30 35.343]表达纬度，用L=[104 2.441]表达经度。为了将经纬度数据转化为以度为单位的实数，下面正确的语句是（） A）P=B(1)+B(2)，Q=L(1)+ L(2)；B）P = 60*B(1) + B(2)，Q=60*L(1)+L(2)； C）P=B(1)+B(2)/60，Q=L(1)+L(2)；D）P = B(1) + B(2)/60，Q=L(1)+L(2)/60。 4．用MA TLAB随机产生一个10到20的正整数，应该使用下面的命令（）A）60+fix(40*rand)；B）10+20*rand；C）60+fix(100*rand)；D）10+round(10*rand) 5．用A、B、C表示三角形的三条边，MATLAB表示“任意两条边之和大于第三条边”的逻辑表达式正确的是（） A）A+B>=C | A+C>=B | A+C>=B；B）A+B<=C | A+C<=B | A+C<=B； C）A+B>C | A+C>B | B+C>A；D）A+B>C & A+C>B & B+C>A； 6．在MATLAB命令窗口中，键入命令syms x；y=int(3*x)。屏幕上将出现的结果是（）A）3/2*x^2；B）3x^2/2；C）1.5x^2；D）1.5*x^2； 7．在MATLAB命令窗口中，键入命令A=[1,2,3;4,5,6;7,8,0]；A(1,:)*A(:,3)。屏幕上将出现的结果是（） A）15；B）30；C）36；D）69； 8．正确表达命题A和B都大于C的逻辑表达式应该用下面哪一行（）A）A > C；B）B>C；C）A >C & B >C；D）A >C | B >C； 9．如果已输入方阵A的数据，在MA TLAB中用命令（）可计算出A的行列式的值A）det(A)；B）eig(A)；C）inv(A)；D）diag(A) 10．火炮发射炮弹的初始速度和发射角为已知，由此可估算出炮弹在空中的飞行时间Tfly，精品文档《数学实验》试题答案
数学软件MATLAB实验作业
大学数学数学实验(第二版)第7,8章部分习题答案
《数学软件》实验报告-符号计算基础与符号微积分
重庆大学数学实验 方程模型及其求解算法 参考答案
四川师范大学数学与软件科学学院程序设计实验报告实验九(推荐文档)
《数学建模与数学实验》上机报告
数学软件实验报告实验七
东华大学MATLAB数学实验第二版答案(胡良剑)
数学软件MATLAB实验报告 实验一
数学软件与实验 第一次上机作业
数学实验软件
《大学物理实验》模拟试卷与答案
南京邮电大学数学实验练习题参考答案
重庆大学数学实验一 matlab的基本应用 参考答案
数学软件与实验
大学物理实验及答案
数学实验与数学软件(Mathmaticandmatlab)
最新《数学软件与数学实验》考试题型示例及答案
相关主题
文本预览