伝達関数から周波数特性を書きたい

初心者です。伝達関数から周波数振幅特性を書きたいため、調べていて、わからなかったので質問しました。

伝達関数は上図のようなものです。H’1(z)とH'2(z)の積の周波数振幅特性を求めたいです

ほしいのは、以下のような周波数振幅特性です。

αの値を変えるとどのようになるのかを考えるグラフを作りたいです。

以下のようにコードの見本をまねてやってみました。でもできません。matlabはじめて2時間くらいです。教えてくださるとありがたいです。とりあえず数字はめんどいのであとでにゅうりょくします。

この場合どのようなコードを書けばよいかご教授いただけるとありがたいです。

b0 = ;

b1 = [];

b2 = [];

a1 = [];

a2 = [];

b = b0*conv(b1,b2);

a = conv(a1,a2);

[h,w] = freqz(b,a,'whole',2001);

plot(w/pi,20*log10(abs(h)))

ax = gca;

ax.YLim = [-100 20];

ax.XTick = 0:.5:2;

xlabel('Normalized Frequency (\times\pi rad/sample)')

ylabel('Magnitude (dB)')

のようなコードを見つけて