第6章MATLAB计算结果可视化讲解

第六章MATLAB 计算结果可视化

6.1连续函数和离散函数的可视化

【例6-1】用图形表示离散函数1

)

6(--=n y 。

n=0:12; %产生一组自变量数据 y=1./abs(n-6); %计算相应点的函数值 plot(n,y,'r*','MarkerSize',20) %用红花标出数据点 grid on %画坐标方格

【例6-2】用图形表示连续调制波形)9sin()sin(t t y =。

t1=(0:11)/11*pi; y1=sin(t1).*sin(9*t1);

t2=(0:100)/100*pi; y2=sin(t2).*sin(9*t2);

subplot(2,2,1),plot(t1,y1,'r.'),axis([0,pi,-1,1]),title('子图 (1)') subplot(2,2,2),plot(t2,y2,'r.'),axis([0,pi,-1,1]),title('子图 (2)') subplot(2,2,3),plot(t1,y1,t1,y1,'r.') axis([0,pi,-1,1]),title('子图 (3)') subplot(2,2,4),plot(t2,y2)

6.2二维曲线绘图的基本操作

6.2.1 plot 的基本调用格式

【例6-3】用图形表示连续调制波形)9sin()sin(t t y 及其包络线。

t=(0:pi/100:pi)'; %长度为101的时间采样列向量 y1=sin(t)*[1,-1]; %包络线函数值，是（101x2）的矩阵 y2=sin(t).*sin(9*t); %长度为101的调制波列向量 t3=pi*(0:9)/9; y3=sin(t3).*sin(9*t3);plot(t,y1,'r:',t,y2,'b',t3,y3,'bo')

6.2.2泛函绘图指令fplot

【例6-4】fplot 与一般绘图指令的绘图效果比较。

[x,y]=fplot('cos(tan(pi*x))',[-0.4,1.4],0.2e-3);n=length(x); subplot(1,2,1),plot(x,y)

title('\fontsize{20}\fontname{隶书}泛函绘图指令效果') t=(-0.4:1.8/n:1.4)';

subplot(1,2,2),plot(t,cos(tan(pi*t)))

6.2.3曲线的色彩、线型和数据点形

【例6-5】用图形演示平面上一个方块四个顶点在仿射投影（Affine Projection ）下的位置、

形状变化。

%平面上的四个点和它们构成的方块

p1=[-0.5,0,1]';p2=[-0.5,1,1]';p3=[0.5,1,1]';p4=[0.5,0,1]';

Sq=[p1,p2,p3,p4,p1];

%平移投影：沿x轴移动0.5 , 沿y轴移动1 。

dx=0.5;dy=1;T=[1,0,dx;0,1,dy;0,0,1];

%旋转投影：逆时针旋转30度。

th=pi/6;R=[cos(th),-sin(th),0;sin(th),cos(th),0;0,0,1];

%刻度投影：x方向放大到2倍，y方向放大到3倍。

alpha=2;beta=3;S=[alpha,0,0;0,beta,0;0,0,1];

E=eye(3,3);%为编程方便，设计一个单位阵。在它投影下，任何形状都不变。

TRS={E,T,R,S}; %用元胞数组存放四个变换矩阵

ss={'r^','rd','rp','rh'}; %用元胞数组存放数据点形设置值

tt={'Original Square','Translation','Rotation','Scaling'};

%用元胞数组存放四张子图的图名

for i=1:4

W=TRS{i}*Sq; %进行仿射变换

subplot(2,2,i)

for k=1:4

plot(W(1,k),W(2,k),ss{k}); %用不同点形标志四个顶点

axis([-3,3,-1,5]),axis equal

hold on %使以后图形画在当前子图上end

plot(W(1,:),W(2,:)) %连接四个顶点使图形封闭

grid on %画过坐标格线

title(tt{i}) %给子图题写图名

hold off %使以后图形不再画在当前子图上

end

6.2.4坐标控制

【例6-6】观察各种轴控制指令的影响。演示采用长轴为3.25，短轴为1.15的椭圆。注意：采用多子图表现时，图形形状不仅受“控制指令”影响，而且受整个图面“宽高比”及“子图数目”的影响。

t=0:2*pi/99:2*pi;

x=1.15*cos(t);y=3.25*sin(t); %y为长轴，x为短轴

subplot(2,3,1),plot(x,y),axis normal,grid on,

title('Normal and Grid on')

subplot(2,3,2),plot(x,y),axis equal,grid on,title('Equal')

subplot(2,3,3),plot(x,y),axis square,grid on,title('Square')

subplot(2,3,4),plot(x,y),axis image,box off,title('Image and Box off') subplot(2,3,5),plot(x,y),axis image fill,box off

title('Image and Fill')

6.2.5刻度、分格线和坐标框

【例6-7】绘制二阶系统阶跃响应。

clf;t=6*pi*(0:100)/100;y=1-exp(-0.3*t).*cos(0.7*t);

tt=t(find(abs(y-1)>0.05));ts=max(tt);

subplot(1,2,1),plot(t,y,'r-','LineWidth',3),grid on

axis([0,6*pi,0.6,max(y)])

title('y=1–exp(-alpha*t)*cos(omega*t)')

text(11,1.25,'alpha=0.3');text(11,1.15,'omega=0.7')

hold on;plot(ts,0.95,'bo','MarkerSize',10);hold off

text(ts+1.5,0.95,['ts=' num2str(ts)])

xlabel('t -->'),ylabel('y -->')

subplot(1,2,2),plot(t,y,'r-','LineWidth',3)

axis([-inf,6*pi,0.6,inf])

set(gca,'Xtick',[2*pi,4*pi,6*pi],'Ytick',[0.95,1,1.05,max(y)]) grid on

title('\it y = 1 - e^{ -\alphat}cos{\omegat}')

text(13.5,1.2,'\fontsize{12}{\alpha}=0.3')

text(13.5,1.1,'\fontsize{12}{\omega}=0.7')

hold on;plot(ts,0.95,'bo','MarkerSize',10);hold off

cell_string{1}='\fontsize{12}\uparrow';

cell_string{2}='\fontsize{16} \fontname{隶书}镇定时间';

cell_string{3}='\fontsize{6} ';

cell_string{4}=['\fontsize{14}\rmt_{s} = ' num2str(ts)];

text(ts,0.85,cell_string)

xlabel('\fontsize{14} \bft \rightarrow')

图 6-7 二阶阶跃响应图用MATLAB4.x 和5.x 版标识时的差别

6.2.6多次叠绘

【例6-8】利用hold 绘制离散信号通过零阶保持器后产生的波形。

t=2*pi*(0:20)/20;y=cos(t).*exp(-0.4*t);

6.2.7双纵坐标图

【例6-9】受热压力容器的期望温度是120度，期望压力是0.25Mpa 。在同一张图上画出它们的阶跃响应曲线。

S1=tf([1 1],[1 3 2 1]); %温度的传递函数对象模型。 S2=tf(1,[1 1 1]); %压力的传递函数对象模型。 [Y1,T1]=step(S1);

%计算阶跃响应 [Y2,T2]=step(S2); %计算阶跃响应 6.2.8多子图

【例6-10】演示subplot 指令对图形窗的分割。

clf;t=(pi*(0:1000)/1000)';

y1=sin(t);y2=sin(10*t);y12=sin(t).*sin(10*t); subplot(2,2,1),plot(t,y1);axis([0,pi,-1,1]) subplot(2,2,2),plot(t,y2);axis([0,pi,-1,1]) subplot('position',[0.2,0.05,0.6,0.45])

6.2.9三维线图指令plot3

【例6-11】简单例题。

t=(0:0.02:2)*pi;x=sin(t);y=cos(t);z=cos(2*t);

','宝石')

6.2.10网线图、曲面图基本指令格式

【例6-12】用曲面图表现函数22y x z +=。

clf,x=-4:4;y=x;[X,Y]=meshgrid(x,y); %生成 x-y 坐标“格点”矩阵

Z=X.^2+Y.^2; %计算格点上的函数值 surf(X,Y,Z);hold on,colormap(hot)

图 6-12 曲面图和格点

6.2.11各种直方图bar, barh, bar3, bar3h

【例6-13】二维直方图有两种图型：垂直直方图和水平直方图。而每种图型又有两种表现模式：累计式：分组式。本例选其两种加以表现。

x=-2:2; %注意：自变量要单调变化

Y=[3,5,2,4,1;3,4,5,2,1;5,4,3,2,5]; %各因素的相对贡献份额

subplot(1,2,1),bar(x',Y','stacked') %“累计式”直方图

xlabel('x'),ylabel('\Sigma y'),colormap(cool)%控制直方图的用色

legend('因素A','因素B','因素C')

subplot(1,2,2),barh(x',Y','grouped') %“分组式”水平直方图

6.2.12饼图pie, pie3

【例6-14】饼图指令pie , pie3用来表示各元素占总和的百分数。该指令第二输入宗量为与第一宗量同长的0-1向量，1使对应扇块突出。

a=[1,1.6,1.2,0.8,2.1];

subplot(1,2,1),pie(a,[1 0 1 0 0]),legend({'1','2','3','4','5'})

6.3图形窗的图形编辑功能

【例6-15】运行以下指令

clf;shg,t=(pi*(0:1000)/1000)';y1=sin(t);y12=sin(t).*sin(10*t);

plot(t,y12,'b-',t,[y1,-y1],'r:'),axis([0,pi,-1,1])

图6-15 MATLAB 5.3版图形窗

【例6-16】编辑图形

图 6-16 经交互编辑后的图形

6.4函数绘图的简捷指令

【例6-17】绘制y e t t

=-233

2

2cos 和它的积分?=t dt t y t s 0)()(在]*4,0[pi 间的图形。

syms t tao;

y=2/3*exp(-t/2)*cos(sqrt(3)/2*t); %定义符号函数 s=subs(int(y,t,0,tao),tao,t); %获得积分函数

subplot(1,2,1),ezplot(y,[0,4*pi]);grid %一元函数简捷绘图指令 subplot(1,2,2),ezplot(s,[0,4*pi]);grid

【例6-18】画出2

2

22)

sin(y

x y x z ++=

所表示的三维曲面。y x ,的取值范围是]8,8[-。

解：clear; %清内存

x=-8:0.5:8; %定义自变量x 的一维刻度向量

y=x'; %定义自变量y 的一维刻度向量，“’”为矩阵转置 X=ones(size(y))*x; %计算自变量平面上取值点x 坐标的二维数组 Y=y*ones(size(x)); %计算自变量平面上取值点y 坐标的二维数组 R=sqrt(X.^2+Y.^2)+eps; %计算中间变量22y x R +=

Z=sin(R)./R;

%计算与自变量二维数组相应的函数值R

R

z sin =

mesh(Z); %绘制三维网格图

colormap(hot)

%指定网格图用hot 色图绘制

图 6-18 三维网线图

习题

1．绘制函数x xe y -=在10≤≤x 时的曲线。 2. 绘制函数y=sin(t) 的曲线。

3.编写程序作图：函数y=２*sin(3*x)+cos(x)，区间：0≤x ≤4π, -5≤y ≤5。（参考：fplot ）

4.利用peaks 函数产生矩阵x 、y 、z ，在一个图形窗中作出曲面图和等值线图（曲面图位于界面上部，等值线图位于界面下部）。（参考：peaks,subplot,surf,contour ）

潮流计算(matlab)实例计算

潮流例题：根据给定的参数或工程具体要求（如图），收集和查阅资料；学习相关软件（软件自选：本设计选择Matlab进行设计）。 2.在给定的电力网络上画出等值电路图。 3.运用计算机进行潮流计算。 4.编写设计说明书。 一、设计原理 1．牛顿-拉夫逊原理 牛顿迭代法是取x0 之后，在这个基础上，找到比x0 更接近的方程的跟，一步一步迭代，从而找到更接近方程根的近似跟。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0 的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。电力系统潮流计算，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡节点外）可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知节点赋电压初值，一般为额定电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、节点电压不平衡量（未知的）构成了误差方程，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新

的电压不平衡量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。 牛顿—拉夫逊迭代法的一般步骤： （1）形成各节点导纳矩阵Y。 （2）设个节点电压的初始值U和相角初始值e 还有迭代次数初值为0。 （3）计算各个节点的功率不平衡量。 （4）根据收敛条件判断是否满足，若不满足则向下进行。 （5）计算雅可比矩阵中的各元素。 （6）修正方程式个节点电压 （7）利用新值自第（3）步开始进入下一次迭代，直至达到精度退出循环。 （8）计算平衡节点输出功率和各线路功率 2．网络节点的优化 １）静态地按最少出线支路数编号 这种方法由称为静态优化法。在编号以前。首先统计电力网络个节点的出线支路数，然后，按出线支路数有少到多的节点顺序编号。当由n 个节点的出线支路相同时，则可以按任意次序对这n 个节点进行编号。这种编号方法的根据是导纳矩阵中，出线支路数最少的节点所对应的行中非零元素也２）动态地按增加出线支路数最少编号在上述的方法中，各节点的出线支路数是按原始网络统计出来的，在编号过程中认为固定不变的，事实上，在节点消去过程中，每消去一个节点以后，与该节点相连的各节点的出线支路数将发生变化(增加，减少或保持不变)。因此，如果每消去一个节点后，立即修正尚未编号节点的出线支路数，然后选其中支路数最少的一个节点进行编号，就可以预期得到更好的效果，动态按最少出线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中各节点出线支路数目的变动情况。 3．MATLAB编程应用 Matlab 是“Matrix Laboratory”的缩写，主要包括：一般数值分析，矩阵运算、数字信号处理、建模、系统控制、优化和图形显示等应用程序。由于使用Matlab 编程运算与人进行科学计算的思路和表达方式完全一致，所以不像学习高级语言那样难于掌握，而且编程效率和计算效率极高，还可在计算机上直接输出结果和精美的图形拷贝，所以它的确为一高效的科研助手。 二、设计内容 1.设计流程图

MATLAB精通科学计算_偏微分方程求解

一、Maple V 系统 Maple V是由Waterloo大学开发的数学系统软件，它不但具有精确的数值处理功能，而且具有无以伦比的符号计算功能。Maple V的符号计算能力还是MathCAD和MATLAB等软件的符号处理的核心。Maple提供了2000余种数学函数，涉及范围包括：普通数学、高等数学、线性代数、数论、离散数学、图形学。它还提供了一套内置的编程语言，用户可以开发自己的应用程序，而且Maple自身的2000多种函数，基本上是用此语言开发的。 Maple采用字符行输入方式，输入时需要按照规定的格式输入，虽然与一般常见的数学格式不同，但灵活方便，也很容易理解。输出则可以选择字符方式和图形方式，产生的图形结果可以很方便地剪贴到Windows应用程序内。 二、MATLAB 系统 MATLAB原是矩阵实验室（Matrix Laboratory）在70年代用来提供Linpack和Eispac k软件包的接口程序，采用C语言编写。从80年代出现3.0的DOS版本，逐渐成为科技计算、视图交互系统和程序语言。MATLAB可以运行在十几个操作平台上，比较常见的有基于W indows 9X/NT、OS/2、Macintosh、Sun、Unix、Linux等平台的系统。 MATLAB程序主要由主程序和各种工具包组成，其中主程序包含数百个内部核心函数，工具包则包括复杂系统仿真、信号处理工具包、系统识别工具包、优化工具包、神经网络工具包、控制系统工具包、μ分析和综合工具包、样条工具包、符号数学工具包、图像处理工具包、统计工具包等。而且5.x版本还包含一套几十个的PDF文件，从MATLAB的使用入门到其他专题应用均有详细的介绍。 MATLAB是数值计算的先锋，它以矩阵作为基本数据单位，在应用线性代数、数理统计、自动控制、数字信号处理、动态系统仿真方面已经成为首选工具，同时也是科研工作人员和大学生、研究生进行科学研究的得力工具。MATLAB在输入方面也很方便，可以使用内部的E ditor或者其他任何字符处理器，同时它还可以与Word6.0/7.0结合在一起，在Word的页面里直接调用MATLAB的大部分功能，使Word具有特殊的计算能力。 三、MathCAD 系统 MathCAD是美国Mathsoft公司推出的一个交互式的数学系统软件。从早期的DOS下的1. 0和Windows下的4.0版本，到今日的8.0版本，功能也从简单的数值计算，直至引用Map le强大的符号计算能力，使得它发生了一个质的飞跃。 MathCAD是集文本编辑、数学计算、程序编辑和仿真于一体的软件。MathCAD7.0 Profe ssional（专业版）运行在Win9X/NT下，它的主要特点是输入格式与人们习惯的数学书写格式很近似，采用WYSWYG（所见所得）界面，特别适合一般无须进行复杂编程或要求比较特殊的计算。MathCAD 7.0 Professional 还带有一个程序编辑器，对于一般比较短小，或者要求计算速度比较低时，采用它也是可以的。这个程序编辑器的优点是语法特别简单。 MathCAD可以看作是一个功能强大的计算器，没有很复杂的规则；同时它也可以和Wor d、Lotus、WPS2000等字处理软件很好地配合使用，可以把它当作一个出色的全屏幕数学公式编辑器。 四、Mathematica 系统 Mathematica是由美国物理学家Stephen Wolfram领导的Wolfram Research开发的数学系统软件。它拥有强大的数值计算和符号计算能力，在这一方面与Maple类似，但它的符

实验五_MATLAB计算的可视化

实验五 MATLAB 计算的可视化（一） 实验目的 1. 熟练掌握MATLAB 二维曲线的绘制 2.掌握图形的修饰 3.掌握三维图的绘制 4.了解各种特殊图形的绘制 内容与步骤 1.在同一幅图形窗口中分别绘制y1=sin(t)和y2=cos(t)二条函数曲线，t 的取值范围为[0,10]。y1用红色虚线表示，y2用蓝色实线表示，横坐标轴名称为“时间t ”，纵坐标轴名称为“正弦、余弦”，整个图形的标题为“正弦和余弦曲线”。在坐标(1.7*pi ，-0.3)处添加文字“sin(t)”， 在坐标(1.6*pi ，0.8)处添加文字“cos(t)”，并在右上角添加图例，其运行界面图如下图所示。之后并尝试修改坐标轴刻度。 2.用subplot 命令在同一个窗口的不同子窗口绘制曲线y=sin(t)，y1=sin(t+0.25) y2=sin(t+0.5)，其中t=[0 10]。 3.绘制三维曲线：?? ? ??=≤≤==)cos()sin()200() cos()sin(t t t z t t y t x π （注意：用plot3命令） 4.三维网线图：绘制z=sin(y)cos(x) 三维网线图。 5. 三维曲面图 绘制22y x z +=的三维曲面图，x 在[-5,5]范围，y 在[-5,5]范围。将曲面图颜色用shading 命令连续变化，并用颜色标尺显示色图（使用函数colorbar 生成）。生成的图形如下图所示。

6.请绘制一个饼形图，数据如下表所示 7. 用semilogx命令绘制传递函数为1//(s+1)(0.5s+1)的对数幅频特性曲线，横坐标为w，纵坐标为Lw，w的范围为10-2-103，按对数分布。

科学计算与MATLAB语言(第四课)

第四讲绘图功能

作为一个功能强大的工具软件，Matlab 具有很强的图形处理功能，提供了大量的二维、三维图形函数。由于系统采用面向对象的技术和丰富的矩阵运算，所以在图形处理方面即常方便又高效。

4.1 二维图形 一、plot函数 函数格式：plot(x,y)其中x和y为坐标向量函数功能：以向量x、y为轴，绘制曲线。【例1】在区间0≤X≤2 内，绘制正弦曲线Y=SIN（X），其程序为： x=0:pi/100:2*pi; y=sin(x); plot(x,y)

一、plot函数 【例2】同时绘制正、余弦两条曲线Y1=SIN（X）和Y2=COS（X），其程序为： x=0:pi/100:2*pi; y1=sin(x); y2=cos(x); plot(x,y1,x,y2) plot函数还可以为plot(x,y1,x,y2，x,y3，…)形式，其功能是以公共向量x为X轴，分别以y1，y2，y3，…为Y轴，在同一幅图内绘制出多条曲线。

一、plot函数 （一）线型与颜色 格式：plot(x,y1,’cs’,...) 其中c表示颜色，s表示线型。 【例3】用不同线型和颜色重新绘制例4.2图形，其程序为：x=0:pi/100:2*pi; y1=sin(x); y2=cos(x); plot(x,y1,'go',x,y2,'b-.') 其中参数'go'和'b-.'表示图形的颜色和线型。g表示绿色，o表示图形线型为圆圈；b表示蓝色，-.表示图形线型为点划线。

一、plot函数 （二）图形标记 在绘制图形的同时，可以对图形加上一些说明，如图形名称、图形某一部分的含义、坐标说明等，将这些操作称为添加图形标记。 title(‘加图形标题'); xlabel('加X轴标记'); ylabel('加Y轴标记'); text(X,Y,'添加文本');

流体力学简单计算MATLAB程式

用matlab进行编程计算 第一问： z=30;p1=50*9.8*10^4;p2=2*9.8*10^4;jdc=0.00015;gama=9800;d=0.257;L=50000 ;mu=6*10^(-6); hf=z+(p1-p2)/(0.86*gama) xdc=2*jdc/d; beta=4.15;m=1; Q=(hf*d^(5-m)/(beta*mu^m*L))^(1/(2-m)); v=4*Q/(pi*d^2); Re=v*d/mu; Re1=59.7/xdc^(8*xdc/7)； Re2=(665-765*log(xdc))/xdc； i=hf/L; if Re<3000 Q=Q; elseif 3000

《Matlab与科学计算》作业 2010010099

《Matlab与科学计算》作业 第一章MATLAB环境 1、MATLAB通用操作界面窗口包括哪些？命令窗口、历史命令窗口、当前目录窗口、工作空间窗口各有哪些功能？ 答：MATLAB通用操作界面窗口包括：命令窗口、历史命令窗口、当前目录浏览器窗口、工作空间窗口、变量编辑器窗口、M文件编辑/调试器窗口、程序性能剖析窗口、MATLAB帮助。 命令窗口是MATLAB命令操作的最主要窗口，可以把命令窗口当做高级的“草稿纸”。在命令窗口中可以输入各种MATLAB的命令、函数和表达式，并显示除图形外的所有运算结果。 历史命令窗口用来记录并显示已经运行过的命令、函数和表达式，并允许用户对它们进行选择、复制和重运行，用户可以方便地输入和修改命令，选择多行命令以产生M文件。 当前目录窗口用来设置当前目录，可以随时显示当前目录下的M、MKL等文件的信息，扬文件类型、文件名、最后个修改时间和文件的说明信息等，并可以复制、编辑和运行M文件及装载MAT数据文件。 工作空间窗口用来显示所有MATLAB工作空间中的变量名、数据结构、类型、大小和字节数。 2、熟悉课本中表格1.4、1.5、1.6、1.7、1.8的内容。 3、如何生成数据文件？如何把数据文件中的相关内容输入到工作空间中，用实例进行操作。 生成数据文件：

把数据文件中的相关内容输入到工作空间中： 结果： 4、在工作空间中可以通过哪些命令管理变量，写出每种语法的具体操作过程。答：（1）把工作空间中的数据存放到MAT数据文件。 语法：save filename 变量1 变量2 ……参数。 （2）从数据文件中取出变量存放到工作空间。 语法：load filename 变量1 变量2 ……。

《MATLAB与科学计算》期末论文

盐城师范学院《MATLAB与科学计算》期末论文 2016－2017学年度第一学期 用MATLAB解决解析几何的图形问题 学生姓名吴梦成 学院数学与统计学院 专业信息与计算科学 班级数15(5)信计 学号 15213542

用MATLAB 解决解析几何的图形问题 摘 要 将 MATLAB 的图形和动画功能都用于解析几何教学，可使教学形象生动。以图形问题为例，详细给出了实例的程序编写和动画实现过程 。在解析几何教学中有一定的应用价值。 【关键词】： MATLAB ； 解析几何 ；图形 ； 动 画；编程 1 引 言 在解析几何的教学中，使用传统的教学方法。许多曲线及曲面的形成过程与变换过程只通过传统的教师讲授静态图示就很难形象生动地表示出来 。在解析几何教学中使用MATLAB 软件辅助教学，不仅可以很容易绘制出复杂的立体图形。把曲线、曲面的形成和变化过程准确地模拟出来 ，而且还能够对它们进行翻转 、旋转 ，甚 至还能够轻而易举地实现图形的动画效果!这对提高教学效率和培养学生的空间想象能力可起到事半功倍的效果。下面结合实例从几个方面说明MATLAB 在解析几何画图方面的应用。 2 利用 MATLAB 绘制三维曲线 在空间解析几何中，各种曲线和曲面方程的建立都离不开图形 ，而空间曲线和曲面图形既难画又费时。借助MATLAB 的绘图功能 ，可以快捷 、 准确地绘出图形，使教学变得形象 、生动 。有利于学生观察三维空间图形的形状 ， 掌握图形的性质 。 一 般地 ，MATLAB 可用plot3，ezplot3，comet3等函数来各种三维曲线 。 例如画螺旋曲线的图形，其参数方程设为 ：t at cos x =，t b sin t y -=，ct =z 。使用 plot3语句画螺旋曲线图形的方法如下( 设a =2 ,b=4,c=3)： );*3),sin(*.*4),cos(*.*2(3;*10:50/:0t t t t t plot pi pi t -= MATLAB 用两条简单的语句就可以画出螺旋 曲线(图1)，但上述方法是静态的 ，为了体

实验4、matlab的计算可视化和GUI设计

p345 subplot(2,2,1) t1=0:0.1:2; y1=sin(2*pi*t1); plot(t1,y1); title('y=sin(2\pit)') 练习： subplot(2,2,2) t2=0:0.1:2; y2=[exp(-t2);exp(-2*t2);exp(-3*t2)]; plot(t2,y2) axis([0 2 -0.2 1.2]); title('y=e-t,y=e-2t,y=e-3t') 练习： subplot(2,2,3); t3=[0 1 1 2 2 3 4]; y3=[0 0 2 2 0 0 0]; plot(t3,y3); axis([0 4 -0.5 3]); title('脉冲信号') 练习： subplot(2,2,4); t4=0:0.1:2*pi; plot(sin(t4),cos(t4));

axis([-1.2 1.2 -1.2 1.2]); axis equal; title('圆') 练习： P346 x=0:0.1:20; zeta=0 y1=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt( 1-zeta^2)*x+acos(zeta)); plot(x,y1) zeta=0.3; y2=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt( 1-zeta^2)*x+acos(zeta)); hold on plot(x,y2,'r:') zeta=0.5; y3=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt( 1-zeta^2)*x+acos(zeta)); plot(x,y3,'g*') zeta=0.707; y4=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt( 1-zeta^2)*x+acos(zeta)); plot(x,y4,'m-') title('二阶系统曲线') legend('\zeta=0','\zeta=0.3','\zeta=0.5','\zeta=0. 707') grid on gtext('\zeta=0') gtext('\zeta=0.3') gtext('\zeta=0.5') gtext('\zeta=0.707') ginput(3) zeta = ans = 2.6037 0.9035 13.1106 2.0029 4.2166 1.0380 P347 h_fig=gcf h_axis=gca h_line1=gco h_title=get(gca,'title') h_text2=findobj(h_fig,'string','\zeta=0.3') h_fig = 1 h_axis = 151.0018 h_line1 = 1 h_title = 152.0018 h_text2 = Empty matrix: 0-by-1 set(h_line1,'linewidth',5)

Matlab与科学计算样题(加主观题答案)

Matlab 与科学计算考试样题（客观题） 1 下面的MATLAB 语句中正确的有： a) 2a ＝pi 。 b) record_1=3+4i c) a=2.0, d) c=1+6j 2. 已知水的黏度随温度的变化公式如下，其中a=0.03368，b=0.000221，计算温度t 为20，30，40度时的粘度分别是： 2 1at bt μμ=++0μ为0℃水的黏度，值为31.78510-?；a 、b 为常数，分别为0.03368、0.000221。 3. 请补充语句以画出如图所示的图形： [x,y]=meshgrid(-2:0.1:2, -2:0.1:2)。 Z=x.*exp(-x.^2-y.^2)。 。 a) Plot3(x,y,Z) b) plot3(x,y,Z) c) mesh(x,y,Z) d) plot3(x,y,z) 2 a) 0.4900 1.2501 0.8560 b) 0.8560 1.2501 0.4900 c) -0.6341 3.8189 -3.7749 d) 3.8189 -3.7749 2.8533 解释说明：

>> x=0.5:0.5:3.0。 >> y=[1.75,2.45,3.81,4.80,8.00,8.60]。 >> a=polyfit(x,y,2) a = 0.4900 1.2501 0.8560 >> x1=[0.5:0.25:3.0]。 >> y1=a(1)*x1.^2+a(2)*x1+a(3) >> plot(x,y,'*') >> hold on >> plot(x1,y1,'--r') 5. 求方程在 x=0.5附近的根. 21 x x += a) 0.6180 b) -1.1719e-25 c) －1 d) -1.6180 6. 用Newton-Cotes方法计算如下积分 1 5 x? （a）133.6625 (b)23.8600 (c) 87.9027 (d) -1.6180 7. y=ln(1+x)，求x=1时y" a) -0.25 b) 0.5 c) -0.6137 d) -1.6137 8.某公司用3台轧机来生产规格相同的铝合金薄板。取样测量薄板的 厚度，精确至‰厘M。得结果如下： 轧机1：0.236 0.238 0.248 0.245 0.243 轧机2：0.257 0.253 0.255 0.254 0.261 轧机3：0.258 0.264 0.259 0.267 0.262 计算方差分析结果，并判定各台轧机所生产的薄板的厚度有无显著的差异？ a) p＝1.3431e-005，没有显著差异。

matlab计算结果的可视化

第五讲计算结果的可视化 本节介绍MATLAB 的两种基本绘图功能：二维平面图形和三维立体图形。 5.1 二维平面图形 5.1.1 基本图形函数 plot 是绘制二维图形的最基本函数，它是针对向量或矩阵的列来绘制曲线的。也就是说，使用plot 函数之前，必须首先定义好曲线上每一点的x 及y 坐标，常用格式为：（1）plot(x) 当x 为一向量时，以x 元素的值为纵坐标，x 的序号为横坐标值绘制 曲线。当x 为一实矩阵时，则以其序号为横坐标，按列绘制每列元素值相对于其序号的曲线， 当x 为m× n 矩阵时，就由n 条曲线。 （2）plot(x,y) 以x 元素为横坐标值，y 元素为纵坐标值绘制曲线。 （3）plot(x,y1,x,y2,…) 以公共的x 元素为横坐标值，以y1,y2,… 元素为纵坐标值绘 制多条曲线。 例5.1.1 画出一条正弦曲线和一条余弦曲线。 >> x=0:pi/10:2*pi; >> y1=sin(x); >> y2=cos(x); >> plot(x,y1,x,y2) 图5.1.1 函数plot 绘制的正弦曲线 在绘制曲线图形时，常常采用多种颜色或线型来区分不同的数据组，MATLAB 软件专 门提供了这方面的参数选项（见表5.1.1），我们只要在每个坐标后加上相关字符串，就可实 现它们的功能。 - 2 - 表5.1.1 绘图参数表 色彩字符颜色线型字符线型格式标记符号数据点形式标记符号数据点形式 y 黄- 实线. 点<小于号 m 紫：点线o 圆s 正方形 c 青-. 点划线x 叉号 d 菱形 r 红- - 虚线+ 加号h 六角星 g 绿* 星号p 五角星 b 蓝v 向下三角形 w 白^ 向上三角形 k 黑>大于号 例如，在上例中输入 >> plot(x,y1,'r+-',x,y2,'k*:') 则得图5.1.2 图5.1.2 使用不同标记的plot 函数绘制的正弦曲线 5.1.2 图形修饰 MATLAB 软件为用户提供了一些特殊的图形函数，用于修饰已经绘制好的图形。 表5.1.2 图形修饰函数表

中国矿业大学 实验六 MATLAB数据可视化

实验六MATLAB数据可视化 一、实验目的 掌握MATLAB 二维、三维图形绘制，掌握图形属性的设置和图形修饰；掌握图像文件的读取和显示。 二、实验内容 (1) 二维图形绘制。 (2) 三维曲线和三维曲面绘制。 三、实验步骤 1．二维图形绘制 (1) 二维图形绘制主要使用函数plot。 >> clear all; >> x=linspace(0,2*pi,100); >> y1=sin(x); >> plot(x,y) >> hold on >> y2=cos(x) >> plot(x,y) >> hold off

注：hold on 用于保持图形窗口中原有的图形，hold off解除保持。 (2) 函数plot 的参数也可以是矩阵。 >> close all >> x=linspace(0,2*pi,100); >> y1=sin(x); >> y2=cos(x); >> A=[y1 ; y2]'; >> B=[x ; x]' >> plot(B,A)

(3) 选用绘图线形和颜色。>> close all >> plot(x,y1,'g+',x,y2, 'r:') >> grid on

(4) 添加文字标注。 >> title('正弦曲线和余弦曲线') >> ylabel('幅度') >> xlabel('时间') >> legend('sin(x)', 'cos(x)') >> gtext('\leftarrowsinx')

(5) 修改坐标轴范围。 >> axis equal >> axis normal >> axis([0 pi 0 1.5]) 程序如下： x=linspace(0,2*pi,100); y1=sin(x); y2=cos(x); A=[y1 ; y2]'; B=[x ; x]' plot(B,A) plot(x,y1,'g+',x,y2, 'r:') axis equal axis normal axis([0 pi 0 1.5])

Matlab与工程计算 第二章 Matlab矩阵及其运算

第2章Matlab矩阵及其运算 2.1 Matlab变量 2.2 Matlab数值矩阵 2.3 运算符 2.4 基本数学函数 2.5 稀疏矩阵 2.6 矩阵分析 2.8 字符串 2.9 结构数据 2.10 细胞矩阵

2.1 Matlab变量 1. 变量命名规则 在MATLAB 6.5中，变量名是以字母开头，后接字母、数字或下划线的字符序列，最多63个字符。在MATLAB中，变量名区分字母的大小写。 2．变量赋值 变量=表达式 3．预定义变量 i，j，pi，eps，realmin，realmax，inf，NaN 预定义变量有特定的含义，在使用时，应尽量避免对这些变量重新赋值。

内存变量的管理 1．指令操作法 who whos clear 2. 现场菜单操作法 3. 内存变量文件(.mat) save [文件名] [变量名表] [-append][-ascii] load [文件名] [变量名表] [-ascii] help save help load

数值数据的输出格式 MATLAB用十进制数表示一个常数，具体可采用日常记数法和科学记数法两种表示方法。 在一般情况下，MATLAB内部每一个数值数据元素都是用双精度数来表示和存储的。 数据输出时用户可以用format命令设置或改变数据输出格式。format 命令的格式为： format格式符 其中格式符决定数据的输出格式 help format

2.2 MATLAB数值矩阵 2.2.1 矩阵的建立 1．直接输入法 最简单的建立矩阵的方法是从键盘直接输入矩阵的元素。具体方法如下：将矩阵的元素用方括号括起来，按矩阵行的顺序输入各元素，同一行的各元素之间用空格或逗号分隔，不同行的元素之间用分号分隔。

科学计算与MATLAB

单元测验已完成成绩：100.0分 1 【单选题】 MATLAB一词来自（）的缩写。 ?A、 Mathematica Laboratory ?B、 Matrix Laboratory ?C、 MathWorks Lab ?D、 Matrices Lab 我的答案：B得分：50.0分 2 【单选题】 下列选项中能反应MATLAB特点的是（）。 ?A、 算法最优 ?B、 不需要写程序 ?C、 程序执行效率高 . .

?D、 编程效率高 我的答案：D得分：50.0分 单元测验已完成成绩：96.4分 1 【单选题】 当在命令行窗口执行命令时，如果不想立即在命令行窗口中输出结果，可以在命令后加上（）。 ?A、 冒号（:） ?B、 逗号（,） ?C、 分号（;） ?D、 百分号（%） 我的答案：C得分：7.1分 2 【单选题】 fix(264/100)+mod(264,10)*10的值是（）。 ?A、 86 . .

?B、 62 ?C、 423 ?D、 42 我的答案：D得分：7.1分 3 【单选题】 在命令行窗口输入下列命令后，x的值是（）。 >> clear >> x=i*j ?A、 不确定 ?B、 -1 ?C、 1 ?D、 i*j 我的答案：B得分：7.1分 . .

4 【单选题】 使用语句x=linspace(0,pi,6)生成的是（）个元素的向量。 ?A、 8 ?B、 7 ?C、 6 ?D、 5 我的答案：C得分：7.1分 5 【单选题】 ceil(-2.1)的结果为（）。 ?A、 -2 ?B、 -3 ?C、 1 . .

?D、 2 我的答案：A得分：7.1分 6 【单选题】 eval('sqrt(4)+2')的值是（）。 ?A、 sqrt(4)+2 ?B、 4 ?C、 2 ?D、 2+2 我的答案：B得分：7.1分 7 【单选题】 已知a为3×5矩阵，则执行完a(:,[2,4])=[]后（）。 ?A、 a变成行向量 ?B、 a变为3行2列 . .

【原创】MATLAB实验报告-第二次-用MATLAB实现计算数据可视化-北京交通大学

MATLAB 上机实验报告( 2 ) 实验内容: 一、试用如下几种方法来建立向量，观察结果 ( 1) x=1:5, x=(1:5) ' 实验结果：x=1:5是行向量,x=(1:5)是列向量.且1为初始值，5为终止值，默认的步长为 1. >> x=1:5 1 2 3 4 5 >> x=(1:5)' x = 1 2

3 4 5 ( 2) x=0:pi/4:pi 实验结果：x=0:pi/4:pi 指的是x=(0,0.25*pi,0.50*pi,0.75*pi,pi). 其中pi为圆周率初始值为0,终止值为pi,步长为pi/4. >> x=0:pi/4:pi x = 0 0.7854 1.5708 2.3562 3.1416 (3)x=(0：0.2：3) ', y=e-x)p.(*sin(x) 实验结果：x的初始值为0,终止值为3,步长为0.2.而函数y表示将x向量中的每一个数代入函数y=e%x)*sin(x)得到的函数值组成的向量. >> x=(0:0.2:3)', y=exp(-x).*sin(x)

x = 0.2000 0.4000 0.6000 0.8000 1.0000 1.2000 1.4000 1.6000 1.8000 2.0000 2.2000 2.4000 2.6000 2.8000 3.0000

0.1627 0.2610 0.3099 0.3223 0.3096 0.2807 0.2430 0.2018 0.1610 0.1231 0.0896 0.0613 0.0383 0.0204 0.0070 (4) k=linspace(-pi,pi,5), k=logspace(-3,-1,5) 实验结果：k=linspace(-pi,pi,5)产生的是初始值为-pi，终止值为 pi,元素总数为5的行向量，即k的步长为pi/2. k=logspace(-3,-1,5产生的是初始值为10八(-3),终止值为10八(-1),元素总数为5的列向量.

实验二MATLAB计算的可视化

课程实验报告 学年学期2011-2012学年第1学期 课程名称MATLAB与科学计算 实验名称实验二MATLAB计算的可视化实验室测量测绘实验中心计算机室专业年级热动113 学生姓名白治朋 学生学号2011012106 提交时间2013年10月23日 成绩 任课教师许景辉、牛亚斌 水利与建筑工程学院

实验二 MATLAB 计算的可视化 1、目的和要求 （1）熟练掌握MATLAB 二维曲线、三维图形的绘制。 （2）熟练掌握各种特殊图形的绘制。 （3）熟练掌握三维图形绘制命令。 （4）了解GUI 设计的一般过程和方法。 2、内容和步骤 参见教材实验四。 3、实验报告提交要求 （1） x=[1 2 3]，y=[1 2;2 3;5 8]，z=[2 6 9;3 8 8;1 5 7]，绘制plot （x ，y ）、plot （x ，z ），说明其各 自绘制的内容。 （2） 绘制如下图形，建立figure （2），绘图同样曲线，但标题为“你的姓名（黑体，16号字）”， 在x 坐标和y 坐标上分别标识学号和班级名称，并将网格线打开。 数组X 的列个数与矩阵y 的行个数相同， plot （ x ，y ）绘制的是x 为横坐标y 的每列为纵坐标的图像。如图1。 图1 数组X 的列个数与方阵z 的行列个数相同，plot （x ，z ）绘制的是x 为横坐标z 的每列为纵坐标的图像。如图2。 图2

（3）演示P133页，例题4.17 。

（4）完成课本P336图S 4.1实验，并用.m文件显示其程序内容。 （5）完成P302第四章例题4.

（6）通过绘制二阶系统阶跃响应，综合演示图形标识，请注释每条命令的含义。 clf; %清除图形窗口 t=6*pi*(0:100)/100;y=1-exp(-0.3*t).*cos(0.7*t); % 数据准备 tt=t(find(abs(y-1)>0.05)); %找出符合条件(y-1)的绝对值>0.05的对应t，赋值给tt ts=max(tt); %ts为tt中最大值ts=9.6133 plot(t,y,'r-','LineWidth',3) %画曲线t-y,红色实线，线粗3磅 axis([-inf,6*pi,0.6,inf]) %设置坐标轴范围。x轴下限自动产生,上限为6*pi；y轴下限0.6,上限自动产生 set(gca,'Xtick',[2*pi,4*pi,6*pi],'Ytick',[0.95,1,1.05,max(y)]) %二维坐标刻度设置。x轴刻度线取2*pi,4*pi,6*pi,y轴取0.95,1,1.05,max(y) grid on %显示坐标刻度线 title('\it y = 1 - e^{ -\alphat}cos{\omegat}') %用斜体1书写图名 text(13.5,1.2,'\fontsize{12}{\alpha}=0.3') %图形标识，添加文字注释。在x=13.5,y=1.2处,字体大小12磅，标注α=0.3 text(13.5,1.1,'\fontsize{12}{\omega}=0.7') %图形标识，添加文字注释。在x=13.5,y=1.1处,字体大小12磅，标注ω=0.7 hold on; %保持原有图形 plot(ts,0.95,'bo','MarkerSize',10); %在x=ts，y=0.95处画蓝色的空心圆圈,大小为

第6章MATLAB计算结果可视化讲解

第六章MATLAB 计算结果可视化 6.1连续函数和离散函数的可视化 【例6-1】用图形表示离散函数1 ) 6(--=n y 。 n=0:12; %产生一组自变量数据 y=1./abs(n-6); %计算相应点的函数值 plot(n,y,'r*','MarkerSize',20) %用红花标出数据点 grid on %画坐标方格 【例6-2】用图形表示连续调制波形)9sin()sin(t t y =。 t1=(0:11)/11*pi; y1=sin(t1).*sin(9*t1); t2=(0:100)/100*pi; y2=sin(t2).*sin(9*t2); subplot(2,2,1),plot(t1,y1,'r.'),axis([0,pi,-1,1]),title('子图 (1)') subplot(2,2,2),plot(t2,y2,'r.'),axis([0,pi,-1,1]),title('子图 (2)') subplot(2,2,3),plot(t1,y1,t1,y1,'r.') axis([0,pi,-1,1]),title('子图 (3)') subplot(2,2,4),plot(t2,y2)

6.2二维曲线绘图的基本操作 6.2.1 plot 的基本调用格式 【例6-3】用图形表示连续调制波形)9sin()sin(t t y 及其包络线。 t=(0:pi/100:pi)'; %长度为101的时间采样列向量 y1=sin(t)*[1,-1]; %包络线函数值，是（101x2）的矩阵 y2=sin(t).*sin(9*t); %长度为101的调制波列向量 t3=pi*(0:9)/9; y3=sin(t3).*sin(9*t3);plot(t,y1,'r:',t,y2,'b',t3,y3,'bo') 6.2.2泛函绘图指令fplot 【例6-4】fplot 与一般绘图指令的绘图效果比较。 [x,y]=fplot('cos(tan(pi*x))',[-0.4,1.4],0.2e-3);n=length(x); subplot(1,2,1),plot(x,y) title('\fontsize{20}\fontname{隶书}泛函绘图指令效果') t=(-0.4:1.8/n:1.4)'; subplot(1,2,2),plot(t,cos(tan(pi*t))) 6.2.3曲线的色彩、线型和数据点形 【例6-5】用图形演示平面上一个方块四个顶点在仿射投影（Affine Projection ）下的位置、 形状变化。 %平面上的四个点和它们构成的方块

MATLAB可靠度计算程序(附解释)

应用MATLAB的可靠度计算程序（附语句解释） 以下是matlab中求可靠度计算程序以及各语句的解释，对可靠度编程入门的童鞋可做参考：function bbeta=LJF(muX,sigmaX)%子函数的settingup sLn=sqrt(log(1+(sigmaX(1)/muX(1))^2));mLn=log(muX(1))-sLn^2/2;%求dltLnR和muLnR aEv=sqrt(6)*sigmaX(3)/pi;uEv= -psi(1)*aEv-muX(3);%求极值I型分布的两个参数 muX1=muX;sigmaX1=sigmaX; x=muX;normX=eps; while abs(norm(x)-normX)/normX > 1e-6%精度控制 normX=norm(x); g=x(1)-x(2)-x(3);%极限方程 gX=[1;-1;-1]; cdfX=[logncdf(x(1),mLn,sLn);1-evcdf(-x(3),uEv,aEv)];%分布函数的当量正态化 pdfX=[lognpdf(x(1),mLn,sLn);evpdf(-x(3),uEv,aEv)];%分布密度的当量正态化 nc=norminv(cdfX);%求标准正态分布函数的反函数 sigmaX1(1:2:3)=normpdf(nc)./pdfX;%求得等效正态变量的均方差 muX1(1:2:3)=[x(1:2:3)-nc.*sigmaX1(1:2:3)];%求得等效正态变量的均值 gs=gX.*sigmaX1;alphaX= -gs/norm(gs);%求得敏感度系数cos(sita) bbeta=(g+gX'*(muX1-x))/norm(gs);%求得可靠度指标 x=muX1+bbeta*sigmaX1.*alphaX;%将结果迭代得新的x end 关于可靠度程序的参考书籍，具体可参考张明先生编著的《结构可靠度分析》。 作为一本结构可靠度分析方法的专著，本书系统介绍了结构可靠度分析的基本概念和主要方法。首先简要介绍了结构随机可靠度的基本概念以明确可靠度分析的目的和意义；之后对结构可靠度分析的几种重要方法进行了详细的阐述，包括一次二阶矩方法、二次二阶矩方法、二次四阶矩方法、渐近积分方法、响应面方法、Monte Carlo方法，还研究了结构体系可靠度分析方法、基于人工神经网络的结构可靠度分析方法；最后对结构模糊随机可靠度分析方法作了阐述。每章介绍一类方法，每节介绍其中一个相对独立的方法，对于每一个方法均给出了典型的例题和用MATLAB软件编写的计算机程序。附录按照字母顺序列出了本书程序中所采用的标识符和MATLAB函数，以方便读者阅读本书和使用书中的程序。 兼顾可靠度方法理论和方法实施，理论和实践并举，是本书的重要特色，也使本书具有很强的实用性。本书可供科技工作者、大专院校教师、研究生和高年级本科生使用，也可供工程技术人员参考使用。

MATLAB可视化方法和技巧1_3复数的计算和图示

复数的计算和图示 表3 MATLAB关于复数运算的函数 {范例3_1}复数的加减法 设有两个复数z1 = 1 + 2i和z2 = 4 + 3i，其中i是虚数单位。求两个复数的和z1 + z2和差z2–z1。 [解析]复数有三种表示形式 (1)代数式 z = x + i y(3_1_1) (2)三角式 z = r(cosθ + isinθ) (3_1_2) 其中r是复数的模，θ是复角。代数式与三角式的换算关系是 r=arctan y θ(3_1_3) x x = r cosθ，y = r sinθ(3_1_4) (3)指数式 z = r e iθ(3_1_5) 其中利用了欧拉公式 e iθ= cosθ + isinθ(3_1_6) 设有两个复数 z1 = x1 + i y1，z2 = x2 + i y2(3_1_7) 复数加法是 z= z1 + z2 = (x1 + x2) + i(y1 + y2) (3_1_8) 复数减法是 z= z1 - z2 = (x1 - x2) + i(y1 - y2) (3_1_9) [程序]P3_1plus.m如下。 %复数的加减法 clear %清除变量 z1=1+2i; %第1个复数(1) x1=real(z1); %取第1个复数的实部(2) y1=imag(z1); %取第1个复数的虚部(2) x2=4; %第2个复数的实部 y2=3; %第2个复数的虚部 z2=x2+i*y2; %形成第2个复数(3) z=z1+z2; %两复数之和(4) x=real(z); %取复数的实部 y=imag(z); %取复数的虚部 figure %创建图形窗口 quiver(0,0,x1,y1,0) %在复平面画第1个复数(5)