12 【例132-2】简单矩阵A=456的输入步骤。 在键盘上输入下列内容 [1,23;4,5,6:7,8,9 (2)按【 Enter】键,指令被执行。 (3)在指令执行后, MATLAB指令窗中将显示以下结果: 369 【例13.2-3】矩阵的分行输入。 A=[1,2,3 4,5,6 7,8,9 【例1.3.2-4】指令的续行输入 s=1-1/2+1/3-1/4+ 1/5-1/6+1/7-1/8 0.6345 1.3.3 数值、变量和表达式 数值的记述 二变量命名规则 三 MATLAB默认的预定义变量 四运算符和表达式 五复数和复数矩阵 【例133-1】复数z=3+4,2=1+2,3=2e6表达,及计算x=~二y (1) z1=3+4 Z1 (2)
【例 1.3.2-2】简单矩阵 ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎣ ⎡ = 7 8 9 4 5 6 1 2 3 A 的输入步骤。 （1）在键盘上输入下列内容 A = [1,2,3; 4,5,6; 7,8,9] （2）按【Enter】键，指令被执行。 （3）在指令执行后，MATLAB 指令窗中将显示以下结果： A = 1 2 3 4 5 6 7 8 9 【例 1.3.2-3】矩阵的分行输入。 A=[1,2,3 4,5,6 7,8,9] A = 1 2 3 4 5 6 7 8 9 【例 1.3.2-4】指令的续行输入 S=1-1/2+1/3-1/4+ ... 1/5-1/6+1/7-1/8 S = 0.6345 1.3.3 数值、变量和表达式 一 数值的记述 二 变量命名规则 三 MATLAB 默认的预定义变量 四 运算符和表达式 五 复数和复数矩阵 【例 1.3.3-1】复数 i z i z i z e 6 1 3 4 , 2 1 2 , 3 2 π = + = + = 表达，及计算 3 1 2 z z z z = 。 （1） z1= 3 + 4i z1 = 3.0000 + 4.0000i （2）
1+2★⊥ =2*exp(主*pi/6) z=z1*z2/z3 1.7321+1.00001 0.3349+5.5801i 【例13.3-2】复数矩阵的生成及运算 A=[1,3;2,4]-[5,8:6,9]* B=[1+51,2+61;3+8*1,4+9*11 C=A*B 1.0000-5.0000i3.0000-8.0000i 2.0000-6.0000i4.0000-9.0000i 1.0000+5.0000i2.0000+6.0000i 3.0000+8.0000i4.0000+9.0000i 1.0e+002 1.1600-0.0900i 【例133-3】求上例复数矩阵C的实部、虚部、模和相角。 c real=real(c) c imag=imag(C) C magnitude=abs(C) c phase=angle(C)*180/pi 以度为单位计算相角 C real 99116 116137 imag C magnitude 99.0000116.3486 116.3486137.0000 c phase 0-4.4365 【例1334】用 MATLAB计算√-8能得到-2吗? r=a^(1/3) 1.0000+1.7321 R=abs(a)^(1/3); Theta=(angle(a)+2*pi*m)/3: rrr=R*exp (i*Theta) 1= 1.0000+1.7321i-2.0000+0.0000i1.0000-1.7321i (3)
z2 = 1 + 2 * i z3=2*exp(i*pi/6) z=z1*z2/z3 z2 = 1.0000 + 2.0000i z3 = 1.7321 + 1.0000i z = 0.3349 + 5.5801i 【例 1.3.3-2】复数矩阵的生成及运算 A=[1,3;2,4]-[5,8;6,9]*i B=[1+5i,2+6i;3+8*i,4+9*i] C=A*B A = 1.0000 - 5.0000i 3.0000 - 8.0000i 2.0000 - 6.0000i 4.0000 - 9.0000i B = 1.0000 + 5.0000i 2.0000 + 6.0000i 3.0000 + 8.0000i 4.0000 + 9.0000i C = 1.0e+002 * 0.9900 1.1600 - 0.0900i 1.1600 + 0.0900i 1.3700 【例 1.3.3-3】求上例复数矩阵 C 的实部、虚部、模和相角。 C_real=real(C) C_imag=imag(C) C_magnitude=abs(C) C_phase=angle(C)*180/pi %以度为单位计算相角 C_real = 99 116 116 137 C_imag = 0 -9 9 0 C_magnitude = 99.0000 116.3486 116.3486 137.0000 C_phase = 0 -4.4365 4.4365 0 【例 1.3.3-4】 用 MATLAB 计算3 − 8 能得到 –2 吗？ （1） a=-8; r=a^(1/3) r = 1.0000 + 1.7321i （2） m=[0,1,2]; R=abs(a)^(1/3); Theta=(angle(a)+2*pi*m)/3; rrr=R*exp(i*Theta) rrr = 1.0000 + 1.7321i -2.0000 + 0.0000i 1.0000 - 1.7321i （3）
t=0: pi/20: 2*pi;x=R*sin(t)iy=R*cos(t); plot(x, y, 'b: ' ),grid hold or plot(rrr(1),' Markersize,50, 'Color!,'r) plot(rrr([2, 3]),'o','MarkerSize,15,'Color','b) axis([-3,3,-3,3]),axis square hold off 图1.3-2 1.3.4 计算结果的图形表示 【例1.341】画出衰减振荡曲线y=e3sin3及其它的包络线y0=e3。t的取值范围是 [0.4x]。(图 t=0:pi/50:4*pi; yO=exp(-t/3) exp(-t/3). *sin(3*t t,y0,':b',t, d 图1.3-3 【例1342】画出二=m√x2+y 所表示的三维曲面(图1.3-4)。x,y的取值范围是 c1ear;x=-8:0.5:8 ★ Y=Y*ones(size(x)) R=sqrt(X.^2+Y.^2)+eps;8<5> Z=sin(R)/R 8<6 urf(x, Y, z)
t=0:pi/20:2*pi;x=R*sin(t);y=R*cos(t); plot(x,y,'b:'),grid hold on plot(rrr(1),'.','MarkerSize',50,'Color','r') plot(rrr([2,3]),'o','MarkerSize',15,'Color','b') axis([-3,3,-3,3]),axis square hold off -3 -2 -1 0 1 2 3 -3 -2 -1 0 1 2 3 图 1.3-2 1.3.4 计算结果的图形表示 【例 1.3.4-1】画出衰减振荡曲线 y e t t sin 3 3 − = 及其它的包络线 3 0 t y e − = 。t 的取值范围是 [0,4π ]。（图 1.3-3） t=0:pi/50:4*pi; y0=exp(-t/3); y=exp(-t/3).*sin(3*t); plot(t,y,'-r',t,y0,':b',t,-y0,':b') grid 0 2 4 6 8 10 12 14 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 图 1.3-3 【例 1.3.4-2】画出 2 2 2 2 sin( ) x y x y z + + = 所表示的三维曲面（图 1.3-4）。 x, y 的取值范围是 [−8,8]。 clear;x=-8:0.5:8; y=x'; X=ones(size(y))*x; Y=y*ones(size(x)); R=sqrt(X.^2+Y.^2)+eps; %<5> Z=sin(R)./R; %<6> surf(X,Y,Z); %
co1。rmap(coo1) xlabel('x'),ylabel('y') zlabel('z') 图1.3-4 14 Command window操作要旨 1.4.1 指令窗显示方式的操作 缺省显示方式 二显示方式的设置 14.2 数值计算结果的显示格式 1.4.3 指令行中的标点符号 1.4.4 指令窗的常用控制指令 14.5 指令窗中指令行的编辑 【例145-1】指令行操作过程示例 15 Command history和实录指令dary 1.5.1 Command history历史指令窗简介
colormap(cool) % xlabel('x'),ylabel('y'),zlabel('z') 图 1.3-4 1.4Command Window 操作要旨 1.4.1 指令窗显示方式的操作 一 缺省显示方式 二 显示方式的设置 1.4.2 数值计算结果的显示格式 1.4.3 指令行中的标点符号 1.4.4 指令窗的常用控制指令 1.4.5 指令窗中指令行的编辑 【例 1.4.5-1】指令行操作过程示例。 1.5Command History 和实录指令 diary 1.5.1 Command History 历史指令窗简介
Command History 回区 File Edit iew Web Window Help A(:)=-4:5 ‰运用“全元素”法向A赋值 0=2*(rand0s,M).5)-1 disp(i), di sp gj〕 S0=2*rand0s,M).5)-1 图1.5-1 1.52 历史指令行的再运行 【例1.52-1】再运行图1.5-2所示历史指令窗中的三行指令 V MATLAB File Edit View Web Window Help D哆易啪舞|? Current Directory: D: MATLAB6pswv圆 Comma nd Hitory 7 a=0.9e-323 ‖A= zeros②2,5) %预页生成一个(2*5)全零数△ 10/22/023:38PA()=-45 运用“全元素”法向A值 A=zeros(2,5 Labs(A)3%产生 logical l) 5 Astart 图1.5-2 1.5.3 指令窗实录指令dary 16 Current directory、路径设置器和文件管理 1.6.1 Current Directory当前目录浏览器简介
图 1.5-1 1.5.2 历史指令行的再运行 【例 1.5.2-1】再运行图 1.5-2 所示历史指令窗中的三行指令。 图 1.5-2 1.5.3 指令窗实录指令 diary 1.6Current Directory、路径设置器和文件管理 1.6.1 Current Directory 当前目录浏览器简介