第4章--MATLAB的绘图.pptVIP

第4章 MATLAB绘图 已知周期方波信号傅里叶展开，用matlab验证，其中，周期为2，即w*t=2*pi,w=3.14 思路：逐步增加w项，叠加信号，使信号趋近于方波信号。 思考(作业p73)：对定义于0到1的方波，其它区间为0，用傅里叶展开来验证。 clear t0=0:7;xz=zeros(size(t0)); t1=[0,0,1,1,1,2,2,2,3,3,3,4,4,4,5,5,5,6,6,6]; ft=[0,1,1,0,-1,-1,0,1,1,0,-1,-1,0,1,1,0,-1,-1,0,1]; subplot(2,1,1) plot(t1,ft,LineWidth,2); hold on plot(t0,xz) axis([0,7,-1.2,1.2]); grid on t=0:0.01:7; sumb=zeros(size(t)); for k=1:200 sumb=sumb+4/pi*sin(3.14*(2*k-1)*t)/(2*k-1); subplot(2,1,2) hold on plot(t,sumb) grid on end subplot(2,1,1) hold on plot(t,sumb,r); clf reset % 清图，所有对象属性设为缺省值 t=(0:100)/100*2*pi; H_line=line(Xdata,t, Ydata,sin(t)) % 以t为横坐标数据，sin(t)为纵坐标数据画一条曲线，并返回该曲线的句柄 text(pi,0.8, \fontsize{14}sin(t)); % 在坐标为[pi,0.8]处添加14磅大小的字符串sin(t) H_c=get(get(H_line, Parent), Children) % 轴之所有“子”对象的句柄 T=get(H_c, Type) /////////////////////// clf reset t=(0:pi/100:2*pi); tt=t*[1,1]; yy=sin(tt)*diag([0.5,1]); plot(tt,yy); Hb=findobj(gca, Color, b) h1=line([0:10],[0:10]); % 创建一个线对象，并返回该对象的句柄值h1 set(h1, Color,[1 0 0]); % 设置线的颜色为红色 pause for t=1:10 set(h1, LineWidth, t); % 设置线宽为l pause(1) % 暂停1s end h2=gca; % 获得当前轴的句柄并赋值给h2 pause set(h2, Visible, off) % 将轴设置为不可见 x=0:2*pi/180:2*pi; y1=sin(x); y2=cos(x); plot(x,y1,x,y2,b); pause h=get(gca,Children); % 获得当前轴对象下所有子对象（sin(x)及cos(x)曲线）的句柄 yy=get(h(2), Ydata); % 取得第二个子对象，即sin(x)曲线的纵坐标数据，赋值给yy yy2=yy/2; set(h(1), Ydata,yy2); % 将第一个子对象，即cos(x)曲线的纵坐标数据设置为sin(x)/2 set(h(2), LineStyle, :); % 将sin(x)曲线的线型设置为虚线 已知m=1kg, L=1m,b=0.1kg/m.s,g=10（其中b为阻力系数）,根据力学原理分析写出 t=0～10s内单摆的位移角度θ(注意θ的方向，第一象限为正)的变化规律的表达式（微分方程）。 有公式-mg*sin(θ)-b*f=ma即-mg*sin(θ)-b*v=ma 在角度很小的情况下sin(θ)≈θ，(也可考虑不化简)则上式可以写成-mg*θ-b*v=ma v=L*dθ，a=d^2θ，所以最终公式为-mg*θ-b*L*dθ=m*L*d^2θ,再代入 m=1kg, L=1m,b=0.1kg/m.s,g=10 function xs=nonel