Re: Forming Cell of systems

[Thomas D. Dean]

I am getting closer, I think.

A0 = 10;
b = 1 / A0;
R1 = 10000;
R2 = R1 * (1/b - 1);
K = R1/(R1+R2);
C = [1:.2:3]*1e-12;
a0 = 1e5;
w1 = 1e4;
w2 = 1e6;
s = tf('s');
a = a0/(1+s/w1)/(1+s/w2);
A = feedback(a,b);

#######################################
## still need a cellfun to replace this loop
b_array=cell(size(C'));
for idx=1:length(C)
  b_array{idx,1}=tf([K*R2*C(idx) K],[K*R2*C(idx) 1]);
endfor;
##
#######################################

A_array=cellfun(@feedback,{a},b_array, 'uniformoutput', false);

function [x]=sys_prod(S1,S2)
  x=S1*S2;
endfunction;

L_array=cellfun(@sys_prod,{a},b_array, 'uniformoutput', false);
[Gm,Pm,Wcg,Wcp] = cellfun(@margin,L_array);

figure 1
step(A,'r',A_array{:});
figure 2
bode(A,A_array{:})
figure 3
plot(C,Pm)

Tom Dean