bodeplot load iddata1 sys1 = n4sid(z1, 2) % discrete-time IDSS model of order 2 sys2 = n4sid(z1, 6) w = linspace(8,10*pi,256); h = bodeplot(sys1,sys2,w);
那你的把实验数据说明啊,如果你有传递函数,那就直接用bode()命令即可,如果不是,那楼主你得要细说清楚了
两个图的坐标系不一样啊.虽然横坐标都是频率omiga,但第一个是实数坐标,第二个是对数坐标.
1、a为矩阵norm(a)/norm(a,2),返回的是矩阵a的二范数,(二范数j就是矩阵a的2范数就是 a的转置矩阵乘以a特征根 最大值的开根号)norm(a,1),返回矩阵的1泛数,就是最大一列的和;norm(a,'inf') 返回矩阵的无穷泛数,也就是最大一行的和;norm(a,'fro') 返回矩阵的frobenius范数.2、a为向量norm(a,p) 返回向量a的p范数.即返回 sum(abs(a).^p)^(1/p),对任意 1 评论0 0 0
直接bode(sys)就可以了啊.比如:num=[1 3 5];den=[2 5 7 8 3]; sys=tf(num,den); bode(sys);grid on;
一般画图用plot(x,y)x轴用对数坐标画图的用semilogx(x,y)
w=logspace(-1,1,32); %w范围和点数ngw=10./((i*w).^2+2*(i*w)+10);% 计算模值mag=abs(gw);l=20*log(mag); % 模取对数semilogx(w,l); % 半对数作图grid on% 5 画网格线
^simple(S),采用多种方法化简多项式,包括simplify,expand、factor、combine、radsimp、convert、collect采用radsimp结果就是x估计程序问题>> syms x>> s=x^(2)^(1/2);>> simple(s)simplify:csgn(x)*xradsimp:xcombine(trig):(x^2)^(1/2)factor:(x^2)
在matlab帮助文档中 help bode有实例
利用matlab中的函数bode()或margin(),用margin()可以直接显示出截止频率、相角穿越频率、幅值裕度和相角裕度的结果,用bode()可以由画出的伯德图读出上述结果.