On Sun, Jun 14, 2015 at 2:49 PM, Thomas D. Dean <address@hidden> wrote:
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
_______________________________________________
Help-octave mailing list
address@hidden
https://lists.gnu.org/mailman/listinfo/help-octave
did you look here for examples