《MATLAB语言与应用》实验课报告
学院:信息学院
班级:测控0902
姓名:陈白杨
学号:20092352
《MATLAB 语言与应用》实验课程任务书
一、 实验教学目标与基本要求
上机实验是本课程重要的实践教学环节;实验的目的不仅仅是验证理论知识,更重要的是通过上机实验,加强学生的实验手段与实践技能,掌握应用MATLAB 语言求解问题的方法,培养学生分析问题、解决问题、应用知识的能力和创新精神,全面提高学生的综合素质。
上机实验共8学时。主要实验内容是基于理论课所学知识对课后典型习题进行MATLAB 求解,基本掌握常见数学问题的求解方法与命令调用,更深入地认识和了解MATLAB 语言强大的计算功能。
上机实验最终以书面报告的形式提交,并作为期末成绩考核内容的一部分。
二、 实验内容(8学时)
第一部分MATLAB 语言编程、科学绘图与基本数学问题求解(4学时)
主要内容:掌握MATLAB 语言编程基础、科学绘图方法、微积分问题、线性代数问题等基本数学问题的求解与应用。 练习题:
1、 安装MATLAB 软件,应用demo 命令了解主要功能,熟悉基本功能,会用help 命令。
2、 用MATLAB 语句输入矩阵A 和B
?????????
???=14
23
143212344321
A , ??
???
??
??
???++++++++++++++++=4j 11j
43j 22j
34j 11j 42j 33j 24j 13j 22j 31j 41j 42j 33j 24j 1B 前面给出的是44?矩阵,如果给出5)6,5(=A 命令将得出什么结果?
>> A=[1 2 3 4;4 3 2 1;2 3 4 1;3 2 4 1]
A =
1 2 3 4 4 3 2 1 2 3 4 1 3 2 4 1
>> B=[1+4j,2+3j,3+2j,4+1j;4+1j,3+2j,2+3j,1+4j;2+3j,3+2j,4+1j,1+4j;3+2j,2+3j,4+1j,1+4j] 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
>> A(5,6)=5
A =
1 2 3 4 0 0
4 3 2 1 0 0
2 3 4 1 0 0
3 2
4 1 0 0
0 0 0 0 0 5
3、假设已知矩阵A,试给出相应的MATLAB命令,将其全部偶数行提取出来,赋给B矩
阵,用magic(8)
A=命令生成A矩阵,用上述命令检验一下结果是不是正确。
A=magic(8)
A =
64 2 3 61 60 6 7 57
9 55 54 12 13 51 50 16
17 47 46 20 21 43 42 24
40 26 27 37 36 30 31 33
32 34 35 29 28 38 39 25
41 23 22 44 45 19 18 48
49 15 14 52 53 11 10 56
8 58 59 5 4 62 63 1
>> B=A(2:2:end,:)
B =
9 55 54 12 13 51 50 16
40 26 27 37 36 30 31 33
41 23 22 44 45 19 18 48
8 58 59 5 4 62 63 1
4、用数值方法可以求出∑
=
+
+
+
+
+
+
=
=
63
63
622
2
8
4
2
1
2
i i
S ,试不采用循环的形式求出和
式的数值解。由于数值方法是采用double 形式进行计算的,难以保证有效位数字,所以结果不一定精确。试采用运算的方法求该和式的精确值。 >> format long ;sum(2.^[0:63])
ans =
1.844674407370955e+019
5、 选择合适的步距绘制出下面的图形。 (1))/1sin(t ,其中)1,1(-∈t ; t=[-1:0.02:1];
y=sin(1./t);
plot(x,y)
(2))tan(sin )sin(tan t t -,其中),(ππ-∈t 。 x=[-pi:0.05:pi];
y=sin(tan(x))-tan(sin(x)); plot(x,y)
6、试绘制出二元函数2
2
2
2
)1(1)1(1),(y
x y
x y x f z +++
+-=
=的三维图和三视图。
xx=[-2:0.1:-1.2,-1.1:0.02:-0.9,-.8:0.1:0.8,0.9:0.02:1.1,1.2:0.1:2]; yy=[-1:0.1:-0.2,-0.1:0.02:0.1,0.2:0.1:1]; [x,y]=meshgrid(xx,yy);
z=1./(sqrt((1-x).^2+y.^2))+1./(sqrt((1+x).^2+y.^2)); surf(x,y,z),shading flat;zlim([0,15])
7、试求出如下极限。 (1)x
x x x 1)93(lim +∞
→;
(2)1
1lim
0-+→→xy xy y x ; (3)2
2)()cos(1lim
2
2
220
0y x y x e
y x y x +→→++-。
8、已知参数方程?
??-==t t t y t x sin cos cos ln ,试求出x y d d 和3
/2
2d d π=t x y
。
syms t;y=cos(t)-sin(t)*t;x=log(cos(t));f1=diff(y,t)/diff(x,t)
f1 =(cos(t)*(2*sin(t) + t*cos(t)))/sin(t)
>> f=diff(y,t,2)/diff(x,t,2); subs(f,t,sym(pi)/3) ans =
3/8 - (pi*3^(1/2))/24
9、 假设?
-=xy
t t e y x f 0
d ),(2
,试求2222
22y
f
y x f x f y x ??+???-??。 10、 试求出下面的极限。
(1)??????-++-+-+-∞→1)2(1
161141121lim 2222n n ; syms n;s=symsum(1/((2*n)^2-1),n,1,inf) s = 1/2
(2))1
31211(
lim 2222π
πππn n n n n n n ++++++++∞→ 。 >> syms k n
limit(n*symsum(1/(n^2+k*pi),k,1,n),n,inf) ans = 1
11、 试求出以下的曲线积分。
(1)?+l
s y x d )(22,l 为曲线)sin (cos t t t a x +=,)cos (sin t t t a y -=,
)20(π≤≤t 。
(2)?-+++l
y y y xe x e yx )dy 2(xy d )(33,其中l 为22222c y b x a =+正向上半椭圆。 >>
syms t ;syms a b c positive;
x=c*cos(t)/a;y=c*sin(t)/b;f=[y*x^3+exp(y),x*y^3+x*exp(y)-2*y];ds=[diff(x,t);diff(y,t)];I=int(f
*ds,t,pi,0)
I =
(2*c*(15*b^4 - 2*c^4))/(15*a*b^4)
12、 (2*c*(15*b^4 - 2*c^4))/(15*a*b^4)试求出Vandermonde 矩阵???????
?
???????
?=1e e e e 1d d d d 1c c c c 1b b b b
1a a a a 2
3
4234234
234234A 的
行列式,并以最简的形式显示 结果。
syms a b c d e;A=[a^4 a^3 a^2 a 1;b^4 b^3 b^2 b 1;c^4 c^3 c^2 c 1; d^4 d^3 d^2 d 1;e^4 e^3 e^2 e 1];simple(det(A)) ans =
(a - b)*(a - c)*(a - d)*(b - c)*(a - e)*(b - d)*(b - e)*(c - d)*(c - e)*(d - e)
13、 试对矩阵????
??
???
???-------=22120.54.50.520.50.51.500.50.50.52A 进行Jordan 变换,并得出变换矩阵。 14、 试用数值方法和解析方法求取下面的Sylvester 方程,并验证得出的结果。
??
?
??
??
?
????????-------=??????????-----+?????
??
?
???
?????-------36644616521411291229212343040
01101013376364224150463X X
>>A=[3,-6,-4,0,5;1,4,2,-2,4;-6,3,-6,7,3;-13,10,0,-11,0;0,4,0,3,4];B=[3,-2,1;-2,-9,2;-2,-1,9];C=-[-2,1,-1;4,1,2;5,-6,1;6,-4,-4;-6,6,-3];X=lyap(A,B,C),norm(A*X+X*B+C) X =
4.0569 14.5128 -1.5653 -0.0356 -2
5.0743 2.7408 -9.4886 -25.9323 4.4177 -2.6969 -21.6450 2.8851 -7.7229 -31.9100 3.7634
ans =
4.0149e-013 >> A*X+X*B+C ans =
1.0e-012 *
0.0231 -0.0284 0.0071 -0.0924 -0.2274 0.0213 -0.0782 -0.1990 0.0213 -0.0497 -0.1137 0 -0.0995 -0.1705 0.0213
15、 假设已知矩阵A 如下,试求出At e ,At sin ,)sin(2t e A e At At 。
????
?
????
???---------=3110 1.52.511.50.50.540.5 1.50.50
4.5A A=[-4.5,0,0.5,-1.5;-0.5,-4,0.5,-0.5;1.5,1,-2.5,1.5;0,-1,-1,-3];B=expm(A) B =
0.0283 0.0283 0.0498 -0.0215 0.0034 0.0283 0.0249 0.0034 0.0464 0.0215 0.0747 0.0464 -0.0249 -0.0498 -0.0747 0.0249 > C=funm(A,'sin') C =
0.9745 -0.4400 -0.4244 1.1156 0.0550 0.4089 -0.4950 0.0550 -1.0450 -0.5500 -0.6361 -1.0450 -0.0706 0.9900 0.9194 -0.2117
>> D=expm(A)*funm((A^2)*expm(A),'sin') D =
0.0133 0.0181 0.0303 -0.0082 0.0039 0.0113 0.0141 0.0039 0.0244 0.0102 0.0357 0.0244 -0.0161 -0.0283 -0.0444 0.0054
第二部分 数学问题求解与数据处理(4学时)
主要内容:掌握代数方程与最优化问题、微分方程问题、数据处理问题的MATLAB 求解方法。 练习题:
1、对下列的函数)(t f 进行Laplace 变换。 (1)t
t
t f a αsin )(=
;(2)t t t f b αsin )(5=;(3)t t t f c αcos )(8=。 2、对下面的)(s F 式进行Laplace 反变换。 (1))
)((1)(222b s a s s s F a +-=
;(2)b s a s s F b ---=)(;(3)b
s a
s s F c --=ln
)(。 3、试求出下面函数的Fourier 变换,对得出的结果再进行Fourier 反变换,观察是否能得出
原来函数。
(1)ππ20),23()(2≤≤-=x x x x f ;(2)ππ20,)2()(22≤≤-=t t t t f 。 4、请将下述时域序列函数)(kT f 进行Z 变换,并对结果进行反变换检验。 (1))cos()(kaT kT f a =;(2)akT b e kT kT f -=2)()(;(3))1(1
)(akT c e akT a
kT f -+-=
。 5、用数值求解函数求解下述一元和二元方程的根,并对得出的结果进行检验。 (1))25sin(2
/)1()(2+++-=x x e
x f π;(2)xy
y x
e xy y x y x
f ---++=22
)(),(22。
6、试求出使得?-1
2d )(x cx e x 取得极小值的c 值。
>> syms x c; d=(exp(x)-c*x)^2;F=simple(int(d,x,0,1)) F =
c^2/3 - 2*c + exp(2)/2 - 1/2
>> syms x c;F=c^2/3 - 2*c + exp(2)/2 - 1/2;f=diff(F) f = (2*c)/3 - 2
>> syms c;c=solve('(2*c)/3 - 2') c = 3
7、试求解下面的非线性规划问题。
min )12424(2212
2211++++x x x x x e x
x ???????≤≤--≥≥++-≤+10
,10105.10
.s.t 21212
12121x x x x x x x x x x
function y= MyFun( x )
y=exp(x(1))*(4*x(1)^2+2*x(2)^2+4*x(1)*x(2)+2*x(2)+1); end
function [ c,ce ] = MyFun2( x ) ce=[];
c=[x(1)+x(2); x(1)*x(2)-x(1)-x(2)+1.5; -10-x(1)*x(2)];
end
>> A=[]; B=[]; Aeq=[]; Beq=[]; xm=[-10; -10]; xM=[10; 10]; x0=(xm+xM)/2;
ff=optimset; ff.TolX=1e-10; ff.TolFun=1e-20;i=1; x=x0; while (1)
[x,a,b]=fmincon('MyFun',x,A,B,Aeq,Beq,xm,xM,'MyFun2',ff); if b>0, break; end i=i+1; end x,i x =
1.1825 -1.7398 i = 5
8、求解下面的整数线性规划问题。
max )23374855273381592(7654321x x x x x x x ++++++
x ???≤++++++≥119567
3044515285891767235635340
.s.t 7654321x x x x x x x x
9、试求出微分方程x e x x y x
x y x x y 52)()1
1()()12()(-=-+-- 的解析解通解,并求出满足边界条
件1)(,)1(==ππy y 的解析解。
syms x; y=dsolve('D2y-(2-1/x)*Dy+(1-1/x)*y=x^2*exp(-5*x)','x') y =
C9*exp(x) + (30*x + 6*log(x) + 108*x^2*log(x) + 216*x^3*log(x) - 6*Ei(-6*x)*exp(6*x) + 36*x*log(x) + 36*x^2 + 11)/(1296*exp(5*x)) - (log(x)*(216*x^3 + 108*x^2 + 36*x + 6))/(1296*exp(5*x)) + C10*exp(x)*log(x)
syms x;y=dsolve('D2y-(2-1/x)*Dy+(1-1/x)*y=x^2*exp(-5*x)','y(1)=pi','y(pi)=1','x') y =
(30*x + 6*log(x) + 108*x^2*log(x) + 216*x^3*log(x) - 6*Ei(-6*x)*exp(6*x) + 36*x*log(x) + 36*x^2 + 11)/(1296*exp(5*x)) - (log(x)*(216*x^3 + 108*x^2 + 36*x + 6))/(1296*exp(5*x)) + (exp(x)*(1296*pi
-
77/exp(5)
+
6*Ei(-6)*exp(1)))/(1296*exp(1))
-
(exp(x)*log(x)*(1296*pi*exp(pi) - 1296*exp(1) - (77*exp(pi))/exp(5) + (11*exp(1))/exp(5*pi) + (30*pi*exp(1))/exp(5*pi) - 6*exp(1)*exp(pi)*Ei(-6*pi) + (36*pi^2*exp(1))/exp(5*pi) + 6*Ei(-6)*exp(1)*exp(pi)))/(1296*exp(1)*exp(pi)*log(pi)) 10、 试求出下面微分方程的通解。
(1)1)()(2)(2+=++t t x t t x t t x
;(2)2
)(2)(x xe x xy x y -=+ 。 11、 考虑著名的ssler o R 化学反应方程组???
??-+=+=--=z
c x b z
ay x y z y x
)( ,选定2.0==b a ,7.5=c ,且
)0()0()0(321x x x ==,绘制仿真结果的三维相轨迹,并得出其在x-y 平面上的投影。在实际求解中建议将c b a ,,作为附加参数,同样的方程若设2.0=a ,5.0=b ,10=c 时,绘制出状态变量的二维图和三维图。
>> f=inline('[-x(2)-x(3); x(1)+a*x(2); b+(x(1)-c)*x(3)]',... 't','x','flag','a','b','c');[t,x]=ode45(f,[0,100],[0;0;0],[],0.2,0.2,5.7); plot3(x(:,1),x(:,2),x(:,3));grid
>>plot(x(:,1),x(:,2));grid
>> f=inline('[-x(2)-x(3); x(1)+a*x(2); b+(x(1)-c)*x(3)]',... 't','x','flag','a','b','c');[t,x]=ode45(f,[0,100],[0;0;0],[],0.2,0.5,10); plot3(x(:,1),x(:,2),x(:,3)); grid
>> plot(x(:,1),x(:,2)); grid
12、 试选择状态变量,将下面的非线性微分方程组转换成一阶显式微分方程组,并用
MATLAB 对其求解,绘制出解的相平面或相空间曲线。
???
???
?==-===----=+++---=-6)1(,7)1(,2)1(4)1(,2)1(26)()3()
3(32y y y x
x t e x y y t y y x y x x x
>>f=inline(['[x(2);-x(1)-x(3)-(3*x(2))^2+(x(4))^3+6*x(5)+2*t;','x(4);x(5);-x(5)-x(2)-exp(-x(1))-t]'],'t','x');[t1,x1]=ode45(f,[1,0],[2, -4, -2, 7, 6]');[t2,x2]=ode45(f,[1,2],[2, -4, -2, 7, 6]'); t=[t1(end:-1:1); t2]; x=[x1(end:-1:1,:); x2]; plot(t,x)
figure;plot(x(:,1),x(:,3))
13、 考虑简单的线性微分方程)3/4sin(246553)3()4(π++=++++--t e e y y y
y y t t ,且方程的初值为1)0(=y ,2/1)0()0(==y y
,2.0)0()3(=y ,试用Simulink 搭建起系统的仿真模
型,并绘制出仿真结果曲线。
[t,x,y]=sim('Mymdl',[0,10]); plot(t,x)
>> [t,x,y]=sim('Mymdl',[0,10]); plot(t,y)
数学实验答案 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 = 1.2500
MATLAB数学实验练习题
————————————————————————————————作者:————————————————————————————————日期:
“MATLAB”练习题 要求:抄题、写出操作命令、运行结果,并根据要求,贴上运行图。 1、求230x e x -=的所有根。(先画图后求解)(要求贴图) >> solve('exp(x)-3*x^2',0) > ezplot('exp(x)-3*x^2') >> grid on 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 -.199936 -1.00450-1.06095*i 1.10447-1.05983*i 2)1sin 02 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,x,0) 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)
实验三 MATLAB 绘图 一、实验目的 1.掌握二维图形的绘制。 2.掌握图形的标注 3.了解三维曲线和曲面图形的绘制。 二、实验的设备及条件 计算机一台(带有MATLAB7.0以上的软件环境)。 设计提示 1.Matlab 允许在一个图形中画多条曲线:plot(x1,y1,x2,y2,……) 指令绘制y 1 = f 1(x 1), y 2 = f 2 (x 2 )等多条曲线。Matlab 自动给这些曲线以不同颜色。标注可用text 函数。 2.绘图时可以考虑极坐标和直角坐标的转换。 3.三维曲线绘图函数为plot3,注意参考帮助中的示例。 三、实验内容 1.生成1×10 维的随机数向量a ,分别用红、黄、蓝、绿色绘出其连线图、 杆图、阶梯图和条形图,并分别标出标题“连线图”、“杆图”、“阶梯图”、“条形图”。 2、绘制函数曲线,要求写出程序代码。 (1) 在区间[0:2π]均匀的取50个点,构成向量t (2) 在同一窗口绘制曲线y1=sin(2*t-0.3); y2=3cos(t+0.5);要求y1曲 线为红色点划线,标记点为圆圈;y2为蓝色虚线,标记点为星号。 (3) 分别在靠近相应的曲线处标注其函数表达式。 3.将图形窗口分成两个绘图区域,分别绘制出函数: ???+-=+=1 352221x x y x y 在[0,3]区间上的曲线,并利用axis 调整轴刻度纵坐标刻度,使1y 在[0,12] 区间上,2y 在[-2,1.5]区间上。 4.用mesh 或surf 函数,绘制下面方程所表示的三维空间曲面,x 和y 的
取值范围设为[-3,3]。 10102 2y x z +-= 思考题: 1. 编写一个mcircle(r)函数,调用该函数时,根据给定的半径r ,以原点 为圆心,画一个如图所示的红色空心圆。(图例半径r=5);左图参考 polar 函数的用法,右图绘制圆形的参数方程为x=sin (t ),y=cos (t )。其中,t 的区间为0~2*pi ,步长为0.1。 2.(1)绘一个圆柱螺旋线(形似弹簧)图。圆柱截面直径为10,高度为5, 每圈上升高度为1。如左图所示。 (2)利用(1)的结果,对程序做少许修改,得到如右图所示图形。
注意:在下面的题目中m 为你的学号的后3位(1-9班)或4位(10班以上). 第一次练习题 1.求解下列各题: 1)30sin lim x mx mx x ->- 2)(4)cos ,1000.0=x mx y e y 求 3)21/2 0mx e dx ?(求近似值,可以先用inline 定义被积函数,然后用quad 命令) 4)4 224x dx m x +? 5 0x =展开(最高次幂为8). 2.对矩阵21102041A m -?? ?= ? ?-?? ,分别求逆矩阵,特征值,特征向量,行列式,并求矩阵,P D (D 是对角矩阵),使得1A PDP -=。 3. 已知2 1(),()2f x e x μσ=--分别在下列条件下画出)(x f 的图形: (1)/600m σ=,μ分别为0,1,1-(在同一坐标系上作图); (2)0μ=,σ分别为1,2,4,/100m (在同一坐标系上作图). 4.画 (1)sin 020cos 02100x u t t y u t u t z m ??=≤≤?=?≤≤??=?
(2) sin()03,03z mxy x y =≤≤≤≤ (3)sin()(/100cos )02cos()(/100cos )02sin x t m u t y t m u u z u π π=+?≤≤?=+?≤≤?=? 的图(第4题只要写出程序). 5.对于方程50.10200 m x x --=,先画出左边的函数在合适的区间上的图形,借助于软件中的方程求根的命令求出所有的实根,找出函数的单调区间,结合高等数学的知识说明函数为什么在这些区间上是单调的,以及该方程确实只有你求出的这些实根。最后写出你做此题的体会. 第二次练习题 判断迭代收敛速度的程序 x0=1;stopc=1;eps=10^(-8);a=1;c=1;b=2*c;d=a;k=0; f=inline('(a*x+b)/(c*x+d)'); kmax=100; while stopc>eps&k result5 = ( 1 ) a = 1 : 2 : 5 a = 1 3 5 ( 2 ) b = [ a' , a' , a' ;a ] b = 1 1 1 3 3 3 5 5 5 1 3 5 ( 3 ) c = a + b ( 2 , : ) c = 4 6 8 2、下列运算是否合法,为什么如合法, 结果是多少 >> result2=a*b Error using * Inner matrix dimensions must agree. >> result3=a+b result3 = 3 6 2 58 11 >> result4=b*d result4 = 31 22 22 40 49 13 31 22 22 40 49 13 -5 -8 7 >> result6=a.*b result6 = 2 8 -3 415 30 >> result7=a./b result7 = >> result8= Attempt to reference field of non-structure array. >> result9=a.\b result9 = >> result10=a92 result10 = 1 4 9 16 25 36 >> resultl 1=29a result11 = 2 4 8 16 32 64 >> result5=[b;c']*d 3、用MATLAB求解下面的的方程组。 1、求以下变量的值,并在MATLAB^验证。 1 2 x1 3 2 x2 11 5 x3 2 1 3 x4 >> A=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> B=[4 7 -1 0] >> B=B' >> x=inv(A)*B >> A1=[1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6] >> B2=[1;8;3;5] >> x2=inv(A1)*B2 7 2 1 2 9 15 3 2 2 2 11 5 1 3 2 13 (1)求矩阵A的秩(rank) (2)求矩阵 A 的行列式(determinant) (3)求矩阵 A 的逆(inverse) (4)求矩阵 A 的特征值及特征向量 (eigenvalue and eigenvector) >> A3=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> r=rank(A3) >> b=inv(A3) >> a=det(A3) >> [V,D]=eig(A3) 10 n 10 查看y 的值) m1=0; for m=-10:10 m仁m1+2^m; end m1 m1 = 6、求分段函数的值。 用if 语句实现,算出下列表中x 对应的y 值。 x=input('enter x='); if x<0 y=x A2+x-6; elseif x>=0&&x<5 y=xA2-5*x+6; else y=xA2-x-1; end y 7、分别用if 和switch 语句实现,将百分 制成绩转换为成绩等级A、B、C、D、E。 其中90~1 00分为A,80~89 分为B,70~79 分为C,60~69 分为D,60 分以下为E。 对超出百分制范围的成绩,给出错误提示 信息。 if 结构程序: x=input('please enter score='); if x>=90&&x<=100 9 2 10 disp('A') 7 2 9 15 (1) 2 2 1 3 4 7 1 0 A 4、已知 2n 2 10 29 复习题 1、写出3 2、i nv(A)表示A的逆矩阵; 3、在命令窗口健入 clc,4、在命令窗口健入clea 5、在命令窗口健入6、x=-1:0.2:17、det(A)表示计算A的行列式的值;8、三种插值方法:拉格朗日多项式插值,分段线性插值,三次样条插值。 9、若A=123456789?? ????????,则fliplr (A)=321654987?????????? A-3=210123456--??????????A .^2=149162536496481?????????? tril(A)=100450789?????????? tri u(A,-1)=123456089??????????diag(A )=100050009?????????? A(:,2),=2 58A(3,:)=369 10、nor mcd f(1,1,2)=0.5%正态分布mu=1,s igm a=2,x =1处的概率 e45(@f,[a,b ],x0),中参数的涵义是@fun 是求解方程的函数M 文 件,[a,b ]是输入向量即自变量的范围a 为初值,x0为函数的初值,t 为输出指定的[a,b],x 为函数值 15、写出下列命令的功能:te xt (1,2,‘y=s in(x)’ hold on 16fun ction 开头; 17 ,4) 3,4) 21、设x 是一向量,则)的功能是作出将X十等分的直方图 22、interp 1([1,2,3],[3,4,5],2.5) Ans=4.5 23、建立一阶微分方程组? ??+='-='y x t y y x t x 34)(3)(2 的函数M 文件。(做不出来) 二、写出运行结果: 1、>>ey e(3,4)=1000 01000010 2、>>s ize([1,2,3])=1;3 3、设b=ro und (unifrnd(-5,5,1,4)),则=3 5 2 -5 >>[x,m]=min(b);x =-5;m=4 ,[x,n ]=sort(b ) -5 2 3 5 4 3 1 2 mea n(b)=1.25,m edian(b)=2.5,range(b)=10 4、向量b如上题,则 >>an y(b),all(b<2),all(b<6) Ans =1 0 1 5、>>[5 6;7 8]>[7 8;5 6]=00 11 6、若1234B ??=???? ,则 7、>>diag(d iag (B ))=10 04 8、>>[4:-2:1].*[-1,6]=-4 12 9、>>acos(0.5),a tan(1) ans = 1.6598 ans= 课程实验报告 学年学期2015-2016学年第2学期 课程名称MATLAB与科学计算 实验名称实验一 MATLAB环境与数值计算实验室测量测绘实验中心计算机室专业年级土木13级 学生姓名 学生学号 20130 提交时间 2016-06-02 成绩 任课教师王增红 水利与建筑工程学院 实验一 MATLAB环境与数值计算 1、目的和要求 (1)熟练掌握MATLAB的启动和退出、MATLAB的命令窗口;常用选单和工具栏,通过实例初步认识对MATLAB的命令和程序运行等。 (2)熟练掌握MATLAB变量的使用; (3)熟练掌握矩阵的创建、掌握MATLAB的矩阵和数组的运算; (4)熟悉MATLAB多项式的运用,掌握MATLAB的拟合和插值方法、积分运算和统计运算。2、内容和步骤 参见教材实验一、二。 3、实验报告提交要求 (1)用MALTAB命令和M程序两种方式求解三元一次方程组。 答:以求解三元一次方程组{2a?b+3c=5 3a+b?5c=5 4a?b+c=9 为例: (2)用MALTAB命令和M程序两种方式求解3×3复数矩阵的模和相角。答: (3)创建一个4╳4矩阵,求逆矩阵、产生对角矩阵、求矩阵的秩等。 (5)创建一个一元四次多项式,进行多项式求值,求根计算。 432 (6)创建一个一元三次多项式,取十个不同自变量值,计算多项式的值,并分别进行一阶、二阶拟合,图示原始数据和拟合结果 (7)自建一个一元三次多项式,取十个不同自变量值,计算多项式的值,并用线性、相邻点和三次插值,图示其插值结果。 解: (8)给出多组温度的测定值(内含平均温度、最低、最高温度),并求其平均温度、最低、最高温度的平均值、标准差等。 解: >> a=fix(40*rand(6,3)+1) 数学实验报告 姓名: 班级: 学号: 第一次实验任务 过程: a=1+3i; b=2-i; 结果: a+b =3.0000 + 2.0000i a-b =-1.0000 + 4.0000i a*b = 5.0000 + 5.0000i a/b = -0.2000 + 1.4000i 过程: x=-4.5*pi/180; y=7.6*pi/180; 结果: sin(abs(x)+y)/sqrt(cos(abs(x+y))) =0.2098 心得:对于matlab 中的角度计算应转为弧度。 (1)过程: x=0:0.01:2*pi; y1=sin(x); y2=cos(x); y3=exp(x); y4=log(x); plot(x,y1,x,y2,x,y3,x,y4) plot(x,y1,x,y2,x,y3,x,y4) 结果: (2)过程:>> subplot(2,2,1) >> plot(x,y1) >> subplot(2,2,2) >> plot(x,y2) ./,,,,2,311b a b a b a b a i b i a ?-+-=+=计算、设有两个复数 6,7,5.4) cos()sin(2=-=++y x y x y x ,其中、计算的图形。 下分别绘制)同一页面四个坐标系)同一坐标系下(、在( x y e y x y x y x ln ,,cos ,sin 213==== >> subplot(2,2,3) >> plot(x,y3) >> subplot(2.2.4) >> subplot(2,2,4) >> plot(x,y4) 结果: 心得:在matlab中,用subplot能够实现在同一页面输出多个坐标系的图像,应注意将它与hold on进行区别,后者为在同一坐标系中划出多条曲线。 5、随机生成一个3x3矩阵A及3x2矩阵B,计算(1)AB,(2)对B中每个元素平方后得到的矩阵C,(3)sinB,(4)A的行列式,(5)判断A是否可逆,若可逆,计算A的逆矩阵,(6)解矩阵方程AX=B,(7)矩阵A中第二行元素加1,其余元素不变,得到矩阵D,计算D。 过程:A=fix(rand(3,3).*10) ; B=fix(rand(3,3).*10); 1、求以下变量的值,并在MATLAB 中验证。 ( 1 ) a = 1 : 2 : 5 a = 1 3 5 ( 2 ) b = [ a' , a' , a' ;a ] b = 1 1 1 3 3 3 5 5 5 1 3 5 ( 3 ) c = a + b ( 2 , : ) c = 4 6 8 2、下列运算是否合法,为什么?如合法,结果 是多少? >> result2=a*b Error using * Inner matrix dimensions must agree. >> result3=a+b result3 = 3 6 2 5 8 11 >> result4=b*d result4 = 31 22 22 40 49 13 >> result5=[b;c']*d result5 = 31 22 22 40 49 13 -5 -8 7 >> result6=a.*b result6 = 2 8 -3 4 1 5 30 >> result7=a./b result7 = 0.5000 0.5000 -3.0000 4.0000 1.6667 1.2000 >> result8=a.c Attempt to reference field of non-structure array. >> result9=a.\b result9 = 2.0000 2.0000 -0.3333 0.2500 0.6000 0.8333 >> result10=a.^2 result10 = 1 4 9 16 25 36 >> result11=2.^a result11 = 2 4 8 16 32 64 3、用MATLAB 求解下面的的方程组。 (1)????????????-=?????????????????????????----01741323151122231592127 4321x x x x >> A=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> B=[4 7 -1 0] >> B=B' >> x=inv(A)*B (2)???????=-++=--=-++=++56533332821 w z y x w y x w z y x z y x >> A1=[1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6] >> B2=[1;8;3;5] >> x2=inv(A1)*B2 4、已知????????????----=1323151122231592127A (1)求矩阵A 的秩(rank) 实验十周期函数 【实验目的】 1.了解周期函数的基本概念。 2.了解周期函数经过四则运算、复合运算、求导运算、积分运算后的周期性。 3.学习、掌握MATLAB软件有关命令。 【实验内容】 从图形上观察六个三角函数的周期性。 【实验准备】 1.周期函数的基本概念 函数() f x是以T为周期的周期函数是指对任何x,有 += ()() f x T f x 使得上式成立的最小正数T称为函数的最小正周期。 2.周期函数的四则运算 3.周期函数的最小正周期 【实验重点】 1、周期函数的四则运算与复合 2、周期函数的求导与积分运算 【实验难点】 1、最小正周期的确定 【实验方法与步骤】 练习 1 图形上观察六个三角函数 sinx,cosx,tanx,cotx,secx,cscx的周期性。绘制正弦函数y=sinx 在区间[-6π,6π]的图形,相应的MATLAB代码为 >>x=-6*pi:2*pi/30:6*pi; >>y=sin(x); >>plot(x,y); >>xlabel('x');ylabel('y'); 运行结果见图10.1。 从图形中可以看出y=sinx为周期函数,最小正周期T≈6。实际上,最小正周期T=2π=6.28…。同样,可以画出余弦函数y=cosx的图形,见图10.2,其最小正周期也为T=2π。 画正切函数y=tanx 的图形时,要注意函数在,0,1,2,...2x k k π π=+=±±是不连续,所以我们只能分别绘出函数在区间(,),0,1,2,...22k k k ππ ππ-+=±±的图形。相应的MATLAB 代码为 >>x=-1.5:0.01:1.5; >>x1=x-pi;x2=x+pi; >>y=tan(x);y1=tan(x1);y2=tan(x2); >>plot(x,y,x1,y1,x2,y2); >>xlabel('x');ylabel('y'); 运行结果见图10.3。 第一次练习 教学要求:熟练掌握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 lim x mx mx x →-与3 sin lim x mx mx x →∞- syms x limit((902*x-sin(902*x))/x^3) ans = 366935404/3 limit((902*x-sin(902*x))/x^3,inf) ans = 0 1.2 cos 1000 x mx y e =,求''y syms x diff(exp(x)*cos(902*x/1000),2) ans = (46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算 22 11 00 x y e dxdy +?? dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1) ans = 2.1394 1.4 计算4 2 2 4x dx m x +? syms x int(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 x diff(exp(x)*cos(902*x),10) ans = -356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x) 1.6 0x =的泰勒展式(最高次幂为4). 数学实验—实验报告(免积分) 一、实验项目:Matlab实验三—迭代 二、实验目的和要求 a.熟悉MATLAB软件的用户环境,掌握其一般目的命令和MATLAB数组操作与 运算函数; b.掌握MATLAB软件的绘图命令,能够熟练应用循环和选择结构实现各种循环 选择功能; c.借助MATLAB软件的绘图功能,对函数的特性进行探讨,广泛联想,大胆猜 想,发现进而证实其中的规律。 三、实验内容 问题一:将方程53 x x x +-+=改写成各种等价的形式进行迭代 5210 观察迭代是否收敛,并给出解释。 问题二:迭代以下函数,分析其收敛性。 4 f(x)=x-a 使用线性连接图、蛛网图或费根鲍姆图对参数a进行讨论和观察,会得到什么结论? 问题一: (1)画图 x1=-6:0.01:6; x2=-3:0.01:3; x3=-1:0.01:1; x4=-0.8:0.01:-0.75; y1=x1.^5 +5*x1.^3-2*x1+1; y2=x2.^5 +5*x2.^3-2*x2+1; y3=x3.^5 +5*x3.^3-2*x3+1; y4=x4.^5 +5*x4.^3-2*x4+1; subplot(2,2,1),plot(x1,y1) ,title('图(1)') ,grid on, subplot(2,2,2),plot(x2,y2) ,title('图(2)'),grid on, subplot(2,2,3),plot(x3,y3) ,title('图(3)'),grid on, subplot(2,2,4),plot(x4,y4) ,title('图(4)') ,grid on, 第1章MATLAB基础 1.4 MATLAB操作桌面有几个窗口?如何使某个窗口脱离桌面成为独立窗口?又如何将脱离出去的窗口重新放置到桌面上? 答:在MATLAB操作桌面上有五个窗口。在每个窗口的右上角有两个小按钮,一个是关闭窗口的Close 按钮,一个是可以使窗口成为独立窗口的Undock按钮,点击Undock按钮就可以使该窗口脱离桌面成为独立窗口。在独立窗口的Desktop菜单中选择Dock...项就可以将独立的窗口重新放置到桌面上。 1.5 如何启动M文件编辑/调试器? 答:在操作桌面上选择“建立新文件”或“打开文件”操作时,M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗?如何操作? 答:存储在工作空间的数组可以通过数组编辑器进行编辑:在工作空间浏览器中双击要编辑的数组名打开数组编辑器,再选中要修改的数据单元,输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外,还有什么用途? 答:命令历史窗口除了用于查询以前键入的命令外,还可以直接执行命令历史窗口中选定的内容、将选 定的内容拷贝到剪贴板中、将选定内容直接拷贝到M 文件中。 1.8 如何设置当前目录和搜索路径,在当前目录上的文件和在搜索路径上的文件有什么区别? 答:当前目录可以在当前目录浏览器窗口左上方的输入栏中设置,搜索路径可以通过选择操作桌面的file 菜单中的Set Path 菜单项来完成。在没有特别说明的情况下,只有当前目录和搜索路径上的函数和文件能够被MATLAB 运行和调用,如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件,如果没有特别说明,数据文件将存储在当前目录上。 1.9 在MATLAB 中有几种获得帮助的途径? 答:(1)帮助浏览器:选择view 菜单中的Help 菜单项或选择Help 菜单中的MATLAB Help 菜单项可以打开帮助浏览器。 (2)help 命令:在命令窗口键入“help ” 命令可以列出帮助主题,键入“help 函数名”可以得到指定函数的在线帮助信息。 (3)lookfor 命令:在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数。 (4)模糊查询:输入命令的前几个字母,然后按Tab 键,就可以列出所有以这几个字母开始的命令和函数。 注意:lookfor 和模糊查询查到的不是详细信息,通常还需要在确定了具体函数名称后用help 命令显示详细信息。 第2章 MATLAB 矩阵运算基础 2.1 在MATLAB 中如何建立矩阵?? ? ? ??194375,并将其赋予变量a ? 答:在Command Window 窗口输入操作: 实验一Matlab使用方法和程序设计 一、实验目的 1、掌握Matlab软件使用的基本方法; 2、熟悉Matlab的数据表示、基本运算和程序控制语句 3、熟悉Matlab绘图命令及基本绘图控制 4、熟悉Matlab程序设计的基本方法 二、实验内容: 1、帮助命令 使用help命令,查找sqrt(开方)函数的使用方法; 解:sqrt Square root Syntax B = sqrt(X) Description B = 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. Remarks See sqrtm for the matrix square root. Examples sqrt((-2:2)') ans = 0 + 1.4142i 0 + 1.0000i 1.0000 1.4142 2、矩阵运算 (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)矩阵的转置及共轭转置 1、求以下变量的值,并在MATLAB中验证。( 1 ) a = 1 : 2 : 5 a = 1 3 5 ( 2 ) b = [ a' , a' , a' ;a ] b = 1 1 1 3 3 3 5 5 5 1 3 5 ( 3 ) c = a + b ( 2 , : ) c = 4 6 8 2、下列运算是否合法,为什么?如合法, 结果是多少? >> result2=a*b Error using * Inner matrix dimensions must agree. >> result3=a+b result3 = 3 6 2 5 8 11 >> result4=b*d result4 = 31 22 22 40 49 13 >> result5=[b;c']*d result5 = 31 22 22 40 49 13 -5 -8 7 >> result6=a.*b result6 = 2 8 -3 4 1 5 30 >> result7=a./b result7 = 0.5000 0.5000 -3.0000 4.0000 1.6667 1.2000 >> result8=a.c Attempt to reference field of non-structure array. >> result9=a.\b result9 = 2.0000 2.0000 -0.3333 0.2500 0.6000 0.8333 >> result10=a.^2 result10 = 1 4 9 16 25 36 >> result11=2.^a result11 = 2 4 8 16 32 64 3、用MATLAB求解下面的的方程组。 (1) ? ? ? ? ? ? ? ? ? ? ? ? - = ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? - - - - 1 7 4 13 2 3 1 5 11 2 2 2 3 15 9 2 1 2 7 4 3 2 1 x x x x >> A=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> B=[4 7 -1 0] >> B=B' >> x=inv(A)*B (2) ? ? ? ? ? ? ? = - + + = - - = - + + = + + 5 6 5 3 3 3 3 2 8 2 1 w z y x w y x w z y x z y x >> A1=[1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6] >> B2=[1;8;3;5] >> x2=inv(A1)*B2 4、已知 ? ? ? ? ? ? ? ? ? ? ? ? - - - - = 13 2 3 1 5 11 2 2 2 3 15 9 2 1 2 7 A MATLAB 数学实验报告 指导老师: 班级: 小组成员: 时间:201_/_/_ Matlab 第二次实验报告 小组成员: 1 题目:实验四,MATLAB 选择结构与应用实验 目的:掌握if 选择结构与程序流程控制,重点掌握break,return , pause语句的应用。 问题:问题1:验证“哥德巴赫猜想” ,即:任何一个正偶数(n>=6)均可表示为两个质数的和。要求编制一个函数程序,输入一个正偶数,返回两个质数的和。 问题分析:由用户输入一个大于6 的偶数,由input 语句实现。由if 判断语句判断是否输入的数据符合条件。再引用质数判断函数来找出两个质数,再向屏幕输出两个质数即可。 编程:function [z1,z2]=gede(n); n=input('please input n') if n<6 disp('data error'); return end if mod(n,2)==0 for i=2:n/2 k=0; for j=2:sqrt(i) if mod(i,j)==0 k=k+1; end end for j=2:sqrt(n-i) if mod(n-i,j)==0 k=k+1; end end if k==0 fprintf('two numbers are') fprintf('%.0f,%.0f',i,n- i) break end end end 结果分析 如上图,用户输入了大于6的偶数返回两个质数5和31,通过 不断试验,即可验证哥德巴赫猜想。 纪录:if判断语句与for循环语句联合嵌套使用可使程序结构更加明晰,更快的解决问题。 2题目:实验四,MATLAB选择结构与应用实验 目的:用matlab联系生活实际,解决一些生活中常见的实际问 题。 Matlab 数学实验报告 一、实验目的 通过以下四组实验,熟悉MATLAB的编程技巧,学会运用MATLAB的一些主要功能、命令,通过建立数学模型解决理论或实际问题。了解诸如分岔、混沌等概念、学会建立Malthu模型和Logistic 模型、懂得最小二乘法、线性规划等基本思想。 二、实验内容 2.1实验题目一 2.1.1实验问题 Feigenbaum曾对超越函数y=λsin(πx)(λ为非负实数)进行了分岔与混沌的研究,试进行迭代格式x k+1=λsin(πx k),做出相应的Feigenbaum图 2.1.2程序设计 clear;clf; axis([0,4,0,4]); hold on for r=0:0.3:3.9 x=[0.1]; for i=2:150 x(i)=r*sin(3.14*x(i-1)); end pause(0.5) for i=101:150 plot(r,x(i),'k.'); end text(r-0.1,max(x(101:150))+0.05,['\it{r}=',num2str(r)]) end 加密迭代后 clear;clf; axis([0,4,0,4]); hold on for r=0:0.005:3.9 x=[0.1]; for i=2:150 x(i)=r*sin(3.14*x(i-1)); end pause(0.1) for i=101:150 plot(r,x(i),'k.'); end end 运行后得到Feigenbaum图 2.2实验题目二 2.2.1实验问题 某农夫有一个半径10米的圆形牛栏,长满了草。他要将一头牛拴在牛栏边界的桩栏上,但只让牛吃到一半草,问拴牛鼻子的绳子应为多长? 2.2.2问题分析 如图所示,E为圆ABD的圆心,AB为拴牛的绳子,圆ABD为草场,区域ABCD为牛能到达的区域。问题要求区域ABCD等于圆ABC 的一半,可以设BC等于x,只要求出∠a和∠b就能求出所求面积。先计算扇形ABCD的面积,2a÷π×πx2=2aπ2,再求AB的面积,用扇形ABE的面积减去三角形ABE的面积即可。 matlab练习题和答案 控制系统仿真实验 Matlab部分实验结果 目录 实验一 MATLAB基本操 作 ............................................................................................ 1 实验二 Matlab编程 .................................................................................................... 5 实验三Matlab底层图形控制 (6) 实验四控制系统古典分析.............................................................................................12 实验五控制系统现代分析 . (15) 实验六 PID控制器的设 计 ...........................................................................................19 实验七系统状态空间设计.............................................................................................23 实验九直流双闭环调速系统仿真 . (25) 实验一 MATLAB基本操作 1 用MATLAB可以识别的格式输入下面两个矩阵 1233,,1443678,i,,,,2357,,,,2335542,i,,,, A,1357B,,,2675342, i,,3239,,,,189543,,,,1894,, 再求出它们的乘积矩阵C,并将C矩阵的右下角2×3子矩阵赋给D矩阵。赋值完成后,调 用相应的命令查看MATLAB工作空间的占用情况。 A=[1,2,3,3;2,3,5,7;1,3,5,7;3,2,3,9;1,8,9,4]; 实验一 MATLAB 运算基础 1. 先求下列表达式的值,然后显示MATLAB 工作空间的使用情况并保存全部变量。 (1) 0 12 2sin851z e =+ (2) 221ln(1)2 z x x =++,其中2 120.45 5i x +??=? ?-?? (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+= ++=-- (4) 2242011 122123t t z t t t t t ?≤ =-≤?-+≤ ,其中t =0:0.5:2.5 解: 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 =MATLAB实验题答案
matlab数学实验复习题(有标准答案)
matlab实验(西北农林科技大学)
浅析Matlab数学实验报告
MATLAB实验题答案
matlab实验十--周期函数
南邮MATLAB数学实验答案(全)
matlab 数学实验 迭代 _ 蛛网图(免积分)
matlab基础练习题及答案讲解
Matlab实验第一次实验答案
MATLAB实验题答案
matlab实验报告
MATLAB数学实验报告
matlab练习题和答案
MATLAB)课后实验答案[1]