数学软件与实验第一次上机作业
上机时间：2013-4-10 地点：E204
班级:071111 学号:07111014 姓名:曹红兴xdhjtang@https://www.doczj.com/doc/3215693674.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
4.生成以x 的元素为对角线的矩阵A，并输出A 的行数
>> A=diag(x)
a=length(A)
A =
1 0 0 0 0 0 0 0 0
0 3 0 0 0 0 0 0 0
0 0 10 0 0 0 0 0 0
0 0 0 15 0 0 0 0 0
0 0 0 0 20 0 0 0 0
0 0 0 0 0 25 0 0 0
0 0 0 0 0 0 30 0 0
0 0 0 0 0 0 0 35 0
0 0 0 0 0 0 0 0 40
a =
9
5.生成一个与A 同阶的正态分布的随机矩阵B, 要求B矩阵元素
为介于0~99之间整数
>> B=round(0+rand(9)*(99-0))
B =
75 10 79 65 25 83 47 38 46
74 82 19 16 50 25 35 56 1
39 69 48 12 69 81 82 8 33
65 31 44 49 88 24 58 5 16
17 94 64 95 95 92 54 53 79
70 3 70 34 54 35 91 77 31
3 43 75 58 1
4 19 28 92 52
27 38 27 22 15 25 75 13 16
5 7
6 6
7 74 25 61 75 56 60
6.输出A 与B 的乘积矩阵C
>> C=A*B
>> C=A*B
C =
5 99
6 10 89 30 50 82 72
54 84 156 222 21 135 225 99 210
660 410 330 730 240 100 620 290 770
495 690 255 840 75 1485 135 1110 435
1780 1520 420 360 880 660 160 20 1380
300 2025 2250 1475 25 725 1925 125 1375
2940 300 2010 900 2670 180 **** **** 1170
1855 630 1610 455 665 1050 1855 2100 210
2800 1440 3600 840 360 200 440 2080 3080
7.生成由A 与B 点乘得到的矩阵D
D=A.*B
D =
5 0 0 0 0 0 0 0 0
0 84 0 0 0 0 0 0 0
0 0 330 0 0 0 0 0 0
0 0 0 840 0 0 0 0 0
0 0 0 0 880 0 0 0 0
0 0 0 0 0 725 0 0 0
0 0 0 0 0 0
2700 0 0
0 0 0 0 0 0 0 2100 0
0 0 0 0 0 0 0 0 3080
8.生成一个由D 的第8、4、2、5 行和第7、1、6、9、2 列组成
的子矩阵 E
>> E=D([8 4 2 5],[7 1 6 9 2])
E =
0 0 0 0 0
0 0 0 0 0
0 0 0 0 246
0 0 0 0 0
9.求出矩阵E 的最大元素
>> m=max(max(E))
m =
246
10.计算的值
>> 1.369^2+sin(7*pi/10)*sqrt(26.48)/2.9
ans =
3.3097
11.产生一个5阶魔术方阵，并执行如下操作：
a)将矩阵的第2行3列元素赋值给变量c
b)将由矩阵第2，3，4行第3，5列构成的子矩阵赋值给变量d
>> Q=magic(5)
c=Q(2,3)
d=[Q(2:4,3),Q(2:4,5)]
Q =
17 24 1 8 15
23 5 7 14 16
4 6 13 20 22
10 12 19 21 3
11 18 25 2 9
c =
7
d =
7 16
13 22
19 3
12.给出区间[0，1]上的6个等分点数据。
>> L=0:1/5:1
L =
0 0.2000 0.4000 0.6000 0.8000 1.0000
13.建立如下矩阵
1)错误！未找到引用源。
2)错误！未找到引用源。
>> R=3.14*eye(6)
R =
3.1400 0 0 0 0 0
0 3.1400 0 0 0 0
0 0 3.1400 0 0 0
0 0 0 3.1400 0 0
0 0 0 0 3.1400 0
0 0 0 0 0 3.1400 >> A1=round(10+rand(6)*(10-10))
A1(1,1)=0,A1(2,2)=0,A1(3,3)=0,A1(4,4)=0,A1(5,5)=0,A1(6,6)=0
A=A1
A1 =
10 10 10 10 10 10
10 10 10 10 10 10
10 10 10 10 10 10
10 10 10 10 10 10
10 10 10 10 10 10
10 10 10 10 10 10
A1 =
0 10 10 10 10 10
10 10 10 10 10 10
10 10 10 10 10 10
10 10 10 10 10 10
10 10 10 10 10 10
10 10 10 10 10 10
A1 =
0 10 10 10 10 10
10 0 10 10 10 10
10 10 10 10 10 10
10 10 10 10 10 10
10 10 10 10 10 10
10 10 10 10 10 10 A1 =
0 10 10 10 10 10
10 0 10 10 10 10
10 10 0 10 10 10
10 10 10 10 10 10
10 10 10 10 10 10
10 10 10 10 10 10 A1 =
0 10 10 10 10 10
10 0 10 10 10 10
10 10 0 10 10 10
10 10 10 0 10 10
10 10 10 10 10 10
10 10 10 10 10 10 A1 =
0 10 10 10 10 10
10 0 10 10 10 10
10 10 0 10 10 10
10 10 10 0 10 10
10 10 10 10 0 10
10 10 10 10 10 10 A1 =
0 10 10 10 10 10
10 0 10 10 10 10
10 10 0 10 10 10
10 10 10 0 10 10
10 10 10 10 0 10
10 10 10 10 10 0
A =
0 10 10 10 10 10
10 0 10 10 10 10
10 10 0 10 10 10
10 10 10 0 10 10
10 10 10 10 0 10
10 10 10 10 10 0
14.已知函数错误！未找到引用源。,计算f(-1), f(0.5), f(1.5)。
function y=f(x)
n=length(x);
for k=1:n
if x(k)>=1
y(k)=x(k)^2;
elseif x(k)>=0
y(k)=1;
else
y(k)=x(k)+1;
end
end
>> x=[-1,0.5,1.5] y=f(x)
x =
-1.0000 0.5000 1.5000
y =
0 1.0000 2.2500
15.用for-end循环语句求错误！未找到引用源。和100!。
16.用while-end循环语句求不超过1000的偶数之和与奇数之和。
偶数之和：
>> s=0;
>> i=0;
>> while i<=1000
s=s+i;
i=i+2;
end
>> s
s =
250500
>>
奇数之和：
>> s=0;
>> i=1;
>> while i<=1000
s=s+i;
i=i+2;
end
>> s
s =
250000
>>
17.建立一个命令M-文件：求所有的“水仙花数”，所谓“水仙花数”
是指一个三位数，其各位数字的立方和等于该数本身。例如，153是一个水仙花数，因为错误！未找到引用源。
数学软件与数学实验作业 一.《数学软件》练习题（任选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)
“”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> ('(x)-3*x^2',0) = -2*(-1/6*3^(1/2)) -2*(-11/6*3^(1/2)) -2*(1/6*3^(1/2)) 3、求解下列各题： 1）30 sin lim x x x x ->- >> x;
>> (((x))^3) = 1/6 2） (10)cos ,x y e x y =求 >> x; >> ((x)*(x),10) = (-32)*(x)*(x) 3）2 1/2 0(17x e dx ?精确到位有效数字） >> x; >> ((((x^2),0,1/2)),17) =
0.54498710418362222 4）4 2 254x dx x +? >> x; >> (x^4/(25^2)) = 125*(5) - 25*x + x^3/3 5）求由参数方程arctan x y t ??=? =??dy dx 与二阶导 数22 d y dx 。 >> t; >> ((1^2))(t); >> ()() = 1
6）设函数(x)由方程e所确定，求y′(x)。>> x y; *(y)(1); >> ()() = (x + (y)) 7） sin2 x e xdx +∞- ? >> x; >> ()*(2*x); >> (y,0) = 2/5
8） 08x =展开（最高次幂为） >> x (1); taylor(f,0,9) = - (429*x^8)/32768 + (33*x^7)/2048 - (21*x^6)/1024 + (7*x^5)/256 - (5*x^4)/128 + x^3/16 - x^2/8 + 2 + 1 9） 1sin (3)(2)x y e y =求 >> x y; >> ((1)); >> ((y,3),2) =
2017春季数学实验报告 班级：计算机系61 姓名：赵森学号：2160500026（校内赛编号506）班级：计算机系61 姓名：冯丹妮学号：2160500002（校内赛编号327）班级：计算机系63 姓名：郝泽霖学号：2160500054
第一次上机作业 实验8: 练习1： 4．某棉纺厂的原棉需从仓库运送到各车间。各车间原棉需求量、单位产品从各仓库运往各车间的运输费以及各仓库的库存容量如表8.5所列，问如何安排运输任务使得总运费最小？ 设仓库1运往车间1,2,3,的原棉量为x1,x2,x3, 仓库2运往车间1,2,3,的原棉量为x4,x5,x6, 仓库3运往车间1,2,3,的原棉量为x7,x8,x9。 2x1+x2+3x3<=50 2x4+2x5+4x6<=30 3x7+4x8+2x9<=10 X1+x4+x7=40 X2+x5+x8=15 X3+x6+x9=35 程序： c=[2,1,3,2,2,4,3,4,2]; a(1,:)=[1,1,1,0,0,0,0,0,0]; a(2,:)=[0,0,0,1,1,1,0,0,0]; a(3,:)=[0,0,0,0,0,0,1,1,1]; aeq(1,:)=[1,0,0,1,0,0,1,0,0]; aeq(2,:)=[0,1,0,0,1,0,0,1,0]; aeq(3,:)=[0,0,1,0,0,1,0,0,1]; b=[50;30;10]; beq=[40;15;35]; vub=[]; vlb=zeros(9,1); [x,fval]=linprog(c,a,b,aeq,beq,vlb,vub) 结果： x = 10.0000 15.0000 25.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);
学生实验报告 实验课名称： C++程序设计 实验项目名称：综合大作业——学生成绩管理系统专业名称：电子信息工程 班级： 学号： 学生： 同组成员： 教师：
2011 年 6 月 23 日 题目：学生成绩管理系统 一、实验目的： （1）对C++语法、基础知识进行综合的复习。 （2）对C++语法、基础知识和编程技巧进行综合运用，编写具有一定综合应用价值的稍大一些的程序。培养学生分析和解决实际问题的能力，增强学生的自信心，提高学生学习专业课程的兴趣。 （3）熟悉掌握C++的语法和面向对象程序设计方法。 （4）培养学生的逻辑思维能力，编程能力和程序调试能力以及工程项目分析和管理能力。 二、设计任务与要求： （1）只能使用/C++语言，源程序要有适当的注释，使程序容易阅读。 （2）至少采用文本菜单界面（如果能采用图形菜单界面更好）。 （3）要求划分功能模块，各个功能分别使用函数来完成。 三、系统需求分析： 1.需求分析: 为了解决学生成绩管理过程中的一些简单问题，方便对学生成绩的管理 （录入，输出，查找，增加，删除，修改。） 系统功能分析: (1):学生成绩的基本信息：学号、、性别、C++成绩、数学成绩、英语成绩、 总分。 (2):具有录入信息、输出信息、查找信息、增加信息、删除信息、修改信息、 排序等功能。 2.系统功能模块（要求介绍各功能） （1）录入信息(Input)：录入学生的信息。 （2）输出信息(Print)：输出新录入的学生信息。 （3）查找信息(Find)：查找已录入的学生信息。 （4）增加信息(Add)：增加学生信息。 （5）删除信息(Remove)：在查找到所要删除的学生成绩信息后进行删除并输出删除后其余信息。 （6）修改信息(Modify)：在查到所要修改的学生信息后重新输入新的学生信息从而进行修改，然后输出修改后的所有信息。 （7）排序(Sort)：按照学生学号进行排序。 3.模块功能框架图
练习2﹒1 画出下列常见曲线的图形（其中a=1，b=2,c=3）。 1. 立方抛物线y = 解： x=-4:0.1:4; y=x.^(1/3); plot(x,y) -4 -3-2-101234 0.20.40.60.811.21.4 1.6 2.高斯曲线2 x y e -= 解： fplot('exp(-x^2)',[-4,4])
-4 -3 -2 -1 1 2 3 4 00.10.20.30.40.50.60.70.80.9 1 3、笛卡儿曲线23 3 2 2 33,(3)11at at x y x y axy t t = = +=++ 解：ezplot('x^3+y^3-3*x*y',[-4,4])
-4 -3-2-1 01234 -4-3-2-10123 4x y x 3+y 3-3 x y = 0 或：t=-4:0.1:4; x=3*t./(1+t.^2); y=3*t.^2./(1+t.^2); plot(x,y)
-1.5 -1-0.500.51 1.5 00.5 1 1.5 2 2.5 3 4、蔓叶线233 2 2 2 ,()11at at x x y y t t a x = = = ++- 解：t=-4:0.1:4; x=t.^2./(1+t.^2); y=t.^3,/(1+t.^2); y=t.^3./(1+t.^2); plot(x,y)
00.10.20.30.40.50.60.70.80.91 -4 -3-2-10123 4 或： ezplot('y .^2-x.^3/(1-x)',[-4,4])
《数学建模与数学实验》上机报告（第 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]
“MATLAB”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> solve('exp(x)-3*x^2',0) ans = -2*lambertw(-1/6*3^(1/2)) -2*lambertw(-1,-1/6*3^(1/2)) -2*lambertw(1/6*3^(1/2)) 2、求下列方程的根。 1) 5510x x ++= a=solve('x^5+5*x+1',0);a=vpa(a,6)
1.10447+1.05983*i -1.00450+1.06095*i -.199936 -1.00450-1.06095*i 1.10447-1.05983*i 2） 1 sin0 2 x x-=至少三个根 >> fzero('x*sin(x)-1/2', 3) ans = 2.9726 >> fzero('x*sin(x)-1/2',-3) ans = -2.9726 >> fzero('x*sin(x)-1/2',0) ans = -0.7408
3）2sin cos 0x x x -= 所有根 >> fzero('sin(x)*cos(x)-x^2',0) ans = >> fzero('sin(x)*cos(x)-x^2',0.6) ans = 0.7022 3、求解下列各题： 1）30sin lim x x x x ->- >> sym x; >> limit((x-sin(x))/x^3) ans = 1/6 2） (10)cos ,x y e x y =求 >> sym x; >> diff(exp(x)*cos(x),10) 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函数的功能
重庆科技学院 数学实验与数学软件课程设计 课程名称：菜单与对话框设计 开课学期：_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函数将出现我们所需要的图像，如图：
《数学实验与MATLAB》 ——综合实验报告 实验名称：不同温度下PDLC薄膜的通透性 与驱动电压的具体关系式的研究学院：计算机与通信工程学院 专业班级： 姓名： 学号： 同组同学： 2014年 6月10日
一、问题引入 聚合物分散液晶(PDLC)是将低分子液晶与预聚物Kuer UV65胶相混合，在一定条件下经聚合反应，形成微米级的液晶微滴均匀地分散在高分子网络中，再利用液晶分子的介电各向异性获得具有电光响应特性的材料，它主要工作在散射态和透明态之间并具有一定的灰度。聚合物分散液晶膜是将液晶和聚合物结合得到的一种综合性能优异的膜材料。该膜材料能够通过驱动电压来控制其通透性，可以用来制作PDLC型液晶显示器等，具有较大的应用范围。已知PDLC薄膜在相同光强度及驱动电压下，不用的温度对应于不同的通透性，不同温度下的阀值电压也不相同。为了尽量得到不同通透性的PDLC薄膜，有必要进行温度对PDLC薄膜的特性的影响的研究。现有不同温度下PDLC 薄膜透过率与驱动电压的一系列数据，试得出不同温度下PDLC薄膜通透性与驱动电压的具体关系式，使得可以迅速得出在不同温度下一定通透性对应的驱动电压。 二、问题分析 想要得到不同温度下PDLC薄膜通透性与驱动电压的具体关系式可以运用MATLAB多项式农合找出最佳函数式，而运用MATLAB多项式插值可以得出在不同温度下一定通透性所对应的驱动电压。 三、实验数据 选择10、20、30摄氏度三个不同温度，其他条件一致。
（1）、10摄氏度 实验程序： x=2:2:40; y=[5.2,5.4,5.8,6.4,7.2,8.2,9.4,10.8,12.2,14.0,16.6,22.0, 30.4,39.8,51.3,55.0,57.5,58.8,59.6,60.2]; p3=polyfit(x,y,3); p5=polyfit(x,y,5); p7=polyfit(x,y,7); disp('三次拟合函数'),f3=poly2str(p3,'x') disp('五次拟合函数'),f5=poly2str(p5,'x') disp('七次拟合函数'),f7=poly2str(p7,'x') x1=0:1:40; y3=polyval(p3,x1); y5=polyval(p5,x1); y7=polyval(p7,x1); plot(x,y,'rp',x1,y3,'--',x1,y5,'k-.',x1,y7); legend('拟合点','三次拟合','五次拟合','七次拟合') 实验结果：
练习6.7 1.有两个煤厂A,B，每月进煤不少于60t,100t，它们担负供应三个居 民区的用煤任务，这三个居民区每月用煤量分别为45t,75t和45t.A 厂离这三个居民区的距离分别为10km,5km,6km,B厂离这三个居民区的距离分别为4km,8km,15km.问这两个煤厂如何分配供煤量能使总运输量（t.km）最小。 解：设甲对三个居民区的供煤量分别为：x1,x2,x3,乙对三个居民区的供煤量分别为x4,x5,x6.由已知有： y=10x1+5x2+6x3+4x4+8x5+15x6 -x1-x2-x3<=-60, -x4-x5-x6<=-100, x1+x4=45,x2+x5=75,x3+x6=40, X1>=0,x2>=0,x3>=0,x4>=0,x5>=0,x6>=0. 输入命令： > c=[10 5 6 4 8 15];A=[-1 -1 -1 0 0 0;0 0 0 -1 -1 -1;0 0 0 0 0 0;0 0 0 0 0 0;0 0 0 0 0 0;0 0 0 0 0 0]; >> b=[-60;-100;0;0;0;0];Aeq=[1 0 0 1 0 0;0 1 0 0 1 0;0 0 1 0 0 1;0 0 0 0 0 0;0 0 0 0 0 0;0 0 0 0 0 0]; >> beq=[45 75 40 0 0 0]; >> lb=ones(6,1); >> [x,fval]=linprog(c,A,b,Aeq,beq,lb) Optimization terminated.
结果为： x = 1.0000 20.0000 39.0000 44.0000 55.0000 1.0000 fval =975.0000 这说明甲乙两个煤厂分别对三个居民区输送1t 20t 39t,44t 55t 1t的煤才能使总运输量最小，且总运输量为975t.km 2．某银行经理计划用一笔资金进行有价证券的投资，可供购进的证券及其信用等级、到期年限、税前收益如下表所示。按照规定，市政证券的收益可以免税，其他证券的收益需按40%的税率纳税。此外还有以下限制： （1）政府及待办机构的证券总共至少购进400万元； （2）所构证券的平均信用等级不超过1.4（信用等级数字越小，信用程度越高）； （3）所构证券的平均到期年限不超过5年。
《数学实验》报告 实验名称李萨如图模拟（Matlab大作业） 2011年11月8日
一、【实验目的】 运用数学知识与MATLAB相结合，运用数学方法，建立数学模型，用MATLAB软件辅助求解模型，解决实际问题。 二、【实验任务】 一个质点沿 X轴和 Y轴的分运动都是简谐运动,分运动的表达式分别为: x=Acos ( w1t+beta ) , y=Acos(w2t+beta ) 。如果二者的频率有简单的整数比, 则相互垂直的简谐运动合成的运动将具有封闭的稳定的运动轨迹, 这种图称为李萨如图。 1，用matlab分别画出同一方向的传播波频率之比为2,3,4/5,1/2,1/3,5/4的图像（未合成）2，用matlab画出同一方向的传播波频率之比为2,3,4/5,1/2,1/3,5/4的合成图像 3，用matlab画出x轴方向和y轴方向传播波频率之比为2,3,4/5,1/2,1/3,5/4的合成图像。（李萨如图） 三、【实验分析及求解】 1，设两个波的振幅为1，他们的beta为pi/5,我们可以根据波的传播公式，y =Acos ( w1t+beta ) 分别画出两个波的传播图像。 2，设两个波的振幅为1，他们的beta为pi/5,我们可以根据波的传播公式，y =Acos ( w1t+beta )， 用matlab画出同一方向的传播波频率之比为2,3,4/5,1/2,1/3,5/4的合成图像。
3，设两个波的振幅为1，他们的beta为pi/5,我们可以根据波的传播公式，画出x轴方向和y 轴方向传播波频率之比为2,3,4/5,1/2,1/3,5/4的合成图像。（李萨如图）。
4. 问题： 某公司将3种不同含硫量的液体原料（分别记为甲、乙、丙）混合生产两种产品（分别记为A,B ）。按照生产工艺的要求，原料甲、乙必须首先导入混合池中混合，混合后的液体再分别与原料丙混合生产A,B 。一直原料甲、乙、丙的含硫量分别是3%，1%，2%，进货价格分别为6千元/t ，16千元/t ，10千元/t ；产品A,B 的含硫量分别不能超过2.5%，1.5%，售价分别为9千元/t ，15千元/t 。根据市场信息，原料甲、乙、丙的供应量都不能超过500t ；产品A,B 的最大市场需求量分别为100t ，200t 。 (1) 应如何安排生产？ (2) 如果产品A 的最大市场需求量增长为600t ，应如何安排生产？ (3) 如果乙的进货价格下降为13千元/t ，应如何安排生产？分别对(1)、(2)两种情况进 行讨论。 模型： （只考虑问题1，问题2,3只需改变一些约束条件） 设生产时使用原料甲、乙分别为12,x x t ，分别取混合后的液体34,x x t 再加入原料丙 56,x x t 生产产品A,B 。 有质量守恒，可得 1234x x x x +=+ 甲乙混合后的液体的含硫量可表示为 12 12 3%x x x x ++，根据含硫量的要求，可得 12 353512 124646 12 3%*2%* 2.5%*()3%*2%* 1.5%*() x x x x x x x x x x x x x x x x +?+≤+?+?? +?+≤+?+? 根据市场的限制，易得 12563546500 500500100200 x x x x x x x x ≤?? ≤?? +≤??+≤??+≤? 当然还有非负约束 123456,,,,,0x x x x x x ≥ 公司的净利润为（单位：千元）： 35461256123456 9()15()61610()6169155z x x x x x x x x x x x x x x =+++---+=--++-+
数学软件与数学实验_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)
“M A T L A B ”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> solve('exp(x)-3*x^2',0) ans = -2*lambertw(-1/6*3^(1/2)) -2*lambertw(-1,-1/6*3^(1/2)) -2*lambertw(1/6*3^(1/2)) 2、求下列方程的根。 1) 5510x x ++= a=solve('x^5+5*x+1',0);a=vpa(a,6) a = 1.10447+1.05983*i -1.00450+1.06095*i -. -1.00450-1.06095*i
1.10447-1.05983*i 2） 1 sin0 2 x x-=至少三个根 >> fzero('x*sin(x)-1/2', 3) ans = 2.9726 >> fzero('x*sin(x)-1/2',-3) ans = -2.9726 >> fzero('x*sin(x)-1/2',0) ans = -0.7408 3）2 sin cos0 x x x -=所有根 >> fzero('sin(x)*cos(x)-x^2',0) ans = >> fzero('sin(x)*cos(x)-x^2',0.6)
0.7022 3、求解下列各题： 1）3 0sin lim x x x x ->- >> sym x; >> limit((x-sin(x))/x^3) ans = 1/6 2） (10)cos ,x y e x y =求 >> sym x; >> diff(exp(x)*cos(x),10) ans = (-32)*exp(x)*sin(x) 3）2 1/2 0(17x e dx ?精确到位有效数字） >> sym x; >> vpa((int(exp(x^2),x,0,1/2)),17)
数学软件与实验第一次上机作业 上机时间：2013-4-10 地点：E204 班级:071111 学号:07111014 姓名:曹红兴xdhjtang@https://www.doczj.com/doc/3215693674.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
数学实验大作业——抽象群与应用“RSA加密系统” 合作人：郭元镇尹庆宇杨瑞飞 综述 1）RSA 加密算法的历史 RSA公钥加密算法是1977年由Ron Rivest、Adi Shamirh和LenAdleman在（美国麻省理工学院）开发的。RSA取名来自开发他们三者的名字。RSA是目前最有影响力的公钥加密算法，它能够抵抗到目前为止已知的所有密码攻击，已被ISO推荐为公钥数据加密标准。RSA算法基于一个十分简单的数论事实：将两个大素数相乘十分容易，但那时想要对其乘积进行因式分解却极其困难，因此可以将乘积公开作为加密密钥。这种算法1978年就出现了，它是第一个既能用于数据加密也能用于数字签名的算法。它易于理解和操作，也很流行。早在1973年，英国国家通信总局的数学家Clifford Cocks就发现了类似的算法。但是他的发现被列为绝密，直到1998年才公诸于世。 RSA是被研究得最广泛的公钥算法，从提出到现在已近二十年，经历了各种攻击的考验，逐渐为人们接受，普遍认为是目前最优秀的公钥方案之一。RSA的安全性依赖于大数的因子分解，但并没有从理论上证明破译RSA的难度与大数分解难度等价。即RSA的重大缺陷是无法从理论上把握它的保密性能如何，而且密码学界多数人士倾向于因子分解不是NPC问题。 2）RSA 加密算法的原理 RSA算法是一种非对称密码算法，所谓非对称，就是指该算法需要一对密钥，使用其中一个加密，则需要用另一个才能解密。 RSA的算法涉及三个参数，n、e1、e2。 其中，n是两个大质数p、q的积，n的二进制表示时所占用的位数，就是所谓的密钥长度。 e1和e2是一对相关的值，e1可以任意取，但要求e1与(p-1)*(q-1)互质；再选择e2，要求(e2*e1)mod((p-1)*(q-1))=1。 (n及e1),(n及e2)就是密钥对。
数学软件与数学实验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)