资源描述：
《实验二 绘图及简单科学计算》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、实验二绘图及简单科学计算主要内容：1)绘图简介2)科学计算，曲线拟合与插值绘图简介MATLAB不但擅长于矩阵相关的数值运算，也适合用于各种科学可视化（Scientificvisualization）。下面介绍MATLAB基本二维和三维的各项绘图命令，包含一维曲线及二维曲面的绘制、打印及保存。一、二维绘图：lplot是绘制一维曲线的基本函数，但在使用此函数之前，我们需先定义曲线上每一点的x及y坐标。下例可画出一条正弦曲线：»x=linspace(0,2*pi,100);%100个点的x坐标»y=sin(x);%对应的y坐标»plot(x,y);l
2、若要画出多条曲线，只需将坐标对依次放入plot函数即可：»plot(x,sin(x),x,cos(x));l若要改变颜色，在坐标对后面加上相关字符串即可：»plot(x,sin(x),'c',x,cos(x),'g');l若要同时改变颜色及线型（Linestyle），也是在坐标对后面加上相关字符串即可：»plot(x,sin(x),'co',x,cos(x),'g*');plot绘图函数的参数字符颜色字符图线型态黄色.点k黑色o圆w白色xxb蓝色++g绿色**r红色-实线c亮青色:点线m锰紫色-.点虚线--虚线l图形完成后，我们可用axis([
3、xmin,xmax,ymin,ymax])函数来调整坐标轴的范围：»axis([0,6,-1.2,1.2]);lMATLAB也可对图形加上各种注解与处理：»xlabel('InputValue');%x轴注解»ylabel('FunctionValue');%y轴注解»title('TwoTrigonometricFunctions');%图形标题»legend('y=sin(x)','y=cos(x)');%图形注解»gridon;%显示格线l用subplot来同时画出数个小图形于同一个窗口之中：»subplot(2,2,1);plot(x,
4、sin(x));%把窗口分成2*2个子窗口,在第一个子窗口绘图»subplot(2,2,2);plot(x,cos(x));%在第二个子窗口绘图»subplot(2,2,3);plot(x,sinh(x));%在第三个子窗口绘图»%在第四个子窗口绘图lMATLAB还有其他各种二维绘图函数，以适合不同的应用，详见下表。Bar长条图Errorbar图形加上误差范围Fplot较精确的函数图形Polar极坐标图Hist累计图Rose极坐标累计图Stairs阶梯图以下我们针对每个函数举例。l当数据点数量不多时，条形图是很适合的表示方式：»closeall
5、;%关闭所有的图形窗口»x=1:10;»y=rand(size(x));»bar(x,y);柱状图l如果已知数据的误差量，就可用errorbar来表示。下例以单位标准差来做数据的误差量：»x=linspace(0,2*pi,30);»y=sin(x);»e=std(y)*ones(size(x));»errorbar(x,y,e)l对于变化剧烈的函数，可用fplot来进行较精确的绘图，会对剧烈变化处进行较密集的取样，如下例：»fplot('sin(1/x)',[0.020.2]);%[0.020.2]是绘图范围l若要产生极坐标图形，可用pola
6、r：»theta=linspace(0,2*pi);»r=cos(4*theta);»polar(theta,r);l对于大量的数据，我们可用hist来显示数据的分布情况和统计特性。下面几个命令可用来验证randn产生的高斯随随机数分布：»x=randn(5000,1);%产生5000个m=0，s=1的高斯随机数»hist(x,20);%20代表长条的个数lrose和hist很接近，只不过是将数据大小视为角度，数据个数视为距离，并用极坐标绘制表示：»x=randn(1000,1);»rose(x);lstairs可画出阶梯图：»x=linspa
7、ce(0,10,50);»y=sin(x).*exp(-x/3);»stairs(x,y);1、y=[1234;3456;2358];plot(y)2、x=[2345];y=[1234;3456;2358];plot(x,y)1、练习：还有其他可能情况，请同学们实验2、holdonholdoff介绍3、练习：利用已学知识绘制分段函数图像;4、图轴控制：（1）画出两个刻度不同的y轴x=linspace(0,2*pi);y1=sin(x);y2=exp(-x);plotyy(x,y1,x,y2)（2）axis([xmin,xmax,ymin,yma
8、x])(3)plot(x,y1,'--',x,y2,':o')xlabel('x轴');ylabel('y轴')title('sin(x)与e^{-x
