1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29
|
clc;
clear all;
close all;
%---de=1/2, und do=-1/4---
de=0.5;
do=-0.25;
%--im Fall einer rechtseitige Diskontinuität, d.h einer Transmission---
%--U11/Uo= ½[(1-de²)*(1+de^2+?e^4+de^6+...)]+½[(1-do²)*(1+do^2+do^4+do^6+...)]---
para_1=1/2*(1-de^2);
para_2=1/2*(1-do^2);
y=zeros(5,1);% Matrix null mit 5 Zeile und 1 Spalte
t=zeros(5,1);
t(1:5)=1:2:9; %---hier t/T=1; t/T=3;t/T=5;t/T=7;t/T=9
%---k=0---- wenn t/T=1---kommmt die erste Wert :1 von der Gleichung---
y(1,1)=(para_1*1)+(para_2*1);
%---k=1--wenn t/T=3, wird die zweite Wert von der Gleichung hinzugefügt
y(2,1)=(para_1)*(1+de^2)+((para_2)*(1+do^2));
%---k=2--wenn t/T=5, der dritte Wert von der Gleichung kommt dazu
y(3,1)=(para_1)*(1+de^2+de^4)+((para_2)*(1+do^2+do^4));
%---k=3--für t/T=7, der vierte Wert von der Gleichung wird dazu addiert
y(4,1)=(para_1)*(1+de^2+de^4+de^6)+((para_2)*(1+do^2+do^4+do^6));
%---k=4--für t/T=9, der fünfte Wert wird dazu addiert
y(5,1)=(para_1)*(1+de^2+de^4+de^6+de^8)+((para_2)*(1+do^2+do^4+do^6+do^8));
stairs([1 ; t(:,1)],[0 ; y(:,1)]);
ylim([0 1]);
xlim([0 9]);
xlabel('t/T');
ylabel('U11/U0'); |
Partager