how to calculate a difference of function with any parameter

Question

shiv gaur 2022-1-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1624945-how-to-calculate-a-difference-of-function-with-any-parameter

评论： shiv gaur 2022-1-8

采纳的回答： Torsten

function y=f(x)

for i=10

y=i*sin(x)

end

calculate the change in y/change in x vs i

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Torsten 2022-1-8

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1624945-how-to-calculate-a-difference-of-function-with-any-parameter#answer_870730

编辑：Torsten 2022-1-8

在 MATLAB Online 中打开

function main
  h = 1e-6;
  xleft = -pi;
  xright = pi;
  n = 200;
  dx = (xright-xleft)/n;
  x = xleft:dx:xright;
  change_in_y_by_change_in_x = (f(x+h)-f(x))/h;
  plot(x,change_in_y_by_change_in_x)
  hold on
  plot(x,cos(x))    % for comparison, plot (sin(x))' = cos(x)
end
function y = f(x)
  y = sin(x)
end

7 个评论
显示 5更早的评论隐藏 5更早的评论

shiv gaur 2022-1-8

Dr torsten this the program you help me for solving equation I am little bit stuck in program pl help

program is below

function shiv7

%K0= (2*pi/1530)*1e6:(2*pi/10)*1e6:(2*pi/1630)*1e6;

%K0= linspace((2*pi/1530)*1e6,(2*pi/1630)*1e6,5);

%D1=1e-9:1e-9:3e-7;

D1= linspace(1e-9,3e-7,15);

for j=1:numel(D1)

da = D1(j);

p0 = 1;

p1 = 1.5;

p2 = 2;

TOL = 10^-8;

N0 = 100; format long

h1 = p1 - p0;

h2 = p2 - p1;

DELTA1 = (f(p1,da) - f(p0,da ))/h1;

DELTA2 = (f(p2,da ) - f(p1,da ))/h2;

d = (DELTA2 - DELTA1)/(h2 + h1);

i=3;

while i <= N0

b = DELTA2 + h2*d;

D = (b^2 - 4*f(p2,da )*d)^(1/2);

if abs(b-D) < abs(b+D)

E = b + D;

else

E = b - D;

end

h = -2*f(p2,da)/E;

p = p2 + h;

if abs(h) < TOL

disp(p)

break

end

p0 = p1;

p1 = p2;

p2 = p;

h1 = p1 - p0;

h2 = p2 - p1;

DELTA1 = (f(p1,da ) - f(p0,da ))/h1;

DELTA2 = (f(p2,da ) - f(p1,da ))/h2;

d = (DELTA2 - DELTA1)/(h2 + h1);

i=i+1;

end

%pv(j)=p;

P(j)=real(p)

end

plot(D1,P)

end

function y=f(x,da )

l=6330*(10)^(-10);

k0=2*pi/l;

nc=1.33;

ec=nc^2;

na=1.63;

ea=na^2;

nf=1.59;

ef=nf^2;

nm=0.064+1i*4;

em=nm^2;

ns=1.495;

es=ns^2;

df=0.5e-6;

dm=0.03e-6;

kc=k0*sqrt(x^2-ec);

ka=k0*sqrt(ea-x^2);

kf=k0*sqrt(ef-x^2);

km=k0*sqrt(em-x^2);

ks=k0*sqrt(x^2-es);

rfm=(kf-km)/(kf+km);

rms=(km-ks)/(ks+km);

rfa=(kf-ka)/(kf+ka);

rac=(ka-kc)/(ka+kc);

a=(1-rfm)/(1+rfm);

b=(1-rms*exp(2*1i*dm*km))/(1+rms*exp(2*1i*dm*km));

c=(1-rfa)/(1+rfa);

f=(1-rac*exp(2*1i*da*ka))/(1+rac*exp(2*1i*da*ka));

pfms=2*atan(1i*a*b);

pfac=2*atan(1i*c*f);

y=2*kf*df-pfms-pfac;

end

we have to same find out the delta(y)/delta(na) vs da plot thanks

Torsten 2022-1-8

在 MATLAB Online 中打开

Please include the function f in which you calculate "N" dependent on "da" callable as

N = f(da)

shiv gaur 2022-1-8

Dr torsten N is the root ie p function is dependent on da also the change in the value of p for different value of da for some value of na that is given this is running program so Dr torsten pl help the by varying the value of da the value of changing p ie N so we have to calculate the DP/Dda vs da plot

this is running script so minor mistake in program so pl help for three value of na 1.45 1.49 1.50

DP/Dda vs da thanks fig like att .file

请先登录，再进行评论。

how to calculate a difference of function with any parameter

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

7 个评论
显示 5更早的评论隐藏 5更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

how to calculate a difference of function with any parameter

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

7 个评论 显示 5更早的评论隐藏 5更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

7 个评论
显示 5更早的评论隐藏 5更早的评论