matlab实验题目1

合集下载

Matlab数学实验1简单函数曲线解读

Matlab数学实验1简单函数曲线解读

>> plot(x,cos(x),’bd-.’); 可以只指定其中某几个 >> plot(x,cos(x),’k*-’); 排列顺序任意
蓝色、点划线 离散点为菱形
黑色、实线 离散点用星号
线型
- 实线 : 虚线 -. 点划线 -- 间断线
点标记
.点 o 小圆圈 x 叉子符 + 加号 * 星号 s 方格 d 菱形 ^ 朝上三角 v 朝下三角 > 朝右三角 < 朝左三角 p 五角星 h 六角星
所得图形如下所示:
ezpolar是简易极坐标作图命令。也可以把上面的输入改为: ezpolar(‘3*cos(3*t)’),运行此程序也可得到上面的图形。
课堂上机练习
1. 作出函数y tan x和y cot x的图形观察其周期性和变化趋势。
2. 在区间[1,1]内画出函数y sin 1 的图形。 x
颜色
y 黄色 m 棕色 c 青色 r 红色 g 绿色 b 蓝色 w 白色 k 黑色
[0,2 ]
作图命令: x=linspace(0,2*pi,30); y=sin(x); z=cos(x); plot(x,y,'r',x,z,'go');
注:linspace是Matlab中的一个指令,linspace(x1,x2,N)用 于产生x1,x2之间的N点行矢量。其中x1、x2、N分别为起始值、 终止值、元素个数。若缺省N,默认点数为100 。
ezplot(‘f(x,y)’)
f (x, y) 0 在区间 2 x 2 ,2 y 2 上的图形
ezplot(‘f’,’g’,[a,b])
x f (t), y g(t) 在区间 a t b 上的图形

matlab课堂实验小测1带答案

matlab课堂实验小测1带答案

课堂实验小测1◆1、MATLAB 中包含两个函数upper 和lower ,分别把字符串转化为大字和小写。

创建一个新的函数caps ,让字符串每个单词的第一个字母大写,其余为小写。

str=input('Please enter the words: ','s');str(1)=upper(str(1));str(2:end)=lower(str(2:end));disp(str);◆2、工程师们经常用分贝或dB 来描述两功率之比。

1dB 的定义如下:210110log P dB P = P2 是已测量的功率, P1 代表参考功率。

假设参考功率P1 为1mw ,编写一个程序,接受一个输入功率P2 并把转化成为以1mw 为参考功率的dB.(它在工程上有一个特殊单位dBm)。

P2=input('Please enter the P2(mW): ');P1=1;PdB=10*log10(P2/P1)disp('The power is');disp(PdB);disp('mW');◆3、弹簧中的能量.压缩弹簧弹力的大小可由下面的公式计算出来F = kxF 代表弹力,单位为N 。

k 代表劲度系数单位为N/m.存储在压力弹簧中的势能为212E kx = E 代表势能,单位为焦。

下面是4个可用压缩弹簧的信息.编写一个程序,计算出每一个弹簧的压缩量和弹力势能.哪一个弹簧的弹力最大F=[20,24,22,20];k=[500,600,700,800];x=F./k;disp('压缩量:');disp(x);E=1/2*k.*x.^2;disp('势能:')E[SE,ID]=sort(E);disp('最大弹力弹簧号:')ID(end)◆4、把角度的弧度格式转换为相应的度/分/秒。

编写一程序,从磁盘文件中读取弧度格式的角度,并把它转化为相应的度/分/秒格式。

(完整版)MATLAB)课后实验[1]

(完整版)MATLAB)课后实验[1]

实验一 MATLAB 运算基础1. 先求以下表达式的值,尔后显示 MATLAB 工作空间的使用情况并保存全部变量。

(1)2sin 85 z1 21 e(2) 12z ln( x 1 x ) ,其中22 x2 1 2i5(3)ae e az sin( a 0.3) ln , a 3.0, 2.9, L , 2.9, 32 22t 0 t 1(4) 2z t 1 1 t 242t 2t 1 2 t 3,其中解:M 文件:z1=2*sin(85*pi/180)/(1+exp(2))x=[2 1+2*i;-.45 5];z2=1/2*log(x+sqrt(1+x^2))a=-3.0:0.1:3.0;z3=(exp(0.3.*a)-exp(-0.3.*a))./2.*sin(a+0.3)+log((0.3+a)./2)t=0:0.5:2.5;z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3) .*(t.^2-2*t+1)4. 完成以下操作:(1) 求[100,999] 之间能被 21 整除的数的个数。

(2) 建立一个字符串向量,删除其中的大写字母。

解:(1) 结果:m=100:999;n=find(mod(m,21)==0);length(n)ans =43(2). 建立一个字符串向量比方:ch='ABC123d4e56Fg9'; 那么要求结果是:ch='ABC123d4e56Fg9';k=find(ch>='A'&ch<='Z');ch(k)=[]ch =123d4e56g9实验二 MATLAB矩阵解析与办理1. 设有分块矩阵 A E R3 3 3 2O S2 3 2 2,其中 E、R、O、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵,试经过数值计算考据 2A E R RS2O S。

MATLAB实验一:参考--答案

MATLAB实验一:参考--答案

save wenjian.mat A a x load wenjian.mat A a x 小结及思考 在 matlab 中,A*B 与 A .*B 以及 B./A 与 B.\A 之间有什么区别?
A*B 是 A 与 B 进行通常意义上的矩阵乘法得到的矩阵,而 A.*B 是矩阵 A 的每个 元素与 B 的对应位置的元素相乘得到的矩阵。 B./A 是矩阵 B 的每一个元素除以 A 对应位置的元素得到的矩阵, B.\A 是矩阵 A 的每一个元素除以 B 对应位置的元素得到的矩阵。
(1)A=2015*eye(10) (2)A=8*(ones(8)-eye(8))
0 8 (2) 8
8 0 8

8 8 0 88
9. 设有分块矩阵 A
E33 O23
R32 ,其中 E,R,O,S 分别为单位阵、随机阵、零阵和对角阵, S 22
7. 写出下列函数的作用: clc Clear ones det triu fliplr zeros cat
清除命令窗口的记录 清除已声明的变量和函数的记录 生成全 1 矩阵 求一个方阵的行列式 取一个矩阵的上三角 对矩阵进行左右旋转 产生一个零矩阵 拼接两个数组或者矩阵
8. 建立如下矩阵
0 0 2015 2015 0 0 (1) 0 0 2015 1010
3.
计算 1.369 sin
2
7 26.48 2.9 的值(输入程序) 10
1.369^2+sin(7/10*pi)*26.48^(1/2)/2.9
4.
1 3 3 1 1 2 输入矩阵 A 4 4 6 , B 2 2 2 ,指出下列命令的含义 6 8 9 3 3 3

数学实验matlab练习1

数学实验matlab练习1

《数学实验》第一讲实验习题1.执行下列指令,观察其运算结果,理解其意义:(1)[1 2;3 4]+10-2ians =9 1011 12(2)[1 2;3 4].*[0.1 0.2;0.3 0.4]ans =0.1000 0.40000.9000 1.6000(3)[1 2;3 4].\[20 10;9 2]ans =20.0000 5.00003.0000 0.5000(4)[1 2;3 4].^2ans =1 49 16(5)exp([1 2;3 4])ans =2.7183 7.389120.0855 54.5982(6)log([1 10 100])ans =0 2.3026 4.6052(7)prod([1 2;3 4])ans =3 8(8)[a,b]=min([10 20;30 40])a =10 20b =1 1(9)abs([1 2;3 4]-pi)ans =2.1416 1.14160.1416 0.8584(10) [1 2;3 4]>=[4 3;2 1]ans =0 01 1(11) find([10 20;30 40]>=[40 30;20 10])ans =24(12) [a,b]=find([10 20;30 40]>=[40 30;20 10])a =22b =12(13) all([1 2;3 4]>1)ans =0 1(14) linspace(3,4,5)ans =3.0000 3.2500 3.5000 3.7500(15) A=[1 2;3 4];A(:,2)ans =242. 执行下列指令,观察其运算结果、变量类型和字节数,理解其意义:(1) clear;a=1,b=num2str(a),c=a>0,a= =b,a= =c,b= =cans =(2) clear;fun=’abs(x)’,x=-2,eval(fun),double(fun)ans =2ans =97 98 115 40 120 413. 本金K 以每年n 次,每次%p 的增值率(np 为每年增值额的百分比)增加,当增加到rK 时所花费的时间为(单位:年) 1001ln ln(.)T r n p += ln()ln(10.01)r T n p =+ 用MATLAB 表达式写出该公式并用下列数据计算:20512,.,r p n ===。

matlab实验1

matlab实验1

一、书上命令练习:用diary 记录下所以练习内容A=[1,2,3;4,5,6]A =1 2 34 5 6A=[1 2 3;4 5 6]A =1 2 34 5 6A=[1 2 3]A =1 2 3A=[1 2 3 ]A =1 2 3A=[1 2 34 5 6]A =1 2 34 5 6 a=A(2,1)a =4A(2,3)ans =6b=A(3)b =2A(4)ans =5A(2,1)=7A =1 2 37 5 6A(3,4)=1A =1 2 3 07 5 6 00 0 0 1W=zeros(2,3)W =0 0 00 0 0u=onea(3)??? Undefined function or variable 'onea'. u=ones(3)u =1 1 11 1 11 1 1v=eye(3,4)v =1 0 0 00 1 0 00 0 1 0x=rand(1,3)x =0.9501 0.2311 0.6068 A(3,:)ans =0 0 0 1A(:,2)ans =25b=A(1:2,:)b =1 2 3 07 5 6 0c=B(:,2:4)??? Undefined function or variable 'B'. C=b(:,2:4)C =2 3 05 6 0D=A(2:end,[2,4])D =5 00 1d(:,1)=[]??? Index of element to remove exceeds matrix dimensions.D(:,1)=[]D =1E=[C,ones(2,3)]E =2 3 0 1 1 15 6 0 1 1 1 F=[A(1:2,:);eye(1,4)]F =1 2 3 07 5 6 01 0 0 0G=[C,zeros(2);9,F(2,:)]G =2 3 0 0 05 6 0 0 09 7 5 6 0H=C(:)H =2536E=E+3E =5 6 3 4 4 48 9 3 4 4 4 A=[1,0,2;3,4,0]A =1 0 23 4 0 B=E(:,1:3)B =5 6 38 9 3A.*Bans =5 0 624 36 0B.^Aans =5 1 9512 6561 1A.\BWarning: Divide by zero.(Type "warning off MATLAB:divideByZero" to suppress this warning.)ans =5.0000 Inf 1.50002.6667 2.2500 InfB.\Aans =0.2000 0 0.66670.3750 0.4444 02.^Aans =2 1 48 16 1A.^2ans =1 0 49 16 0a=1:5a =1 2 3 4 5 b=1:2:7b =1 3 5 7c=6:-3:-6c =6 3 0 -3 -6b=[0:2:8,ones(1,3)]b =0 2 4 6 8 11 1linspace(0,1,9)ans =Columns 1 through 60 0.1250 0.2500 0.3750 0.5000 0.6250Columns 7 through 90.7500 0.8750 1.0000x=0:pi/4:pix =0 0.7854 1.5708 2.3562 3.1416x=linspace(0,pi,5)x =0 0.7854 1.5708 2.3562 3.1416a=[1 2 3 4 5];b=[1 3 5 7 9]b =1 3 5 7 9a=[1 2 3 4 5];b=[1 3 5 7 9];...c=a.*b,d=a*b',e=a'*bc =1 6 15 28 45d =95e =1 3 5 7 92 6 10 14 183 9 15 21 274 12 20 28 365 15 25 35 45a=[0 1 0],b=[1 0 0],c=a./ba =0 1 0b =1 0 0Warning: Divide by zero.(Type "warning off MATLAB:divideByZero" to suppress this warning.)c =0 Inf NaNs1='Hello';s2='every';s3='body';s=[s1,',',s2s 3]'ss=ss=s(1:5)???s1='Hello';s2='every';s3='body';s=[s1,',',s2s 3]'ss=ss=s(1:5)|Error: Missing operator, comma, or semicolon.s1='Hello';s2='every';s3='body';s=[s1,',',s2s 3]'ss=ss=s(1:5)???s1='Hello';s2='every';s3='body';s=[s1,',',s2s 3]'ss=ss=s(1:5)|Error: Missing operator, comma, orsemicolon.s1=;s2=;s3=;s=[s1,,s2,s3]'ss=ss=s(1:5) ??? s1=;s2=;s3=;s=[s1,,s2,s3]'ss=ss=s(1:5)|Error: Expected a variable, function, or constant, found ";".s1='Hello';s2='every';s3='body';s=[s1,',',s2,' ',s3],ss=s(1:5)s =Hello,every bodyss =Hellox=(0:0.2:1)*pi;y=sin(x)y =0 0.5878 0.9511 0.9511 0.5878 0.0000a=[-3.5 4.6];...b=round(a),c=floor(a);d=ceil(a),e=fix(a),f=r ats(a)b =-4 5d =-3 5e =-3 4f =-7/2 23/5x=(0:0.2:1)*pi;y=feval('sin',x)y =0 0.5878 0.9511 0.9511 0.5878 0.0000x=(0:0.2:1)*pi;y=feval('sin',x)y =0 0.5878 0.9511 0.9511 0.5878 0.0000x=(0:0.2:1)*pi;y=feval(@sin,x)y =0 0.5878 0.9511 0.9511 0.5878 0.0000x=(0:0.2:1)*pi;y=feval(@sin,x)y =0 0.5878 0.9511 0.9511 0.5878 0.0000x=(0:0.2:1)*pi;y=feval(@sin,x)y =0 0.5878 0.9511 0.9511 0.5878 0.0000x=(0:0.2:1)*pi;y=feval(inline("sin(x)+2'),x) ???x=(0:0.2:1)*pi;y=feval(inline("sin(x)+2'),x)| Error: Missing variable or function.x=(0:0.2:1)*pi;y=feval(inline('sin(x)+2'),x)y =2.0000 2.5878 2.9511 2.9511 2.5878 2.0000a=[4 3.1 -1.2 0 6];...b=min(a),c=sum(a),d=median(a),e=sort(a)b =-1.2000c =11.9000d =3.1000e =-1.2000 0 3.1000 4.0000 6.0000f=[1:3;4:6;7:9];f1=prod(f),f2=prod(f1)f1 =28 80 162f2 =362880a=[1 2 3;4 5 6;7 8 9;10 11 12];b=reshape(a,2,6)b =1 72 83 94 105 116 12sa=size(a),sb=size(b)sa =4 3sb =2 6c=reshape(1:9,3,3)c =1 4 72 5 83 6 9help sqtrtsqtrt.m not found.help sqrtSQRT Square root.SQRT(X) is the square root of the elements of X. Complexresults are produced if X is not positive.See also SQRTM.Overloaded methodshelp sym/sqrt.mwhoYour variables are:A D G a c f s s3 ss xB E H ans d f1 s1 sa u yC F W b e f2 s2 sb vdisp??? Error using ==> dispNot enough input arguments.disp(x)0 0.6283 1.25661.88502.51333.1416cleardiary off> A=[1,2,3;4,5,6]??? diary off> A=[1,2,3;4,5,6]|Error: Missing operator, comma, or semicolon.diary off二、MatLab 练习11, (5*2+1.3-0.8)*10^2/25 答案:5*2+1.3-0.8)*10^2/25ans =422,y =23log )10sin(4)43.0(2+⋅+-e π答案:y=sin(10*pi)*exp(-0.3+4^2)+1/2*log2(23)y =2.26183,z = 10sin(π/3) ...cos(π/3)答案:z==10*sin(pi/3)... *cos(pi/3) z =4.33014,x = sin(223/3), y = x^2,z = y*10 ; 求x+2y-5z 答案x=sin(223/3);y=x^2;z=y*10; x+2*y-5*zans =-37.60255,ie z i z i z ⋅==+=π23215,2,72,计算3221z z z z z +=答案:z1=2+7*i,z2=2*i,z3=5*exp(2*pi-i),z=(z1*z2)\z2+z3 z1 =2.0000 + 7.0000i z2 =0 + 2.0000i z3 =1.4466e+003 -2.2530e+003i z =1.4467e+003 -2.2531e+003i6, 建立起始值=3,增量值=5.5,终止值=44的一维数组x (填写程序语句) 答案:=[3:5.5:44]x =3.0000 8.5000 14.0000 19.5000 25.0000 30.5000 36.0000 41.50007, 建立等差一维数组x :首项为0,末项为π,项数为15 (填写程序语句) 答案:x=linspace(0,pi,15) x =Columns 1 through 90 0.2244 0.4488 0.6732 0.8976 1.1220 1.3464 1.5708 1.7952Columns 10 through 152.0196 2.2440 2.4684 2.6928 2.91723.1418, 计算行列式的值187624323=A答案:A=[3 2 3;4 2 6;7 8 1];det(A)ans =-89, 矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=187624323A ,矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=333222111B ;分别求出A x B 及A 与B 中对应元素之间的乘积的程序语句A=[3 2 3;4 2 6;7 8 1];B=[1 1 1;2 2 2;3 3 3]; C=A*B C =16 16 16 26 26 26 26 26 26D=A.*BD =3 2 38 4 1221 24 3答案10,写出计算 Sin(36o)的程序语句.答案:=sin(pi/5)x =0.5878三、MatLab 练习2-矩阵及数组运算1.输入矩阵。

matlab 实验教程 实验一

matlab 实验教程 实验一

MATLAB语言及其应用实验教程实验一实验目的:1.熟悉MATLAB的界面,菜单,会使用DEMO,学会使用帮助(help)。

2.学习MATLAB的基本语法实验内容:1.打开MATLAB,点击各个菜单以了解各个子菜单项。

2.更改当前路径,在指定路径下保存所作的实验。

实验具体步骤:如图所示,在D盘根目录下建立“MATLAB实验”文件夹,并在此文件夹下以各自的学号建立子文件夹。

如图所示,点击工具栏上“Current Directory”栏的浏览按钮。

选择刚刚建立的以各自学号命名的文件夹,“Current Directory”栏显示如下路径格式。

注意!以后所有的实验均保存在此目录中。

3.学习打开和关闭命令窗口(command window),工作空间(workspace)和命令历史窗口(command history)。

实验具体步骤:单击命令窗口右上角的“X”标志,关闭命令窗口;在Desktop菜单下选择“Command window”,打开命令窗口;单击工作空间窗口右上角的“X”标志,关闭工作空间;在Desktop菜单下选择“Workspace”,打开工作空间;单击命令历史窗口右上角的“X”标志,关闭命令历史窗口;在Desktop菜单下选择“Command history”,打开命令历史窗口。

4.练习变量的赋值,包括向量赋值,矩阵赋值以及复数的赋值。

实验具体步骤:变量赋值>> a=100a =100>> b=0.1b =0.1000向量赋值>> a=1:1:10a =1 2 3 4 5 6 7 8 9 10>> b=1:2:10b =1 3 5 7 9矩阵赋值>> a=[1 2 3;4 5 6;7 8 9]a =1 2 34 5 67 8 9复数赋值>> c=3+5.2ic =3.0000 + 5.2000i>> z=[1+2i,3+4i;5+6i,7+8i]z =1.0000 +2.0000i3.0000 +4.0000i5.0000 +6.0000i7.0000 +8.0000i5.用变量检查命令who和whos检查工作空间中的变量。

信号与系统实验(MATLAB版) (1)

信号与系统实验(MATLAB版) (1)

《信号与系统MATLAB实现》实验指导书电气信息工程学院2014年2月长期以来,《信号与系统》课程一直采用单一理论教学方式,同学们依靠做习题来巩固和理解教学内容,虽然手工演算训练了计算能力和思维方法,但是由于本课程数学公式推导较多,概念抽象,常需画各种波形,作题时难免花费很多时间,现在,我们给同学们介绍一种国际上公认的优秀科技应用软件MA TLAB,借助它我们可以在电脑上轻松地完成许多习题的演算和波形的绘制。

MATLAB的功能非常强大,我们此处仅用到它的一部分,在后续课程中我们还会用到它,在未来地科学研究和工程设计中有可能继续用它,所以有兴趣的同学,可以对MATLAB 再多了解一些。

MATLAB究竟有那些特点呢?1.高效的数值计算和符号计算功能,使我们从繁杂的数学运算分析中解脱出来;2.完备的图形处理功能,实现计算结果和编程的可视化;3.友好的用户界面及接近数学表达式的自然化语言,易于学习和掌握;4.功能丰富的应用工具箱,为我们提供了大量方便实用的处理工具;MATLAB的这些特点,深受大家欢迎,由于个人电脑地普及,目前许多学校已将它做为本科生必须掌握的一种软件。

正是基于这些背景,我们编写了这本《信号与系统及MATLAB实现》指导书,内容包括信号的MA TLAB表示、基本运算、系统的时域分析、频域分析、S域分析、状态变量分析等。

通过这些练习,同学们在学习《信号与系统》的同时,掌握MATLAB的基本应用,学会应用MATLAB的数值计算和符号计算功能,摆脱烦琐的数学运算,从而更注重于信号与系统的基本分析方法和应用的理解与思考,将课程的重点、难点及部分习题用MATLAB进行形象、直观的可视化计算机模拟与仿真实现,加深对信号与系统的基本原理、方法及应用的理解,为学习后续课程打好基础。

另外同学们在进行实验时,最好事先预习一些MATLAB的有关知识,以便更好地完成实验,同时实验中也可利用MATLAB的help命令了解具体语句以及指令的使用方法。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

----------------------------------------------------------------------------------------------------- 实验题1:平衡结构的梁受力计算(线性方程组的应用)有一个平面结构如下所示,有13条梁(图中标号的线段)和8个铰接点(图中标号的圈)联结在一起。

其中1号铰接点完全固定,8号铰接点竖直方向固定,并在2号、5号和6号铰接点上,分别有图示的10吨、15吨和20吨的负载。

在静平衡的条件下,任何一个铰接点上水平和竖直方向受力都是平衡的;已知每条斜梁的角度都是45º。

(1)请列出由各铰接点处受力平衡方程构成的线性方程组;
(2)用MATLAB软件求解该线性方程组,确定每条梁受力情况。

目的: 1. 了解MATLAB软件, 学会MATLAB软件的一些基本操作.
2. 练习编写简单的MATLAB程序.
3. 利用线性方程组的知识并借助MATLAB软件解决问题。

要求: 1. 实验报告用A4纸打印, 参考附录的格式. 报告上交的截止时间是本学期结束之前。

实验题2:平板的稳态温度分布问题(线性方程组应用)
在热传导的研究中,一个重要的问题是确定一块平板的稳态温度分布。

假定下图中的平板代表一条金属梁的截面,并忽略垂直于该截面方向上的热传导。

已知平板内部有30个节点,每个节点的温度近似等于与它相邻的四个节点温度的
平均值,例如
4292530r
d T T T T T +++=
;设4条边界上的温度分别等于每位同学学
号的后四位的5倍,例如学号为16308209的同学计算时,选择40=l T 、10=u T 、
0=r T 、45=d T 。

求:(1)建立可以确定平板内节点温度的线性方程组;
(2)用MATLAB 软件求解该线性方程组;
(3)用MATLAB 中的函数mesh 绘制三维平板温度分布图(注:
如果不会绘图,可留待讲解绘图命令之后再做此小题)。

目的: 1. 了解MATLAB 软件, 学会MATLAB 软件的一些基本操作.
2. 练习编写简单的MATLAB 程序.
3. 利用线性方程组的知识并借助MATLAB 软件解决问题。

要求: 1. 实验报告用A4纸打印, 参考附录的格式. 报告上交的截止时间是本学
期结束之前。

实验题4:交通网络流量分析问题(线性方程组应用)
城市道路网中每条道路、每个交叉路口的车流量调查,是分析、评价及改善城市交通状况的基础。

问题:某城市有下图所示的交通图,每条道路都是单行线,需要调查每条道
T 1 T 2 T 3 T 4 T 5
T 6
T 10
T 21 T 25 T 26 T 27
T 28 T 29
T 30
………
………
………
┆┆
┆┆
T l
T u
T d
T r
路每小时的车流量。

图中的数字表示该条路段的车流数。

如果每个交叉路口进入和离开的车数相等,整个图中进入和离开的车数相等。

求(1)建立确定每条道路流量的线性方程组; (2)分析哪些流量数据是多余的;
(3)为了唯一确定未知流量,需要增添哪几条道路的流量统计。

目的: 1. 了解MATLAB 软件, 学会MATLAB 软件的一些基本操作.
2. 练习编写简单的MATLAB 程序.
3. 利用线性方程组的知识并借助MATLAB 软件解决问题。

要求: 1. 实验报告用A4纸打印, 参考附录的格式. 报告上交的截止时间是本学
期结束之前。

附录 实验报告参考模板
注:此页是下一页实验报告模板所对应的实验题目
300
500
150
180
350
160
220
300
100
290
400
150
x 1
x 2
x 3
x 4
x 5
x 6
x 7
x 8
x 9
x 10
x 11
x 12
练习
实验1 求解线性方程组内容: 用MATLAB求解如下线性方程组Ax = b, 其中
A =
56000000
15600000
01560000
00156000
00015600
00001560
00000156
00000015
⎡⎤
⎢⎥
⎢⎥
⎢⎥
⎢⎥
⎢⎥
⎢⎥
⎢⎥
⎢⎥
⎣⎦
, b =
1
2
3
4
5
6
7
8
b
b
b
b
b
b
b
b
⎡⎤
⎢⎥
⎢⎥
⎢⎥
⎢⎥
⎢⎥
⎢⎥
⎢⎥
⎢⎥
⎢⎥
⎣⎦
中的8个分量取自各人的学号.
如: 某同学的学号为02107656, 则该同学取b= [0,2,1,0,7,6,5,6]T做这个数学实验.
目的: 1. 了解MATLAB软件, 学会MATLAB软件的一些基本操作.
2. 熟悉MATLAB软件的一些数值计算功能.
3. 练习编写简单的MATLAB程序.
要求: 1. 用3种不同的方案做这个数学实验.
2. 实验报告用A4纸打印, 参考附录的格式.
实验2 研究三个平面的位置关系
内容: 用MATLAB研究下面的3个平面
π1: x + y + z = 1
π2: -x + y = 2
π3: 2x + t2z = t
当t取何值时交于一点? 当t取何值时交于一直线? 当t取何值时没有公共的交点?
并在每一种情形下, 用MATLAB在同一个坐标系内绘制出这3个平面的图形(其中, 没有公共的交点的情况, 只要给t取一个适当的值并绘制出相应的图形即可).
目的: 1. 练习编写简单的MATLAB程序.
2. 掌握用MATLAB软件绘制简单图形的方法.
要求: 1. 实验报告中要附上所绘制的图形.
2. 实验报告用A4纸打印, 参考附录的格式. 附录实验报告参考模板
数学实验报告
学号: , 姓名: , 得分:
实验1 求解线性方程组
实验内容: 用MATLAB求解如下线性方程组Ax = b, 其中
A =
56000000
15600000
01560000
00156000
00015600
00001560
00000156
00000015
⎡⎤
⎢⎥
⎢⎥
⎢⎥
⎢⎥
⎢⎥
⎢⎥
⎢⎥
⎢⎥
⎣⎦
, b = [……] T.
实验目的: 1..
2..
实验原理: 1. 对于…(满足什么条件)…的方程组Ax = b, 根据…法则, 其解为…
2. 当……时, 方程组Ax = b的解为…
3. 当……时, 对增广矩阵[A, b]进行…, 把它化为…, 则…
实验方案: 1. 在MATLAB命令窗口中输入如下命令:
…… %(该命令的含义/作用)
…… %(该命令的含义/作用)
2. 在MATLAB命令窗口中输入如下命令:
…… %(该命令的含义/作用)
…… %(该命令的含义/作用)
3.
实验结果: 1.
2.
3.
对实验结果的分析:
上述3种方案所得的结果一致/不一致, 这是因为…。

相关文档
最新文档