freson ピーク周波数

## 呼び出し手順

    fr=freson(h)

## パラメータ

h

  `syslin` リスト

fr

  ピーク周波数のベクトル (単位: Hz)

## 説明

SISOプラント `h` に関して, ピーク周波数のベクトル(単位: Hz)を返します.

## 例

Continuous time system

    h = syslin('c',-1+%s,(3+2*%s+%s^2)*(50+0.1*%s+%s^2));
    fr = freson(h)
    clf
    bode(h);
    fig=gcf();sca(fig.children(2));
    [phi,db] = phasemag(repfreq(h,fr));
    plot(fr,db(:),'+r')

h=syslin('c',-1+%s,(3+2*%s+%s^2)*(50+0.1*%s+%s^2)); fr=freson(h) clf;bode(h); fig=gcf();sca(fig.children(2)); [phi,db]=phasemag(repfreq(h,fr)); plot(fr,db(:),'+r');

Discrete time system

    h = syslin('c',-1+%s,(3+2*%s+%s^2)*(50+0.1*%s+%s^2));
    hd = dscr(h, 0.05);
    fr = freson(hd)
    clf
    bode(hd);
    fig=gcf();sca(fig.children(2));
    [phi,db] = phasemag(repfreq(hd,fr));
    plot(fr,db(:),'+r');

h=syslin('c',-1+%s,(3+2*%s+%s^2)*(50+0.1*%s+%s^2)); hd=dscr(h,0.05); fr=freson(hd); clf;bode(hd); fig=gcf();sca(fig.children(2)); [phi,db]=phasemag(repfreq(hd,fr)); plot(fr,db(:),'+r');

## 参照

- frep2tf
- zgrid
- h_norm

## 履歴

6.0

handling discrete systems and zpk representation, peak detection improved.