MALTAB DAY 6

clear; Y=rand(1,8); h=plot(Y); syms x;%符号变量必须先写符号表达式，在给自变量赋值 y1=atan(x); y2=diff(y1); x=[-5:0.1:5]; y1=eval(y1);%函数求符号变量的值 y2=1./(1+x.^2);%计算法给符号变量赋值%一个参数必须为方阵，另一个必须为标量。请使用 POWER (.^) 执行按元素求幂 plot(x,y1,’-’,x,y2,’’);%绘制多个图的方式一，较有局限 xlabel(‘a’);%坐标轴标注 ylabel(‘b’); title(‘反正切函数及其微分曲线’); legend(’-atan,diff cure’);%设置图例，实线现实反正切函数
%绘制多个图的方式二 注本例全用函数求符号法来计算 clear; syms x; y=tan(sin(x)); y1=sin(tan(x)); x=[-pi:pi/10:pi]; y=eval(y); y1=eval(y1); plot(x,y,’–rs’,‘linewidth’,2,‘markeredgecolor’,‘k’,‘markerfacecolor’,‘g’,‘markersize’,10); hold on; plot(x,y1,’-ko’,‘linewidth’,2,‘markeredgecolor’,‘Y’,‘markerfacecolor’,‘b’,‘markersize’,10); legend(’–rstan(sin(x)),-kosin(tan(x))’);
%eg5_6 clear; x=[0:0.02:5]; y(1,251)=0; for i=1:101 y(i)=1; end for i=102:251 y(i)=0.2; end h=plot(x,y,’-r’); namearry={‘marker’,‘markersize’}; namevalue={‘o’,1}; set(h,‘marker’,‘o’,‘markersize’,1);%属性设置方法1 set(h,namearry,namevalue)%属性设置方法2 h%在命令行中显示所有属性 a=find(x>=2.0)%将x中》2.0的下表输出到a中
syms x,y; f=x.^2/9+y.2/4-1; %ezplot(‘f’);%注意‘’里面必须为符号函数表达式 h=ezplot(‘x.^2/9+y.2/4-1’); grid on; axis([-3,3,-2,2])%设置坐标轴范围 set(h,‘color’,‘r’,‘linestyle’,’–’);
clear; X=[1.55,1.58,1.6,1.64,1.66,1.68,1.7,1.73,1.78,1.8,1.82,1.85]; Y=[60,57,57,65,63,64,70,65,68,76,72,78]; plot(X,Y,‘b’); hold on; A=polyfit(X,Y,1); syms x y=A[x;1]; x=[1.5:0.1:2]; y=eval(y); plot(x,y,’-r’); grid on; xlabel(‘身高’); ylabel(‘体重’); title(‘身高与体重的关系’)
clear; A=[0.25,0.5,0.75,1.00,2,3,4,5,6,6.5,7;4,3,2,1,0,0,1,0,1,3,6]; t=A(1,:); faunum=A(2,:); t1=0:0.1:7; fauspl=interp1(t,faunum,t1,‘spline’);%用三次样条插值连接故障点 注意interp1（数字） plot(t,faunum,'r’,t1,fauspl,’-b’);%spline为样条曲线
%运用scatter创建散点图 clear; for r=2:20 zeta=0:pi/12:2pi; x=rcos(zeta+rpi/30); y=rsin(zeta+rpi/30); %设置颜色向量【小，大，大】 c(1)=r/20; c(2)=r/20; c(3)=r/20;%颜色会越来越淡 scatter(x,y,15,c,‘field’); hold on end
clear; zeta=0:6pi/600:6pi; r=sin(3zeta); [x,y]=pol2cart(zeta,r); subplot(1,2,1);polar(zeta,r);%将坐标分成1行2列，在左侧绘制极坐标形式的图 subplot(1,2,2);plot(x,y);%在右侧会直角坐标的图 grid on;
clear; Y=magic(4); subplot(1,2,1); bar(Y,1,‘b’); subplot(1,2,2); bar(Y,0.2,‘r’,‘stacked’);%叠加形式的条形图 %barh为水平绘制图形
clear; A=[22,20,33,41;35,19,74,84;81,60,44,52;13,19,46,67]; S=sum(A); [m,i]=max(S);%m为最大值，i为下标 explode=zeros(size(S)); explode(i)=1;%构造分离向量explode pie(S,explode); legend(‘产品1’,‘产品2’,‘产品3’,‘产品4’)%为颜色进行注释 %分左右显示 subplot(1,2,1); pie(S(1:2),explode(1:2)); title(‘12’); subplot(1,2,2); pie(S(3:4),explode(3:4)); title(‘34’);
%对于绘制时间序列样本数据是十分有用的、 X=[800,1000,1200,1500,2000,2500,3500,5000]; Y=[10,13,21,33,12,8,2,1]; h=stairs(X,Y,’–r’); set(h,‘linewidth’,5,‘marker’,‘o’);
%stem(X,Y)当Y为矩阵时，X的长度要与Y的行数相等 clear; A=[5,5,5,5,0,0,0,0,0,0,0,0]; B=[0,0,0,0,1.5,2.5,5,6,6.5,7,7,7]; subplot(2,1,1);axis([0 12 0 8]);stem(A,‘filled’);hold on;%圆点填满 subplot(2,1,2);stem(B,‘filled’);grid on;