Commit 8e0f5322 authored by johankor's avatar johankor
Browse files

conflict

parents d4912053 9624d80e
Pipeline #59 failed with stage
Change Component Designator: Old Designator=R? New Designator=R6
% a = Z1
% b = Z2
% c = Z3
% d = Z4
% e = I1
% f = U1
syms a b c d e f
eqn = e == (f+b*((b*e)/(b-c-d)))/(a+b)
sole = solve(eqn, e)
C1 = 1e-7;
C2 = 1e-11;
L1 = 1e-5;
L2 = 1e-1;
M = 1e-6;
R1 = 0;
R2 = 0;
%G1 = 1e-2;
G1 = 0;
% Time domain parameters
fs = 4e6; % Sampling frequency
dt = 1/fs; % Time resolution
T = 1; % Signal duration
t = 0:dt:T-dt; % Total duration
N = length(t); % Number of time samples
f0=1/(2*pi*(sqrt(L1*C1)))
f0_s=1/(2*pi*(sqrt(L2*C2)))
x2=square(2*pi*f0*t);
a = ((C1*C2*G1*L1*L2)-2*(C1*C2*G1*L1*M)+(C1*C2*G1*M^2));
b = ((C1*C2*G1*L1*R2)+(C1*C2*G1*L2*R1)-2*(C1*C2*G1*M*R1)+(C1*C2*L1));
c = ((C1*C2*G1*R1*R2)+(C1*C2*R1)+(C1*G1*L1)+(C2*G1*L2)-2*(C2*G1*M));
d = ((C1*G1*R1)+(C2*G1*R2)+C2);
e = (G1);
f = (-1)*(C1*C2*M);
g = (-1)*(C1*G1*M);
H = tf([f g 0 0],[a b c d e]);
figure;
lsim(H,x2,t);
axis([0 1e-4 -10 10]);
xlabel('Seconds');
ylabel('Amplitude');
pbaspect([2 1 1]);
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment