MATLAB上机练习(3)
MATLAB上机答案
X
1
2
3
4
5
6
7
8
9
10
Y
16
32
70
142 260 436 682 1010 1432 1960
>> x=1:10;y=[16 32 70 142 260 436 682 1010 1432 1960];
>> p1=polyfit(x,y,1)
>>
p1 =
p2=polyfit(x,y,2),y2=polyval(p2,9.5)
0.01
0.005
0
-0.005
-0.01
-0.015
0
1
2
3
4
5
6
7
（3）大气压强 p 随高度 x 变化的理论公式为
，为验证这一公式，
测得某地大气压强随高度变化的一组数据如表所示。试用插值法和拟合法进行计算并绘图，
看那种方法较为合理，且总误差最小。
高度/m
0
300
600
1000
1500
2000
压强/Pa
equally spaced points between X1 and X2. 以 X1 为首元素，X2 为末元素平均生成 100 个元素的行向量。
LINSPACE(X1, X2, N) generates N points between X1 and X2. For N < 2, LINSPACE returns X2.
ans =
pi =
ans =
5
0
3.1416
答：3 次执行的结果不一样。exist（）函数是返回变量搜索顺序的一个函数。在第一次
MATLAB上机作业提交3
上机练习作业-3 图形功能1．函数22yxexz--⨯=，定义区域为[-2，2]×[-2，2]。
生成网格并计算其网格点上的函数值（meshgrid），并用网格图绘制出来。
>> [X,Y]=meshgrid(-2:0.2:2,-2:0.2:2)Z=X.*exp(-X.^2-Y.^2)mesh(Z)2．在0≤X≤2π区间内，用不同颜色和线型在同一坐标内绘制曲线y1=e-1.5x sin(2πx)+1及曲线y2=2.8e-0.3x sin(x)。
给图形添加图形标注（图形标题，X，Y轴说明，曲线说明和图例）及网络线。
>> x=0:0.1:2.*piy1=exp(-1.5*x).*sin(2*pi*x)+1 y2=2.8*exp(-0.3*x).*sin(x) plot(x,y1,':hb',x,y2,'--g') xlabel('x 轴') ylabel('y 轴') title('两个函数') legend('y1','y2',0)1234567x 轴y 轴两个函数3．在一个图形窗口中以子图形式，分别以条形图、填充图、阶梯图和杆图形式绘制曲线y=3.2e -1.75x^2+3x-1，加上必要的标题和说明。
>> x=0:1:5y=3.2.*exp(-1.75.*x.^2+3.*x-1);subplot(2,2,1);bar(x,y,'g'); %subplot 最后一个1表示正在绘制第一个图像 title('bar(x,y,''g'')');axis([0,6,0,7]); subplot(2,2,2);stairs(x,y,'b');title('stairs(x,y,''b'')');axis([0,6,0,7]); subplot(2,2,3);stem(x,y,'k');title('stem(x,y,''k'')');axis([0,6,0,7]); subplot(2,2,4);fill(x,y,'y');title('fill(x,y,''y'')');axis([0,6,0,7]);012345bar(x,y,'g')0246246stairs(x,y,'b')stem(x,y,'k')0246246fill(x,y,'y')4．用分别用mesh 、surf 、plot3表现函数z=sin(y)cos(x) 三维曲面图；裁掉三维曲面图中z>0.25部分，在同一图形窗口表示。
Matlab 课后上机练习3-Matlab绘图和符号运算
11、绘制曲线13++=x x y ，x 的取值范围为[-5,5]。
clear; x=-5:0.1:5; y=x.^3+x+1; plot(x,y,'k'); title('曲线图像'); xlabel('x') ylabel('y') grid on;hold on;2、有一组测量数据满足-at e =y ，t 的变化范围为0~10，用不同的线型和标记点画出a=0.1、a=0.2和a=0.5三种情况下的曲线。
并添加标题：运动曲线图；添加横坐标：时间 t/s ；添加纵坐标：位移 s/mm ；添加图例。
t=0:0.5:10; y1=exp(-0.1*t); y2=exp(-0.2*t); y3=exp(-0.5*t);plot(t,y1,':*r',t,y2,'-^g',t,y3,'-ob') title('运动曲线图'); xlabel('时间 t/s') ylabel('位移 s/mm') legend('a=0.1','a=0.2','a=0.5')3、22y xxe z --=，当x 和y 的取值范围均为-2到2时，用建立子窗口的方法在同一个图形窗口中绘制出三维线图、网线图、表面图和等高线效果图。
[x,y]=meshgrid([-2:0.2:2]);z=x.*exp(-x.^2-y.^2); mesh(x,y,z) subplot(2,2,1) plot3(x,y,z)title('plot3(x,y,z)') subplot(2,2,2) mesh(x,y,z)title('mesh(x,y,z)') subplot(2,2,3) surf(x,y,z)title('surf(x,y,z)') subplot(2,2,4) surf(x,y,z) shading interptitle('surf(x,y,z) shading interp')shading interp4、在同一坐标内绘制如下曲线：（1）y1=tsin(t)（红色连续线‘—’）;（2）y2=t2-cos(t);（蓝色间断线‘—.’）（3）题头：小车运动学分析曲线；图例：y1曲线、y2曲线（4）x轴：时间t/s；y轴：位移曲线/mm （5）曲线上标注文字说明：该运动曲线良好。
Matlab上机作业
一、 MATLAB 基础实验记录每项练习内容及结果1、最简单的计算器使用法，记录操作过程。
2、复数或复矩阵运算，记录操作过程。
3、利用帮助功能查找一个函数的帮助信息。
记录你都看到了什么。
4、察看一个Matlab 工具箱，记录你都看到了什么。
5、写一段100字，对Matlab 的心得体会。
二、 Matlab 编程实验1、有一组学生的考试成绩（见表），根据规定，成绩在100分时为满分，成绩在90~99之间时为优秀，成绩在80~89分之间时为良好，成绩在60~79分之间为及格，成绩在60分以下时为不及格，编制一个根据成绩划分等级的程序。
2、编程，可输入有10个实数的数组，调用函数f 求他们的平均值V ，然后输出其中与V 之差的绝对值为最小的那个数组元素。
3、 若一个数出现在自己平方数的右边，则称该数为同构数。
如，因有5*5=25，25*25=625，所以5和25都是同构数；请编制程序找出1～100之间的全部同构数。
4、 编程序实现功能：在一个元素值按升序存放的整型数组中插入一个数，使得插入后的数组元素仍然有序5、编程序实现功能：一个正整数与3的和是5的倍数，与3的差是6的倍数，求出符合此条件的最小正整数。
6、编程序实现功能：求出所有各位数字的立方和等于1099的3位整数。
7、编制写一个函数实现，一个给定二维（3x3）数组转置，（禁止直接用Matlab 命令） 8、编制写一个函数实现，用“起泡法” 排列输入的10个数字。
9、编制写一个函数实现，输入16进制，输出10进制，（禁止直接用Matlab 命令） 10、编制写一个函数实现，输入一个3x4矩阵，求其中最大元素和最小元素，（禁止直接用Matlab 命令）三、 Matlab 应用实验1、编程绘函数211xy -=的曲线图，展示其定义域内的特性。
2、编程绘函数2112++-=x xy 的曲线图，展示其定义域内的特性。
3、编程绘函数2312+-=x x y 的曲线图，展示其定义域内的特性。
数学建模作业题+答案
数学建模MATLAB 语言及应用上机作业11. 在matlab 中建立一个矩阵135792468101234501234A ⎡⎤⎢⎥⎢⎥=⎢⎥-----⎢⎥⎣⎦答案：A = [1,3,5,7,9;2,4,6,8,10;-1,-2,-3,-4,-5;0,1,2,3,4]2. 试着利用matlab 求解出下列方程的解（线性代数22页例14）123412423412342583692254760x x x x x x x x x x x x x x +-+=⎧⎪--=⎪⎨-+=-⎪⎪+-+=⎩ 答案：A=[2 ,1,-5,1;1,-3,0,-6;0,2,-1,2;1,4,-7,6]; B=[8;9;-5;0]; X=A\B 或A=[2,1,-5,1;1,-3,0,-6;0,2,-1,2;1,4,-7,6] b=[8,9,-5,0]' X=inv(A)*b3. 生成一个5阶服从标准正态分布的随机方阵，并计算出其行列式的值，逆矩阵以及转置矩阵。
答案:A=randn(5) det(A) inv(A) A'4. 利用matlab 求解出110430002A -⎡⎤⎢⎥=-⎢⎥⎢⎥⎣⎦的特征值和特征向量。
答案：A=[-1,1,0;-4,3,0;0,0,2] [V,D]=eig(A)5.画出衰减振荡曲线3sin3t y et -=在[0,4]π上的图像。
要求，画线颜色调整为黑色，画布底面为白色。
（在实际中，很多打印机时黑白的，因此大多数作图要考虑黑白打印机的效果。
） 给出恰当的x ，y 坐标轴标题，图像x 轴的最大值为4π。
6. 生成一个0-1分布的具有10个元素的随机向量，试着编写程序挑选出向量中大于0.5的元素。
数学建模和Matlab 上机作业2（2016-9-20）跟老师做（不用整合进作业中）：上机演示讲解：函数，递归的两个例子的写法。
附：1. Fibonacci Sequence （斐波那契数列）在数学上，费波那西数列是以递归的方法来定义： F1= 1；F2= 1；F （n ）=F （n-1）+F （n-2） 2. 阶乘举例：数学描述：n!=1×2×……×n ；计算机描述：n!=n*(n-1)!自己做（需要整合进作业中，提交到系统中）：1. 写一个m 文件完成分值百分制到5分制的转换（即输入一个百分制，转换后输出一个5级对应的得分，联系条件控制语句）。
matlab实验3：多项式运算
代数多项式求值
y = polyval(p,x)
计算多项式 p 在 x 点的值
注：若 x 是向量或矩阵，则采用数组运算 (点运算)！ 例：已知 p(x)=2x3-x2+3,分别取 x=2 和一个 22 矩阵，
求 p(x) 在 x 处的每个分量上的值
>> p=[2,-1,0,3]; >> x=2; y = polyval(p,x) >> x=[-1,2;-2,1]; y = polyval(p,x)
例：解方程组
x
2yz xz3
2
x 3y 8
>> A=[1 2 -1; 1 0 1; 1 3 0]; >> b=[2;3;8]; >> x=linsolve(A,b)
b是列向量！
非线性方程的根
Matlab 非线性方程的数值求解
fzero(f,x0)：求方程 f=0 在 x0 附近的根。
符号求解
solve 也可以用来解方程组 solve( f1 , f2 , ... , fN , v , ... , fN 确定的方程组关于 v1 , v2 , ... , vN 的解
例：解方程组
x 2 y z 27
x
z
3
x2 3 y2 28
例：2x3-x2+3 <-> [2,-1,0,3]
特别注意：系数中的零是不能省的！
多项式的符号形式：poly2sym 如，>> poly2sym([2,-1,0,3])
运行结果：ans = 2*x^3-x^2+3
多项式四则运算
多项式加减运算
多项式的加减运算就是其所对应的系数向量的加减运算
Matlab上机题库及详细答案_Tonyxie
ax 2 + bx + c 0.5 ≤ x < 1.5 = y a sin c b + x 1.5 ≤ x < 3.5 c ln b + 3.5 ≤ x < 5.5 x
clc;clear; a=input('Please input a= '); b=input('Please input b= '); c=input('Please input c= '); x=input('Please input x= '); disp('#if#'); if x>=0.5&x<1.5 y=a*x^2+b*x+c elseif x<3.5&x>=1.5 y=a*sin(b)^c+x elseif x>=3.5&x<5.5 y=log(abs(b+(c./x))) else y='ERROR!' end disp('#switch#') switch(x*10/5) case {1 2} y=a*x^2+b*x+c case {3 4 5 6} y=a*sin(b)^c+x case {7 8 9 10} y=log(abs(b+(c./x))) otherwise y='ERROR!' end
பைடு நூலகம்
8、数值与符号计算
ex （1）求极限 lim x →+∞ a + be x
（3）已知线性方程组 Ax=b，其中
（2）求不定积分
∫ xe
ax
2014 第一轮上机练习题答案matlab程序设计练习
练习1 基础练习一、矩阵及数组操作：1．利用基本矩阵产生3×3和15×8的单位矩阵、全1矩阵、全0矩阵、均匀分布随机矩阵（[-1，1]之间）、正态分布矩阵（均值为1，方差为4）。
解eye(3) eye(15,8) ones(3) ones(15,8) zeros(3) zeros(15,8)X1=(1+(1-(-1))*rand(3))X2=(1+(1-(-1))*rand(15，8))Y1=1+sqrt(4)*randn(3)Y2=1+sqrt(4)*randn(15，8)2．利用fix及rand函数生成[0，10]上的均匀分布的10×10的整数随机矩阵a，然后统计a中大于等于5的元素个数。
解：a=fix((10-0+1)*rand(10)+0)x=5i=find(a>=x)n=length(i)3．在给定的矩阵中删除含有整行内容全为0的行，删除整列内容全为0的列。
解：A=[1 2 2 3 4;0 0 0 0 0;4 5 6 2 3;0 0 0 0 0;]；A([2,4],:)=[]B=A'；B(:,[2,4])=[]二、绘图：4．在同一图形窗口画出下列两条曲线图像：y1=2x+5；y2=x^2-3x+1，并且用legend标注。
解：x=0:0.01:10;y1=2*x+5;y2=x.^2-3*x+1;plot(x,y1,x,y2)legend('y1', 'y2')5．画出下列函数的曲面及等高线：z=x^2+y^2+sin(xy)．解：[x,y]=meshgrid(0:0.25:4*pi); z=x.^2+y.^2+sin(x.*y)；contour3(x,y,z);meshc(x,y,z)三、程序设计：6．编写程序计算（x在[-3，3]，间隔0.01）解：x=input('ÇëÊäÈëxµÄÖµ:');if x>=-3&x<-1y=(-x.^2-4*x-3)/2;elseif x>=-1&x<1y=-x.^2+1;elseif x>=1&x<=3y=(-x.^2+4*x-3)/2;elsey='error'endy7．有一列分数序列：求前15项的和。
Matlab期末考试上机题汇总(选做)
1.求解方程组：2.用magic 产生一个5*5的矩阵，将这个矩阵的第2列与第5列互换位置。
3.编写一个M 函数，求和 S= 1+2+3+…+10004.定义一个矩阵，放入相应的数据，然后保存到MAT 文件中，清除内存后，然后再调入内存。
查看数据是否改变。
5.根据下面的方程绘制图形，y=cos(x)+x ，x 取值-pi-pi ，然后用多项式拟合函数进行拟合。
6.创建一个5行5列的magic 矩阵，删除其中的第二、三行与第三、四列。
使其变成3行3列的矩阵。
7.建一个2行2列的元胞数组，a={zeros(2,2,2),’hello ’;1758,int16(1:10)},请在此基础上创建一个新的8行8列元胞数组。
⎪⎩⎪⎨⎧=++-=+-=++-4436322321321321x x x x x x x x x8.要求编写一个函数文件，将一个给定图片的所有颜色都变深（变深多少自己斟酌）,其中黑色为0，白色为2559.根据下面的方程绘制图形，y=cos(2x)+sin(x),然后用5次多项式进行拟合。
10.要求编写一个M 函数，完成下列函数的求解：y ＝sin （x ）；y ＝sin （x ＋pi ）；y ＝sin （x ＋pi ／3）；y ＝sin （x ＋pi/2）自变量x 的变化范围从－pi 到pi ，要求四个图像画在同一个图形窗体中。
11.矩阵的使用多项式107423--+x x x 在[-1 ,3 ]之间求值并画出相应的曲线。
12.元胞和结构的定义创建一个元胞数组A{1,1}=’James Bond ’ ;A{1,2}=[1 2 ;3 4 ;5 6] ;A{2,1}=pi ;A{2,2}=zeros(5) ;请读取A 元胞数组中的数字5。
13.M 文件的编写计算分段函数 ⎪⎩⎪⎨⎧≠+=+++=10101)1cos(2x x x x x x x y14.函数文件的编写编写函数文件，求半径为r 的圆的周长和面积。
MATLAB实验上机练习(三)
一、在matlab 的命令窗口里完成如下计算，其中t 的值分别取-1,0,1，表达式如下：4/3)t y e π-=t=[-1 0 1]y=sqrt(3)/2*exp(-4.*t).*sin(4*sqrt(3.*t)+pi/3)t =-1 0 1y =1.0e+04 *2.0897 + 1.2065i 0.0001 0.0000二、自行产生一个5行5列的数组，分别得到最中间的三行三列矩阵、右下角2行2列矩阵，奇数行矩阵、奇数列矩阵、奇数行奇数列矩阵。
>> t=rand(5)%生成矩阵A=t(2:4,2:4)%中间三行散列矩阵B=t(4:5,4:5)%右下角两行两列矩阵C=t(1:2:end,:)%奇数行矩阵D=t(:,1:2:end)%奇数列矩阵E=t(1:2:end,1:2:end)%奇数行列矩阵t =0.1622 0.6020 0.4505 0.8258 0.1067 0.7943 0.2630 0.0838 0.5383 0.96190.3112 0.6541 0.2290 0.9961 0.0046 0.5285 0.6892 0.9133 0.0782 0.7749 0.1656 0.7482 0.1524 0.4427 0.8173A =0.2630 0.0838 0.53830.6541 0.2290 0.99610.6892 0.9133 0.0782B =0.0782 0.77490.4427 0.8173C =0.1622 0.6020 0.4505 0.8258 0.1067 0.3112 0.6541 0.2290 0.9961 0.00460.1656 0.7482 0.1524 0.4427 0.8173D =0.1622 0.4505 0.10670.7943 0.0838 0.96190.3112 0.2290 0.00460.5285 0.9133 0.77490.1656 0.1524 0.8173E =0.1622 0.4505 0.10670.3112 0.2290 0.00460.1656 0.1524 0.8173三、用magic产生一个5*5的矩阵，将这个矩阵的第二行与第三行互换位置。
Matlab上机题库及详细答案
解：c=input('请输入一个字符','s');
if c>='A' & c<='Z'
disp(setstr(abs(c)+abs('a')-abs('A')));
elseif c>='a'& c<='z'
disp(setstr(abs(c)- abs('a')+abs('A')));
(2)均值为0.6方差为0.1的5阶正态分布随机矩阵
解：>>x=20+(50-20)*rand(5)；
>>y=0.6+sqrt(0.1)*randn(5)
例13：将101~125等25个数填入一个5行5列的表格中，使其每行每列及对角线的和均为565。
解：M=100+magic(5)
M =117 124 101 108 115
-0.4606
0.3848
例2、用简短命令计算并绘制在0x6范围内的sin(2x)、sinx2、sin2x。
解：x=linspace(0,6)
y1=sin(2*x),y2=sin(x.^2),y3=(sin(x)).^2;
plot(x,y1,x, y2,x, y3)
例3：画出指数衰减曲线y1=exp(-t/3)*sin(3*t)和它的包络y2=exp(-t/3),t的取值范围是（0，4pi）。
MATLAB上机模拟考试题(一)
(1)在MATLAB的命令窗口中执行_____命令，将命令窗口的显示内容清空。
（） A.clearB.clcC.echo offD.cd(2)在MATLAB的命令窗口中执行_____命令，使数据输出显示为十六进制表示。
（）A.format longB.format ratC.format hexD.format short e(3)下列变量名中_____是合法的。
（） A.x*y,a,1 B.x\y,a1234 C.end,1 bcx D.char_1,i,j(4)已知x=0:5,则x有_____个元素。
（） A.5 B.6 C.7 D.8(5)一下运算符中哪个的优先级最高_____。
（） A./ B.^ C.~= D.&(6)使用检测函数isnumeric(10)的结果是_____。
（） A.1 B.0 C.false D.true(7)三维图形中默认视角是_____度。
（） A.方位角=0 俯仰角=90 B.方位角=90 俯仰角=0 C.方位角=37.5 仰俯角=30 D.方位角=0 仰俯角=180(8)将符号表达式化简为因式分解因式分解因式分解因式分解形式，使用_____函数。
（）A.collectB.expandC.hornerD.factor(9)运行以下命令，则_____描述是正确的。
（）>>syms a b c d >>A=[a b;c d]A.A占用的内存小于100BB.创建了5个符号变量C.A占用的内存是a b c d的总和D.不存在（10）已知数组a=[1 2 3;4 5 6;7 8 9],则a(:,end)是指_____元素。
（11）运行命令bitor（8,7）的结果是_____。
（12）运行以下命令：>>x=0:10; >>y1=sin(x); >>y2=5*sin(x); >>y3=[10*sin(x );20*sin(x)]; >>plot(x,y1,x,y2,x,y3)则在一个图形窗口中，可以看到_____条曲线。
天津工业大学matlab上机实验题答案3
“MATLAB 及其在通信中的应用” 上机实验三1————MATLAB 程序设计初步 学号 班级1、 编写.m 程序。
打开2*2的图形子窗口，按要求绘制图形：关于z 的函数有22-x y z e --=，其中[22],[22]x y ∈-∈-，子窗口1: 绘制三维螺旋曲线图，曲线颜色红色，线型为 * ，线宽为〔提示：plot3〕子窗口2：绘制三维网格曲线，加入舞台效果子窗口3：绘制三维网格曲面，且绘制等高线子窗口4：绘制三维网格曲面，函数外表颜色连续将所编程序和输出图形截图，填入到下面空白处。
〔20分〕答案：1 完成后以附件形式发送到邮箱commu_matlab2014@sohu 。
邮件主题为“班级”加“下划线”加“”加“下划线”加“实验一”，如“通信1101__实验一”，word 文件名类似，如“通信1101__实验一”。
2、 绘制球体图形，要求球体中心坐标为（4,5,6），球体半径为r=10，球体外表有15*15个小网格。
将程序写在下面空白处即可。
无需对绘制出的球体图形截图。
〔10分〕答案：[x,y,z]=sphere(15);X=x*10+4;Y=y*10+5;Z=z*10+6;mesh(X,Y,Z);3、 编写.m 程序。
打开2*2的图形子窗口，按要求绘制图形：子窗口1：绘制b=sin (a)+3的二维散点图，其中π∈a [0 2]。
指定竖条末端圆圈的填充颜色设置为绿色。
横轴对应的变量为a ，纵轴对应的变量为b子窗口2：绘制b=sin (a)+3的梯形图，其中π∈a [0 2]。
曲线颜色黑色，线型为 - ，线宽为 子窗口3：绘制z=sin(x)+cos(y)的三维散点图，其中x 1,2,3,...,10y ==。
指定竖条末端圆圈的填充颜色设置为红色。
用view 函数将方位角设置为-30。
仰角设置为40。
子窗口4：绘制22-x +y z=**e x y （）的二维等高线，其中x [33]∈-、y [33]∈-。
Matlab上机作业部分参考答案
3. 设A为 数组，B为一个行数大于3的数组，请给出 （1）删除A的第4、8、12三列的命令； （2）删除B的倒数第3行的命令； （3）求符号极限 （4）求 的3阶导数
lim tan( mx ) 的命令集； x 0 nx x3 y arctan ln(1 e 2 x ) 的命令集； x2
1 0.8 0.6 0.4 0.2 0 -0.2 -0.4 -0.6 -0.8
0
2
4
6
8
10
12
14
某校60名学生的一次考试成绩如下: 93 75 83 93 91 85 84 82 77 76 77 95 94 89 91 88 86 83 96 81 79 97 78 75 67 69 68 84 83 81 75 66 85 70 94 84 83 82 80 78 74 73 76 70 86 76 90 89 71 66 86 73 80 94 79 78 77 63 53 55 1)计算均值、标准差、极差、偏度、峰度，画出直方图； 2)检验分布的正态性； 3)若检验符合正态分布，估计正态分布的参数并检验参数。 解答： x=[93 75 83 93 91 85 84 82 77 76 77 95 94 89 91 88 86 83 96 81 79 97 78 75 67 69 68 84 83 81 75 66 85 70 94 84 83 82 80 78 74 73 76 70 86 76 90 89 71 66 86 73 80 94 79 78 77 63 53 55]; mean(x) std(x) range(x) skewness(x) kurtosis(x) hist(x) h=normplot(x) [muhat,sigmahat,muci,sigmaci]=normfit(x) [H,sig,ci]=ttest(x,80.1)
matlab上机实验答案
三、假设已知矩阵A ，试给出相应的MATLAB 命令，将其全部偶数行提取出来，赋给B 矩阵，用magic(8)A =命令生成A 矩阵，用上述命令检验一下结果是不是正确。
>> A=magic(8)>> B=A(2:2:end,:)五、选择合适的步距绘制出下面的图形。
（1）)/1sin(t ，其中)1,1(-∈t ；（2）)tan(sin )sin(tan t t -，其中),(ππ-∈t 。
1．>> t=[-1:0.0001:1];y=sin(1./t);plot(t,y)Warning: Divide by zero.>>2．>> t=[-pi:0.001:pi];y=sin(tan(t))-tan(sin(t));plot(t,y)>>七、试求出如下极限。
（1）x x x x 1)93(lim +∞→；（2）11lim 00-+→→xy xy y x ；（3）22)()cos(1lim 222200y x y x e y x y x +→→++-。
（1）>>symsx;f=(3^x+9^x)^(1/x);limit(f,x,inf)ans =9（2）>>syms x y;f=x*y/(sqrt(x*y+1)-1);limit(limit(f,x,0),y,0)ans =2（3）>>syms x y;f=(1-cos(x^2+y^2))/(x^2+y^2)*exp(x^2+y^2);limit(limit(f,x,0),y,0)ans =0九、假设⎰-=xyt t e y x f 0d ),(2，试求222222y f y x f x f y x ∂∂+∂∂∂-∂∂。
>>syms x y t;f=int(exp(-t^2),t,0,x*y);F=x./y.*diff(f,x,2)+2.*diff(diff(f,x,1),y,1)+diff(f,y,2)F =-6*x^2*y^2*exp(-x^2*y^2)+2*exp(-x^2*y^2)-2*x^3*y*exp(-x^2*y^2)ans =-2*exp(-x^2*y^2)*(3*x^2*y^2-1+x^3*y)十一、试求出以下的曲线积分。
Matlab上机实验答案
Matlab 上机实验答案实验一MATLAB运算基础1.先求以下表达式的值，尔后显示 MATLAB 工作空间的使用状况并保存所有变量。
2sin850(1)z11 e2>>z1=2*sin(85*pi/180)/(1+exp(2)) z1 =(2) z21 1 x2212iln( x) ，此中x5 2>>x=[2 1+2i;-0.45 5];>>z2=1/2*log(x+sqrt(1+x^2))z2 =(3) z3e0.3a e sin(a0.3) ln0.3 a22>>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) (>> z33=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2)可以考据 z3==z33 ，能否都为 1）z3 =Columns 1 through 5Columns 6 through 10Columns 11 through 15Columns 16 through 20 Columns 21 through 25 Columns 26 through 30 Columns 31 through 35 Columns 36 through 40Columns 41 through 45 Columns 46 through 50 Columns 51 through 55 Columns 56 through 60Column 61t 20 t1(4) z4t 211t 2 ，此中t 22t 1 2t3>>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)z4 =002.已知：12344131A 34787, B2033657327求以下表达式的值：(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]>>A=[12 34 -4;34 7 87;3 65 7];>>B=[1 3 -1; 2 0 3;3 -2 7];>>A+6*Bans =1852-10467105215349>>I=eye(3);>>A-B+Ians =1231-3328840671(2)>> A*Bans =684462309-72596154-5241>>A.*Bans =1210246802619-13049(3)>> A^3ans =372262338244860424737014918860076678688454142118820>>A.^3ans =172839304-643930434365850327274625343 (4)>> A/Bans =>>B\Aans =(5)>>[A,B]ans =1234-413-13478720336573-27 >>[A([1,3],:);B^2]ans =1234-436574511101920-5403.设有矩阵 A和B1234530166789101769A 1112131415, B02341617181920970212223242541311(1)求它们的乘积 C。
Matlab上机作业部分参考答案
上机练习二 参考答案
1. 产生一个1x10的随机矩阵，大小位于（-5 5），并 且按照从大到小的顺序排列好！ 【求解】 a=10*rand(1,10)-5; b=sort(a,'descend')
上机练习二 参考答案
2、用MATLAB 语句输入矩阵A 和B
前面给出的是4 ×4 矩阵，如果给出A(5，6) = 5 命令，矩阵A将得出什么 结果？
Matlab 上机课作业
吴梅红 2012.10.15
上机练习一
上机练习一 参考答案
上机练习一 参考答案
上机练习一 参考答案
上机练习二
1. 产生一个1x10的随机矩阵，大小位于（-5 5），并且按 照从大到小的顺序排列好！ 2、用MATLAB 语句输入矩阵A 和B
前面给出的是4 ×4 矩阵，如果给出A(5，6) = 5 命令，矩阵 A将得出什么结果？ 3、假设已知矩阵A ，试给出相应的MATLAB 命令，将其全 部偶数行提取出来，赋给B 矩阵，用A =magic(8) 命令生成A 矩阵，用上述的命令检验一下结果是不是正确。
【求解】用课程介绍的方法可以直接输入这两个矩阵 >> A=[1 2 3 4; 4 3 2 1; 2 3 4 1; 3 2 4 1] A= 1234 4321 2341 3241 若给出A(5,6)=5 命令，虽然这时的行和列数均大于A矩阵当前的维数， 但仍然可以执行该语句，得出 >> A(5,6)=5 A= 123400 432100 234100 324100 000005 复数矩阵也可以用直观的语句输入 3+2i 4+1i; 4+1i 3+2i 2+3i 1+4i; 2+3i 3+2i 4+1i 1+4i; 3+2i 2+3i 4+1i 1+4i]; B= 1.0000 + 4.0000i 2.0000 + 3.0000i 3.0000 + 2.0000i 4.0000 + 1.0000i 4.0000 + 1.0000i 3.0000 + 2.0000i 2.0000 + 3.0000i 1.0000 + 4.0000i 2.0000 + 3.0000i 3.0000 + 2.0000i 4.0000 + 1.0000i 1.0000 + 4.0000i 3.0000 + 2.0000i 2.0000 + 3.0000i 4.0000 + 1.0000i 1.0000 + 4.0000i
下MATLAB言语上机练习内容
注意：每个实验奇数题做实验，偶数题做作业（做实验时一同 带往实验室作为实验预习题提交）
1、指令>>A= [2 2 ; 2 2]; >>A .^2= >>A^2= A=[2 2 ;2 2 ];
A.^2 A^2
ans = 44 44
ans = 88 88
2、数组>>x=10:-2:2; y=1:5; >>(x<5)&(y>2)，ans= >> y==2|x~=2， ans=
准
symsum(f,k,0,inf)
ans = piecewise([abs(x + 1)*(abs(x - 1) - abs(x + 1)) < 0, 2*atanh((x 1)/(x + 1))])
（2）求
k0
xk
的符号解，并进而用该符号解求
确值。 （提示：注意 subs 的使用）
s1=symsum(x^k,k,0,inf) s2=symsum((-1/3)^k,k,0,inf) s3=symsum((1/pi)^k,k,0,inf) s4=symsum(3^k,k,0,inf)
syms t f=abs(sin(t)) s1=diff(f,t) s2=limit(s1,t,0,'left') s3=limit(s1,t,pi/2)
f= abs(sin(t)) s1 = abs(1,sin(t))*cos(t) s2 = -1 s3 =
对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术关，通系电1，力过根保管据护线生高0不产中仅工资2艺料22高试2可中卷以资配解料置决试技吊卷术顶要是层求指配，机置对组不电在规气进范设行高备继中进电资行保料空护试载高卷与中问带资题负料2荷试2，下卷而高总且中体可资配保料置障试时2卷，32调需3各控要类试在管验最路；大习对限题设度到备内位进来。行确在调保管整机路使组敷其高设在中过正资程常料1工试中况卷，下安要与全加过，强度并看工且25作尽52下可22都能护可地1关以缩于正小管常故路工障高作高中；中资对资料于料试继试卷电卷连保破接护坏管进范口行围处整，理核或高对者中定对资值某料，些试审异卷核常弯与高扁校中度对资固图料定纸试盒，卷位编工置写况.复进保杂行护设自层备动防与处腐装理跨置，接高尤地中其线资要弯料避曲试免半卷错径调误标试高方中等案资，，料要编试求5写、卷技重电保术要气护交设设装底备备置。4高调、动管中试电作线资高气，敷料中课并设3试资件且、技卷料中拒管术试试调绝路中验卷试动敷包方技作设含案术，技线以来术槽及避、系免管统不架启必等动要多方高项案中方；资式对料，整试为套卷解启突决动然高过停中程机语中。文高因电中此气资，课料电件试力中卷高管电中壁气资薄设料、备试接进卷口行保不调护严试装等工置问作调题并试，且技合进术理行，利过要用关求管运电线行力敷高保设中护技资装术料置。试做线卷到缆技准敷术确设指灵原导活则。。：对对在于于分调差线试动盒过保处程护，中装当高置不中高同资中电料资压试料回卷试路技卷交术调叉问试时题技，，术应作是采为指用调发金试电属人机隔员一板，变进需压行要器隔在组开事在处前发理掌生；握内同图部一纸故线资障槽料时内、，设需强备要电制进回造行路厂外须家部同出电时具源切高高断中中习资资题料料电试试源卷卷，试切线验除缆报从敷告而设与采完相用毕关高，技中要术资进资料行料试检，卷查并主和且要检了保测解护处现装理场置。设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。
Matlab编程与应用习题和一些参考答案
Matlab 上机实验一、二1．安装Matlab 软件。
2．验证所学内容和教材上的例子。
3.求下列联立方程的解⎪⎪⎩⎪⎪⎨⎧=+-+-=-+=++-=--+41025695842475412743w z y x w z x w z y x w z y x>> a=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10];>> b=[4;4;9;4];>> c=a\bc =5.22264.45701.47181.59944．设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡------=81272956313841A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----=793183262345B ，求C1=A*B’;C2=A’*B;C3=A.*B,并求上述所有方阵的逆阵。
>> A=[1 4 8 13;-3 6 -5 -9;2 -7 -12 -8];>> B=[5 4 3 -2;6 -2 3 -8;-1 3 -9 7];>> C1=A*B'C1 =19 -82 3012 27 3-38 54 29>> C2=A'*BC2 =-15 16 -24 3663 -17 93 -10522 6 117 -6019 46 84 -10>> C3=A.*BC3 =5 16 24 -26-18 -12 -15 72-2 -21 108 -56>> inv(C1)ans =0.0062 0.0400 -0.0106-0.0046 0.0169 0.00300.0168 0.0209 0.0150>> inv(C2)Warning: Matrix is close to singular or badly scaled.Results may be inaccurate. RCOND = 8.997019e-019.ans =1.0e+015 *-0.9553 -0.2391 -0.1997 0.27000.9667 0.2420 0.2021 -0.2732-0.4473 -0.1120 -0.0935 0.1264-1.1259 -0.2818 -0.2353 0.3182>> inv(C3)Error using ==> invMatrix must be square.5．设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y ，把x=0~2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线。