huettern · May 10, 2016 09:46
diff --git a/polyval skin multi-w.m b/polyval skin multi-w.m
 %% Poly with skin
 clear; clc;

 l=1e-6; c0=1e-9;

 r0=100; f0=1000; alpha=0.8;


 w=[2000 3000];
 w=(2*pi).*w;

 w0 = 2*pi*f0;

 %rskin=r0.*(1+(w./w0 ).^alpha);


 n=[0,0,l,r0.*(1+(w(1)./w0 ).^alpha) ;
   0,0,l,r0.*(1+(w(2)./w0 ).^alpha)];

 for i=1:length(w)
    res1(i)=polyval(n(i,:),j.*w(i));
 end

 res1(1)
 res1(2)
	%% Poly with skin
	clear; clc;

	l=1e-6; c0=1e-9;

	r0=100; f0=1000; alpha=0.8;


	w=[2000 3000];
	w=(2pi).w;

	w0 = 2pif0;

	%rskin=r0.*(1+(w./w0 ).^alpha);


	n=[0,0,l,r0.*(1+(w(1)./w0 ).^alpha) ;
	0,0,l,r0.*(1+(w(2)./w0 ).^alpha)];

	for i=1:length(w)
	res1(i)=polyval(n(i,:),j.*w(i));
	end

	res1(1)
	res1(2)