同济matlab选修课第一次实验内容参考答案
matlab、实验一答案
在 command 窗口上输入： >> a=@cos; x=[0,2*pi]; >> fplot(a,x) Enter 后弹出图形窗口：
（ 2 ） 在 同 一 坐 标 系 中 绘 制 余 弦 曲 线 y cos(t 0.25) 和 正 弦 曲 线
y sin(t 0.5), t [0,2 ]
6.基本程序设计 （1）编写命令文件：计算 1 2 n 2000 时的最大 n 值。 function[]=program_6() clear all; sum=0; n=1; i=1;
for i=1:100; sum=sum+i; if sum<2000 n=i; else n=i-1; break; end end
在 command 窗口输入： >> t=0:0.01:2*pi; >> plot(t,cos(t-0.25),'-g') >> hold on; >> plot(t,sin(t-0.5)) Enter 后出现下图结果，其中绿色线条是余弦曲线，蓝色线条是正弦曲线：
5.基本绘图控制，绘制 [0,4 * ] 区间上 x1 10 sin t 曲线，并要求 （1）线形为点划线，颜色为红色，数据点标记为加号。 （2）坐标轴控制：显示范围、刻度线、比例、网格线。 （3）标注控制：坐标轴名称、标题、相应文本。 >> t=0:0.1:4*pi; >> x=10*sin(t); >> plot(t,x,'-.+r'); %设置线形为点划线，数据点为加号，红色 >> axis([0,13,-11,11]) %设置坐标轴显示范围，横轴[0,13]，纵轴[-11,11] >> set(gca,'xtick',[0,0.5*pi,pi,1.5*pi,2*pi,2.5*pi,3*pi,3.5*pi,4*pi]) %设置横轴刻度线 >> set(gca,'yminortick','on'); %设置纵轴刻度线 >> grid on %显示网格线 >> xlabel('t 轴');ylabel('x 轴');title('x=10sint'); %设置坐标轴名称、标题 >> text(2*pi,0,'\leftarrow 10sin(2\pi)','fontsize'2 3;4 5 6;7 8 9] >> a=A(1:2,3) a= 3 6 >> b=A(2:3,:) b= 4 5 6 7 8 9 3.多项式
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． 设要求以0.01秒为间隔，求出y 的151个点，并求出其导数的值和曲线。
clcclearx=0:0.01:1.5;y=sqrt(3)/2*exp(-4*x).*sin(4*sqrt(3)*x+pi/3)y1=diff(y)subplot(2,1,1)plot(x,y)subplot(2,1,2)plot(x(1:150),y1)2绘制极坐标系下曲线（a,b,n 自定数据）clccleara=10;b=pi/2;n=5;theta=0:pi/100:2*pi;rho=a*cos(b+n*theta);polar(theta,rho)3. 列出求下列空间曲面交线的程序clcclearx=[-5:0.5:5];[X,Y]=meshgrid(x);z1=X.^2-2*Y.^2;z2=X.*2-Y.*3;xlabel('x')ylabel('y')zlabel('z')surf(X,Y,z1)hold onsurf(X,Y,z2)k=find(abs(z1-z2)<0.5);x1=X(k)y1=Y(k)z3=x1.^2-2*y1.^2hold onplot3(x1,y1,z3,'*')⎪⎭⎫ ⎝⎛+⋅=-334sin 234πt e y t ()θρn b a +=cos 2212y x z -=y x z 322-=4、设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y 把x=0~2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线，要求有图形标注。
clcclearx=-2*pi:0.1: 2*pi;y=cos(x).*(0.5+sin(x)*3./(1+x.^2));plot(x,y,'b*-');title('绘图');xlabel('x 坐标');ylabel('y 坐标');legend('原函数')gtext('y=cos(x)(0.5+3*sin(x)/(1+x^2))')5、求下列联立方程的解81025695832475412743-=+-+-=-+-=++-=--+w z y x w z x w z y x w z y x clccleara=[3,4,-7,-12;5,-7,4,2;1,0,8,-5;-6,5,-2,10];b=[4,-3,9,-8];c=b/a;x=c(1,1)y=c(1,2)z=c(1,3)w=c(1,4)6. 假设一曲线数据点为x = 0:2:4*pi;y = sin(x).*exp(-x/5);试将x 的间距调成 0.1，采用不同插值方法进行插值，并通过子图的形式将不同插值结果和原始数据点绘制在同一图形窗口。
实验一 MATLAB软件应用复习含具体答案
，其中E,R,O,S分别为单位阵、随机阵、零阵和对角阵，试通过数值计算验证
。 2．某零售店有9种商品的单件进价（元）、售价（元）及一周的销量如表1.1，问哪种商品 的利润最大，哪种商品的利润最小；按收入由小到大，列出所有商品及其收入；求这一周该9种 商品的总收入和总利润。 表1.1
货号
1
2
3. 用两种方法在同一个坐标下作出 这四条曲线的图形，并要求用两种方法在图上加各种标注。 4．用subplot分别在不同的坐标系下作出下列四条曲线，为每幅图形加上标题， 1）概率曲线 ； 2）四叶玫瑰线 ； 3）叶形线
4）曳物线
。 5．作出下列曲面的3维图形， 1） ； 2）环面
（2） 新建一个M文件输入命令如下： v=0:0.01:2*pi; u=v; [U,V]=meshgrid(u,v); X=(1+cos(U)).*cos(V); Y=(1+cos(U)).*sin(V); Z=sin(U); mesh(X,Y,Z) 保存并运行，运行结果如下：
T6 新建一个M文件输入命令如下： function shuixianhuashu for m=100:999; a=fix(m/100); b=rem(fix(m/10),10); c=rem(m,10); if m==a.^3+b.^3+c.^3 disp(m) end end 保存，再在命令窗口输入shuixianhuashu结果如下： >> shuixianhuashu 153
探究实验 自由发挥：自己提出问题，实验探索，广泛联想，发现规律，大胆猜想。比如函数cos(1/x)在 x=0附近的振荡现象，有无规律可寻？
T1
新建一个M文件输入命令如下： E=eye(3,3); R=rand(3,2); O=zeros(2,3); S=diag([3,7]); A=[E,R;O,S]; AA=A*A; B=[E,R+R*S;O,S*S]; if B==AA disp('命题成立'); else disp('命题不成立'); end 保存并运行，结果如下： 命题成立 由运行结果可知，得以验证。
实验一MATLAB基本操作并且附有答案.
实验一、MA TLAB基本操作一、实验目的熟悉MA TLAB软件环境,掌握命令窗口的使用。
二、实验内容及步骤1、命令窗口的简单使用(1简单矩阵的输入1 2 3A = 4 5 67 8 9>> A=[1 2 3 ;4 5 6 ;7 8 9 ](2求[12+2×(7-4]÷32的算术运算结果>>[12+2*(7-4]/3^22、有关向量、矩阵或数组的一些运算(1设A=15;B=20;求C=A+B与c=a+b?(2设A=[1 2 3;4 5 6;7 8 9],B=[9 8 7;6 5 4;3 2 1];求A*B与A.*B? (3设a=10,b=20;求i=a/b与j=a\b?(4设a=[1 -2 3;4 5 -4;5 -6 7]请设计出程序,分别找出小于0的矩阵元素及其位置(单下标、全下标的形式,并将其单下标转换成全下标。
>> a=[1 -2 3;4 5 -4;5 -6 7];>> [i,j]=find(a<0>>b=find(a<0 >> a(b >> [i,j]=ind2sub(size(a,b (5在MA TLAB命令行窗口运行A=[1,2;3,4]+i*[5,6;7,8];看结果如何?如果改成运行A=[1,2;3,4]+i[5,6;7,8],结果又如何?>> >>A=[1,2;3,4]+i*[5,6;7,8]A =1.0000 + 5.0000i2.0000 + 6.0000i3.0000 + 7.0000i4.0000 + 8.0000i当输入A=[1,2;3,4]+i[5,6;7,8]显示出错(6请写出完成下列计算的指令:a=[1 2 3;3 4 2;5 2 3],求a^2=?,a.^2=?>>a=[1 2 3;3 4 2;5 2 3];a^2 代表两个矩阵相乘ans =22 16 1625 26 2326 24 28>> a.^2 代表A矩阵元素的平方ans =1 4 99 16 425 4 9(7有一段指令如下,请思考并说明运行结果及其原因clearX=[1 2;8 9;3 6];X( : 访问向量是列优先。
MATLAB数学实验答案（全）
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 limx mx mx x →-与3sin lim x mx mxx →∞-syms xlimit((902*x-sin(902*x))/x^3) ans =366935404/3limit((902*x-sin(902*x))/x^3,inf)//inf 的意思 ans = 0 1.2 cos1000xmxy e =，求''y syms xdiff(exp(x)*cos(902*x/1000),2)//diff 及其后的2的意思 ans =(46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算221100x y edxdy +??dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1)//双重积分 ans = 2.13941.4 计算4224x dx m x +? syms xint(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 xdiff(exp(x)*cos(902*x),10) ans =-356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x)1.6 0x =的泰勒展式(最⾼次幂为4).syms xtaylor(sqrt(902/1000+x),5,x)//泰勒展式 ans =-(9765625*451^(1/2)*500^(1/2)*x^4)/82743933602 +(15625*451^(1/2)*500^(1/2)*x^3)/91733851-(125*451^(1/2)*500^(1/2)*x^2)/406802 + (451^(1/2)*500^(1/2)*x)/902 +(451^(1/2)*500^(1/2))/500 1.7 Fibonacci 数列{}n x 的定义是121,1x x ==12,(3,4,)n n n x x x n --=+=⽤循环语句编程给出该数列的前20项（要求将结果⽤向量的形式给出）。
Matlab实验1答案
Matlab 曲线绘图练习1. 绘出立方曲线3=。
y xx=-2::2;y=x.^3;plot(x,y)grid on2. 立方抛物线y=y=-2::2;x=y.^3;plot(x,y)grid on3. 高斯曲线2x y e -=。
clear; x=-2::2;y=exp(-x.^2); plot(x,y)以参数方程表示的曲线： 4. 奈尔抛物线2323,()x t y t y x ===clear; t=-2::2;x=t.^3;y=t.^2 plot(x,y)y =Columns 1 through 14Columns 15 through 28Columns 29 through 420 Columns 43 through 56Columns 57 through 70Columns 71 through 815. 半立方抛物线2323,()x t y t y x ===clear; t=-2::2;x=t.^2;y=t.^3 plot(x,y)y =Columns 1 through 14Columns 15 through 28Columns 29 through 420 Columns 43 through 56Columns 57 through 70 Columns 71 through 816. 迪卡尔曲线2332233,(30) 11at atx y x y axyt t==+-= ++clear;a=3;t=-2::2;x=3*a.*t./(1+t.^2);y=3*a.*t.^2./(1+t.^2); plot(x,y);grid on7. 蔓叶线233222,() 11at at xx y yt t a x ===++-clear;a=2;t=-10::10;x=a.*t.^2./(1+t.^2);y=a.*t.^3./(1+t.^2); plot(x,y);grid on8. 摆线(sin),(1cos)x a t t y b t=-=-。
matlab选修课 答案详解 实验指导书Matlab软件应用与开发
《Matlab软件应用与开发》实验指导书实验1 Matlab 操作基础及矩阵运算(一)、实验类型：验证型(二)、实验类别：基础实验(三)、每组人数：1(四)、实验要求：选修(五)、实验学时：3个学时(六)、实验目的：（1）熟悉MATLAB软件中关于矩阵运算的各种命令；（2）学会运用MATLAB软件自定义函数，并求出函数值；（3）学会在MATLAB环境下编写函数。
(七)、预备知识：线性代数中的矩阵运算；高等数学中微积分知识。
本实验所用MATLAB命令●矩阵输入格式：A=[a11 a12;a21 a22];b=初始值:步长:终值●求A的转置：A’●求A加B:A+B●求A减B:A-B●求A乘B:A*B●求A的行列式：det(A)●求A的逆:inv(A)●求A的秩: rank(A)●求函数的极限limit(.)● 求函数的导数diff(.) ● 求函数的积分 int(.) ● 求代数方程的解 solve(.) ● 求微分方程的解 dsolve(.) (八)、内容与要求： 1、 输入矩阵A,B,b;⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=521572215431352134153524852421A ，⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=232556463515148355372182414512B []27252321=b2、 矩阵转置、四则运算。
C1=A ’，C2=A+B,C3=A-B,C4=A*B3、 求行列式。
D1=|A|,D2=|B|4、 求矩阵A 、B 的秩E1,E25、 求极限 11232lim +∞→⎪⎭⎫⎝⎛++x x x x6、 设x xe y =，求)5(y7、 求dx xx ⎰41ln8、求微分方程的通解232++=+'x x y y x9、自定义函数21221221)1()(100),(x x x x x f -+-=，并计算)2,1(f 。
(九)、实验操作1、>> A=[21 24 85 4;2 35 15 34;21 35 31 54;21 72 15 52] A =21 24 85 42 35 15 3421 35 31 5421 72 15 52>> B=[12 45 1 24;18 72 53 35;48 1 15 35;46 56 25 23]B =12 45 1 2418 72 53 3548 1 15 3546 56 25 23>> b=21:2:27b =21 23 25 272、>> C1=A',C2=A+B,C3=A-B,C4=A*BC1 =21 2 21 2124 35 35 7285 15 31 154 34 54 52C2 =33 69 86 2820 107 68 6969 36 46 8967 128 40 75C3 =9 -21 84 -20-16 -37 -38 -1-27 34 16 19-25 16 -10 29C4 =4948 2982 2668 4411 2938 4529 2932 2580 4854 6520 3691 4056 4660 9056 5362 4745 3、>> D1=det(A),D2=det(B)D1 =2181568D2 =-31822764、>> E1=rank(A),E2=rank(B);E1 =4E2 =45、>> syms x; %定义符号变量x>> limit(((2*x+3)/(2*x+1))^(x+1),x,inf) %求函数((2*x+3)/(2*x+1))^(x+1)当x->inf 时的极限 (inf 即英文 infinity “无穷”的缩写) inf表示正无穷大ans =exp(1)6、>> syms y x %定义符号变量>> y=x*exp(x); %定义符号函数>> diff(y,5) %计算符号函数的五阶导数ans =5*exp(x)+x*exp(x)7、>> syms x s %定义符号变量>> s=(log(x))/sqrt(x) %定义符号表达式>> int(s,1,4) %计算符号表达式在区间[1，4]上的定积分ans =8*log(2)-4 8、>> [y]=dsolve('x*Dy+y=x^2+3*x+2','x') %微分或导数的输入是用Dy 、D2y 、D3y 、…来表示y的一阶导数dx dy或y '、二阶导数x d y d 22或y ''、三阶导数x d y d 33或y '''、…。
(完整版)MATLAB)课后实验答案[1]
![(完整版)MATLAB)课后实验答案[1]](https://img.taocdn.com/s3/m/f6c08c5902d276a201292ece.png)
1 + e2 (2) z = 1 ln( x + 1 + x 2 ) ，其中 x = ⎡⎢ 2⎣-0.45 ⎦2 2 ⎪t 2 - 2t + 1 2 ≤ t <3 ⎨实验一MATLAB 运算基础1. 先求下列表达式的值，然后显示 MATLAB 工作空间的使用情况并保存全部变量。
(1) z = 2sin 8501221 + 2i ⎤5 ⎥(3) z = e 0.3a - e -0.3asin(a + 0.3) + ln 0.3 + a ,a = -3.0, - 2.9, L , 2.9, 3.03⎧t 2 0 ≤ t < 1 (4) z = ⎪t 2 - 11 ≤ t <2 ，其中 t=0:0.5:2.5 4⎩解：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;3=(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=⎣O2⨯3⎥，其中E、R、O、S分别为单位矩阵、随机矩阵、零矩S⎦阵和对角阵，试通过数值计算验证A=⎢⎥。
MATLAB)课后实验答案[1]
![MATLAB)课后实验答案[1]](https://img.taocdn.com/s3/m/e0ec95dd10a6f524ccbf85c9.png)
实验一 MATLAB 运算基础1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。
(1) 0122sin 851z e =+(2) 21ln(2z x =+，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦ (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e az a a --+=++=--(4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:: 解：4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。
(2) 建立一个字符串向量，删除其中的大写字母。
)解：(1) 结果：](2). 建立一个字符串向量例如：ch='ABC123d4e56Fg9';则要求结果是：？实验二 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS A OS +⎡⎤=⎢⎥⎣⎦。
解: M 文件如下；5. 下面是一个线性方程组：1231112340.951110.673450.52111456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦(1) 求方程的解。
123d4e56g9(2) 将方程右边向量元素b 3改为再求解，并比较b 3的变化和解的相对变化。
(3) 计算系数矩阵A 的条件数并分析结论。
（解： M 文件如下：实验三 选择结构程序设计1. 求分段函数的值。
2226035605231x x x x y x x x x x x x ⎧+-<≠-⎪=-+≤<≠≠⎨⎪--⎩且且及其他用if 语句实现，分别输出x=,,,,,,时的y 值。
Matlab实验一答案
MATLAB 实验一答案1.3 先自定义一个变量，然后分别用8种不同的数字显示格式显示查看。
>> a=pia =3.1416>> format long>> aa =3.14159265358979>> format short>> aa =3.1416>> format short e>> aa =3.1416e+000>> format long e>> aa =3.141592653589793e+000>> format hex>> aa =400921fb54442d18>> format bank>> aa =3.14>> format +>> aa =+>> format rat>> aa =355/1131.4 下面的语句用于画出函数()0.22x y x e -=在[0,10]区间的值x = 0:0.1:10;y = 2*exp(-0.2*x);plot(x,y)1.5 用Matlab 编辑器创建一个m 文件，把上述语句写入这个m 文件并命名为“test1.m ”，保存在当前路径中，然后在命令窗中键入test1，观察结果和运行程序后工作空间的变化.工作空间和结果均与1.4一样1.6 如何清空工作区间数据？键入 clear ；如何关闭图像窗口？键入close ；除了在命令窗输入文件名，还可以怎样运行一个m 文件程序？点击file ，打开m 文件，点击Run 按钮，运行m 文件程序。
1.7 通过以下两种方式得到关于exp 函数的帮助：(1) 在命令窗中输入help exp 命令；(2) 运用帮助空间窗口。
思考，用什么指令可以直接打开帮助空间中关于exp 函数的说明？键入doc exp1.8 假设x =3，y = 4，用Matlab 计算下列表达式：(1) ()232x y x y - (2) 43x y (3) 24x x π- (4) 33x x x y- >> x=3,y=4;>> x^2*(y^3)/(x-y)^2ans =576>> 4*x/(3*y)ans =1>> 4/x*(pi*x^(-2))ans =0.4654>> x^3/(x^3-y^x)ans =-0.72971.9 在当前目录下创建一个m文件，键入以下程序并保存，先把文件保存为“2.m”，运行后观察结果，总结m文件的文件名（包括Matlab标识符）命名规则。
Matlab实验第一次实验答案
实验一Matlab使用方法和程序设计一、实验目的1、掌握Matlab软件使用的基本方法；2、熟悉Matlab的数据表示、基本运算和程序控制语句3、熟悉Matlab绘图命令及基本绘图控制4、熟悉Matlab程序设计的基本方法二、实验内容：1、帮助命令使用help命令，查找 sqrt（开方）函数的使用方法；解：sqrtSquare rootSyntaxB = sqrt(X)DescriptionB = 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.RemarksSee sqrtm for the matrix square root.Examplessqrt((-2:2)')ans =0 + 1.4142i0 + 1.0000i1.00001.41422、矩阵运算（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）矩阵的转置及共轭转置已知A=[5+i,2-i,1;6*i,4,9-i];求A.', A'解：A=[5+1i,2-1i,1;6*1i,4,9-1i ];A1=A.',A2=A'（4）使用冒号表达式选出指定元素已知： A=[1 2 3;4 5 6;7 8 9];求A中第3列前2个元素；A中所有列第2，3行的元素；方括号[]解：A=[1 2 3;4 5 6;7 8 9 ];B=A([1,2],[3]),C=A(2:end, : )用magic 函数生成一个4阶魔术矩阵，删除该矩阵的第四列 3、多项式（1）求多项式 42)(3--=x x x p 的根解：A=[1 0 -2 -4];B=roots(A)（2）已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ， 求矩阵A 的特征多项式;解：A=[1.2 3 5 .9 ; 5 1.7 5 6 ;3 9 0 1 ;1 2 3 4]; A=poly(A); A=poly2sym(A)把矩阵A作为未知数代入到多项式中；4、基本绘图命令（1）绘制余弦曲线 y=cos(t)，t∈[0，2π]解：t=0:.1:2*pi;y=cos(t);plot(t,y),grid（2）在同一坐标系中绘制余弦曲线y=cos(t-0.25)和正弦曲线y=sin(t-0.5)，t∈[0，2π] t=0:.1:2*pi;y1=cos(t-0.25);y2=sin(t-0.25);plot(t,y1，t,y2)grid5、基本绘图控制绘制[0，4π]区间上的x1=10sint曲线，并要求：（1）线形为点划线、颜色为红色、数据点标记为加号；（2）坐标轴控制：显示范围、刻度线、比例、网络线（3）标注控制：坐标轴名称、标题、相应文本；t=0:.01:4*pi;y=10*sin(t);plot(t,y,'-.',t,y,'r')grid6、基本程序设计（1）编写命令文件：计算1+2+…+n<2000 时的最大n值；解法1：s=0;i=0;while(s<2000) i=i+1;s=s+i; ends=s-i,i=i-1解法2：s=0;for i=1:1000; s=s+i;if(s>2000) ,break;endends=s-i,i=i-1（2）编写函数文件：分别用for和while循环结构编写程序，求2的0到n次幂的和。
MATLAB实验一(课后练习答案)
求： ， ， ， ， ， ，
答：
>> a=[1 2;3 4;5 6],b=[5 6;3 4;1 2],c=[1 2 3;4 5 6]
a =
1 2
3 4
5 6
b =
5 6
3 4
1 2
c =
1 2 3
4 5 6
>> [mA,nA]=size(a),[mB,nB]=size(b)
mA =
3
nA =
-3 6 -5 -9
2 -7 -12 -8
B =
5 4 3 -2
6 -2 3 -8
-1 3 -9 7
>> C1=A*B',C2=A'*B,C3=A.*B
C1 =
19 -82 30
12 27 3
-38 54 29
C2 =
-15 16 -24 36
63 -17 93 -105
22 6 117 -60
19 46 84 -10
2
mB =
3
nB =
2
>>a+b
ans =
6 8
6 8
6 8
>>a-b
ans =
-4 -4
0 0
4 4
>>a*c
ans =
9 12 15
19 26 33
29 40 51
>>a*5
ans =
5 10
15 20
25 30
>>a/2
ans =
0.5000 1.0000
1.5000 2.0000
2.5000 3.0000
(完整版)MATLAB)课后实验答案[1]
![(完整版)MATLAB)课后实验答案[1]](https://img.taocdn.com/s3/m/d27ebaa958fb770bf78a55b8.png)
实验一 MATLAB 运算基础1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。
(1) 0122sin 851z e =+(2) 21ln(2z x =，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦ (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+=++=--L (4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0.5:2.5 解：4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。
(2) 建立一个字符串向量，删除其中的大写字母。
解：(1) 结果：(2). 建立一个字符串向量例如：ch='ABC123d4e56Fg9';则要求结果是：实验二 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS A O S +⎡⎤=⎢⎥⎣⎦。
解: M 文件如下；5. 下面是一个线性方程组：1231112340.951110.673450.52111456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦(1) 求方程的解。
(2) 将方程右边向量元素b 3改为0.53再求解，并比较b 3的变化和解的相对变化。
(3) 计算系数矩阵A 的条件数并分析结论。
解： M 文件如下： 123d4e56g9实验三 选择结构程序设计1. 求分段函数的值。
2226035605231x x x x y x x x x x x x ⎧+-<≠-⎪=-+≤<≠≠⎨⎪--⎩且且及其他用if 语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y 值。
参考答案Matlab实验报告
实验一 Matlab基础知识一、实验目的：1.熟悉启动和退出Matlab的方法。
2.熟悉Matlab命令窗口的组成。
3.掌握建立矩阵的方法。
4.掌握Matlab各种表达式的书写规则以及常用函数的使用。
二、实验内容：1.求[100，999]之间能被21整除的数的个数。
（rem）2.建立一个字符串向量，删除其中的大写字母。
（find）3.输入矩阵，并找出其中大于或等于5的元素。
（find）4.不采用循环的形式求出和式6312ii=∑的数值解。
(sum)三、实验步骤：●求[100，199]之间能被21整除的数的个数。
（rem）1.开始→程序→Matlab2.输入命令：»m=100:999;»p=rem(m,21);»q=sum(p==0)ans=43●建立一个字符串向量，删除其中的大写字母。
（find）1.输入命令：»k=input('’,’s’)；Eie48458DHUEI4778»f=find(k>=’A’&k<=’Z’)；f=9 10 11 12 13»k(f)=[ ]K=eie●输入矩阵，并找出其中大于或等于5的元素。
（find）1.输入命令：»h=[4 8 10;3 6 9; 5 7 3]；»[i,j]=find(h>=5)i=3 j=11 22 23 21 32 3●不采用循环的形式求出和式的数值解。
（sum）1.输入命令：»w=1:63；»q=s um(2.^w)q=1.8447e+019- 1 -实验二 Matlab 基本程序一、 实验目的：1. 熟悉Matlab 的环境与工作空间。
2. 熟悉M 文件与M 函数的编写与应用。
3. 熟悉Matlab 的控制语句。
4. 掌握if,switch,for 等语句的使用。
二、 实验内容：1. 根据y=1+1/3+1/5+……+1/（2n-1）,编程求：y<5时最大n 值以及对应的y 值。
MATLAB)课后实验答案[1]
![MATLAB)课后实验答案[1]](https://img.taocdn.com/s3/m/80dd7d3ad1f34693daef3ec1.png)
实验一 MATLAB 运算基础1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。
(1) 0122sin 851z e =+(2) 21ln(2z x =，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦ (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+=++=--L (4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0.5:2.5 解：4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。
(2) 建立一个字符串向量，删除其中的大写字母。
解：(1) 结果：(2). 建立一个字符串向量例如：ch='ABC123d4e56Fg9';则要求结果是：实验二 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS A OS +⎡⎤=⎢⎥⎣⎦。
解: M 文件如下；5. 下面是一个线性方程组：ch =123d4e56g91231112340.951110.673450.52111456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦(1) 求方程的解。
(2) 将方程右边向量元素b 3改为0.53再求解，并比较b 3的变化和解的相对变化。
(3) 计算系数矩阵A 的条件数并分析结论。
解： M 文件如下：实验三 选择结构程序设计1. 求分段函数的值。
2226035605231x x x x y x x x x x x x ⎧+-<≠-⎪=-+≤<≠≠⎨⎪--⎩且且及其他用if 语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y 值。
同济大学Matlab实验报告(matlab应用实例)
实 验 报 告班级机械三班学号姓名- 1 -- 2 -画出衰减振荡曲线t ey t 3sin 3-=及其它的包络线30t e y -=，31t e y --=。
t 值的范围是');- 3 -4．通过M 脚本文件，画出下列分段函数所表示的曲面，用冷色调。
⎪⎪⎩⎪⎪⎨⎧-≤+≤+<->+=+-------15457.0117575.015457.0),(215.175.375.0216215.175.375.02112122212212122x x e x x ex x e x x p x x x x x x x x [X1,X2]=meshgrid(-1.5:0.1:1.5,-2:0.1:2);P=0.5457*exp(-0.75*X2.^2-3.75*X1.^2-1.5*X1).*(X1+X2>1)... + 0.7575*exp(-X2.^2-6*X1.^2).*(X1+X2>-1&X1+X2<=1)... + 0.5457*exp(-0.75*X2.^2-3.75*X1.^2+1.5*X1).*(X1+X2<=-1); surf(X1,X2,P);colormap(cool);colorbar('horiz'); shading flat;实验体会与总结通过此次实验，把课堂所学的理论知识运用到了实际中，了解了Matlab 的基本功能和用途。
经过4个实验基，本上掌握了Matlab 绘制曲线、曲面和构造分段函数的方法，相信此次实验会为以后Matlab 的学习奠定坚实的基础。
- 1 -实 验 报 告班级 机械三班 学号 姓名的值，输出一元二次方程2ax bx c ++=root方程形式:a*x^2+b*x+c=0 请输入各项系数: a=1 b=1 c=-1ans =-1.6180 0.6180exchange 请输入x ：1 请输入y ：0 x=0, y=1- 2 -',num2str(discount*100),'%']) 请输入商品价格：100折扣： 0%实际价格：100请输入商品价格：300 折扣： 3% 实际价格：291请输入商品价格：700 折扣： 5% 实际价格：665 请输入商品价格：1500 折扣： 8% 实际价格：1380请输入商品价格：3000折扣： 10% 实际价格：2700 请输入商品价格：6000 折扣： 14% 实际价格：5160- 3 -4．在.m 文件中编写程序实现，在同一坐标内，分别用不同线型和颜色绘制曲线0.510.2cos(4)y e x x π-=和0.522cos()y e x x π-=，标记两曲线交叉点，给出每条曲线的图注。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
z1=poly2str(z,'x'); w=polyval(z,x); plot(x,y,'o',x,w,'r')
实验体会与总结
-3-Βιβλιοθήκη 的值，8 3 2 8, b2 3, r2 15
1 10, b1 , r1 28
在三维空间上画出解轨线。 解：
equation.m equation_run.m
function dx=equation(t,x,a,b,c); dx=zeros(3,1); dx=[-a*x(1)+a*x(2);c*x(1)-x(2)-x(1)*x(3);x(1)*x(2)-b*x(3)];
实 验 报 告
班级 实 验 名 称 学号 姓名 Matlab 第一次试验 1．熟悉 Matlab 的操作环境和基本操作； 实 验 目 的 2. 掌握随机矩阵的构建； 3．掌握微分方程数值求解； 4. 掌握曲线拟合 参加实验时间 实 验 地 点 年 月 日 7、8 节课 电信楼 313
实验内容及其结果（写出源程序及其输出结果）
>> x0=[0,1e-10,1]; a1=10;b1=8/3;c1=28; [t,x]=ode45(@equation,[0,100],x0,[],a1,b1,c1); a2=8;b2=3;c2=15; [t,y]=ode45(@equation,[0,100],x0,[],a2,b2,c2); subplot(1,2,1);plot3(x(:,1),x(:,2),x(:,3)); subplot(1,2,2);plot3(y(:,1),y(:,2),y(:,3))
%
[i,j]= find(abs(A)>3);
0 1
2 3
4 5
-1-
x1 x1 x2 2、解方程 x2 rx1 x2 x1 x3 ,初值为[0,10-10,0]，t=[0,100]，方程参数设置成以下两组不同 x x x bx 1 2 3 3
4 2 0 2 4 1．找出数组 A 中所有绝对值大于 3 的元素。 3 1 1 3 5
编程如下： ex_find.m
A=[-4 -2 0 2 4 -3 -1 1 3 5] find(abs(A)>3); ans 输出结果如下： A= -4 -3 ans = 1 9 10 -2 -1
注：如果函数调用有附件参数，函数的调用请用@equation，而非‘equation’ ，前者是用函数 句柄，后者是通过函数名进行引用。若没有附加参数，两种方式都可 3、产生均值为 2，方差为 0.1 的 10x10 正态分布随机阵，然后统计其中大于 1.8，小于 2.1 的元素个数。
A=2+sqrt(0.1)*randn(10,10); B=A>1.8&A<2.1; sum(B(:))
4、求四次多项式满足下列数据点 X Y 0 1 1 7 2 23 3 109 4 307 5 400 6 307 10 123 12 69 15 24 16 6
-2-
输出 X,Y 满足的函数关系式，画出拟合曲线以及原数据点。 解： ex_polyfit.m x=[0,1,2,3,4,5,6,10,12,15,16]; y=[1,7,23,109,307,400,307,123,69,24,6]; z=polyfit(x,y,4);