How can I correctly implement this function (attenuation dependent on the frequency)
in Matlab, where and.
I expect a plot with one graph like this
My code so far is:
%Schirmdämpfung bei longitudinalen magnetischem Feld
clear *;clc;r_0 = 6.2E-4; % [m] Innenradius Schirm
d = 2E-5; % [m] Dicke des Schirms
% Frequenz und Winkelgeschwindigkeit
f1 = 100E6; % [Hz] Frequenz zum Starten
df = 100E6; % [Hz] Schrittweite der Frequenz
fmax = 10E9; % [Hz] letzte Frequenz
f = f1:df:fmax; % [Hz] Frequenz als Vektor ausgeben
omega = 2*pi*f; % [rad/s] Winkelgeschwindigkeit
mu0 = 1.256637061E-6; % [N/A^2] mag. Feldkonstante
mu_r = 1; % realtive mag. Permeabilität
mu = mu0*mu_r; % [N/A^2] effektive mag. Permeabilität
% elektrische Leitfähigkeit
sigma_Al = 37E6; % [S/m] Leitfähigkeit von Aluminium
sigma_Co = 58E6; % [S/m] Leitfähigkeit von Kupfer
% Wirbelstromkonstante
k_w_Al = sqrt(1i*omega*mu*sigma_Al); % Wirbelstromkonstante für Aluminium (3.24)
k_w_Co = sqrt(1i*omega*mu*sigma_Co); % Wirbelstromkonstante für Kupfer (3.24)
% Eindringtiefe
delta_Al = sqrt(1./(pi*f*mu*sigma_Al)); % Eindringtiefe Aluminium
delta_Co = sqrt(1./(pi*f*mu*sigma_Co)); % Eindringtiefe Kupfer
% Abkürzung K
K_Al = k_w_Al*r_0/mu_r; % Abkürzung nach (3.68) für Aluminium
K_Co = k_w_Co*r_0/mu_r; % Abkürzung nach (3.68) für Kupfer
% Dämpfung
a_m_Al = 20.*log10(cosh(k_w_Al.*d)+0.5.*K_Al.*sinh(k_w_Al.*d)); % logarithmische Dämpfung für Aluminium
%a_m_Al_O = cosh(k_w_Al.*d)+0.5.*K_Al.*sinh(k_w_Al.*d);
%a_m_Co = 20*log10(cosh(k_w_Co*d)+0.5*K_Co*sinh(k_w_Co*d)); % logarithmische Dämpfung für Kupfer
%plot(f,a_m_Al)
Best Answer