Question

MATLAB WORK PLEASE

An oscillating current in an electric circuit is given by: I = 9e-t sin(2nt) where I is the current in Amperes and t is the time in seconds. Use fzero to determine the two values of t in the range of 0 <t< 1.5 such that I desired = -1.15A. (What does a negative current physically mean?) The outputs of the function are t, and t2, where tı = first time that I desired = -1.15A and t2 = second time that I desired = -1.15A. Function ® C Reset MATLAB Documentation 1 function (t1, t2] = current(I_desired) 2% Inputs: 3 % I_desired: desired current through the circuit (scalar). 4 % 5 % Outputs: 6 % t1: The first time in the range 0 <= t <1.5 where I=I_desired (scalar). 7 % t2: The second time in the range 0 < t < 1.5 where I=I_desired (scalar). 8 % NOTE: t2 should come after ti (i.e., t2 > t1). 9 % 10 % Write your code below: 11 12 13 14

Answer 1

X 1 Editor - F:\Program Files\MATLAB\R2018a\bin\ current_fn.m current_fn.m X + function [ti, t2] = current_fn (I) 2 - F = @(x) 9. *exp(-x). *sin (2. *pi.*x); $current function 3 - interval = [0,1.5]; $0 to 0.7 interval 4 - N=500; start_pts - linspace (interval(1), interval (2),N); 6 7 - for i=1:numel (start_pts)-1 8 5 - 9- try 10 11 - t(i) fzero (( (x)F(x)-I), [start_pts(i), start_pts(i+1)]); $fzero function 12 13 - 14 15 - end end nonzeroIndices = tx=0; Sextract non-zero elements from matrix 't' output = t (nonzeroIndices); $new matrix with extracted non-zeo elements tl-output(1); t2=output (2); 6 17 - 18 19 20 21 - end 22 23 24

Editor - F:\Program Files\MATLAB\R2018\bin\current_fn.m current_fn.mx function [ti, t2] = current_fn (I) F = @(x) 9.*exp(-x).*sin (2. *pi.*x); $current function interval = [0,1.5]; $0 to 0.7 interval 4 - N=500; 5 - start_pts = linspace (interval (1), interval(2),N); 6 for i=1:numel (start_pts)-1 - try I 10 11 12 13 - t(i) = fzero (((x)F(x) -I), [start_pts(i),start_pts (i+1)]); $fzero function end Command Window >> [ti, t2] = current_fn(-1.15) tl = 0.5350 t2 - 0.9466

Call the function by typing in the command window >> [t1,t2] = current_fn(-1.15)

MATLAB CODE-

function [t1,t2] = current_fn(I)
F = @(x) 9.*exp(-x).*sin(2.*pi.*x); %current function
interval = [0,1.5]; %0 to 0.7 interval
N=500;
start_pts = linspace(interval(1),interval(2),N);

for i=1:numel(start_pts)-1
  
try
  
t(i) = fzero((@(x)F(x)-I),[start_pts(i),start_pts(i+1)]); %fzero function

end
end
nonzeroIndices = t~=0; %extract non-zero elements from matrix 't'
output = t(nonzeroIndices); %new matrix with extracted non-zeo elements
t1=output(1);
t2=output(2);
  

end