矩阵的乘法 已知A1 - 副本

三维图形的绘制 %三维螺旋曲线 t=0:pi/100:10*pi; y1=sin(t);y2=cos(t); figure(4) plot3(y1,y2,t);

title('helix'), %给图形加标题

text(0,0,0,'origin'); %在指定位置创建说明性文字 xlabel('sin(t)'),ylabel('cos(t)'),zlabel('t');

grid; %在所画出的图形坐标中添加栅格，注意用在plot之后

三维网格曲面图 x=[0:0.05:2*pi]; y=[0:0.05:2*pi];

z=sin(y')*cos(x); %矩阵相乘 figure(5) mesh(x,y,z);

实验二

输入典型二阶系统参数，确定系统传递函数模型

T=input('please input T='); T1=input('please input T1='); K1=input('please input K1='); num=K1; den=[T*T1 T K1]; G=tf(num,den)

%判断系统是否稳定 [z,p,k]=tf2zp(num,den); i=find(real(p)>0); n=length(i);

if(n>0)

disp('the system is unstable')

disp('it is no use for getting 动态参数')

else

%调用求取二阶系统阶跃响应动态参数的函数文件 [y,x,t]=step(num,den); figure(1) plot(t,y) xlabel('time-sec')

ylabel('y(t)')

title('step responce') figure(2)

pzmap(num,den)

[mp,ind]=max(y);

%求取时间向量的长度dimt dimt=length(t);

%确定最终的稳定值yss yss=y(dimt); mp=100*(mp-yss)/yss tp=t(ind) % find ts2 i=dimt+1;n=0; while n==0 i=i-1; if i==1 n=1;

elseif y(i)>=1.02 n=1; end

end

t1=t(i);i=dimt+1;n=0; while n==0