Skip to content

Instantly share code, notes, and snippets.

@huettern
Created May 10, 2016 09:46
Show Gist options
  • Save huettern/2f8cc4042ef2d05d53fe6d60aa74d6bd to your computer and use it in GitHub Desktop.
Save huettern/2f8cc4042ef2d05d53fe6d60aa74d6bd to your computer and use it in GitHub Desktop.
%% 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)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment