clc
clear all
clf
n = 0.1;
m = 500;
s = 20;
w0 = sqrt(s/m);
t_tf = 101;
f = logspace(-3,1,t_tf);
t_w = 2*pi*f;
l_w = length(t_w);
l_w2w0 = 101;
t_w2w0 = logspace(-1,1,t_tf);
t_vu1 = (1-(t_w2w0).^2)+(1i*n.*t_w2w0);
t_vd1 = 1+1i*n.*t_w2w0;
t_V = t_vu1./t_vd1;
L_V = 20*log10(abs(t_V));
phi_V = angle(t_V);
n_fig = 1;
figure(n_fig);
for n_w = 1:l_w
w2w0 = t_w2w0(n_w);
subplot(2,2,1)
semilogx(t_w2w0,L_V)
title('Einfügedämmmaß')
grid on
xticks([0.01 0.1 1 sqrt(2) 10 100])
xticklabels({'0.01','0.1','1','\surd 2','10','100'})
ylim = [min(L_V)*1.2 max(L_V)*1.2];
line(w2w0*[1 1],ylim,'color','r')
xlabel('\omega / \omega_0')
ylabel('Einfügedämmmaß R_E [dB]')
axis([min(t_w2w0) max(t_w2w0) min(L_V)*1.2 max(L_V)*1.2]);
subplot(2,2,3)
semilogx(t_w2w0,phi_V)
grid on
yticks([0 pi/4 pi/2 3/4*pi pi])
yticklabels({'0','\pi/4','\pi/2','3/4 \pi','\pi'})
xticks([0.01 0.1 1 sqrt(2) 10 100])
xticklabels({'0.01','0.1','1','\surd 2','10','100'})
title('Phasenverschiebung')
ylabel('Phasenwinkel in [Rad]')
xlabel('\omega / \omega_0')
ylim = [min(phi_V)-0.5 max(phi_V)+0.5];
line(w2w0*[1 1],ylim,'color','r')
axis([min(t_w2w0) max(t_w2w0) min(phi_V)-0.5 max(phi_V)+0.5]);
subplot(1,2,2)
pF = abs(1/(t_V(n_w))).*exp(1i*phi_V(n_w));
plot([0,real(pF)],[0,imag(pF)],'r','linewidth',1.5);
axis equal
axis ([-5 5 -5 5])
grid on
title('Eingeleitete Kraft in das Fundament')
xlabel('')
ylabel('Kraft')
M(n_w) = getframe(gcf);
end
N = 2;
FPS = 6;
movie(gcf,M,N,FPS)