matlab编程怎么收费,求matlab大佬帮忙编程写一个小公式

CODE:

function zzz35A1

clc

clear all

close all

d = 2;

k1 = 0.8;

k2 = 0.8;

r1 = 0.1;

r2 = 0.1;

D1=1000;

D2=500;

L1=pi*D1;

L2=pi*D2;

n1=1.45;

n2=1;

Delte1 = 0.0001;   % 步长

lambdaYY = 1.55:Delte1:1.56;  % lambda取值   lambda =1.55:0.00001:1.56

[EC2, EC3, EC4, EC5, EC6, EC7, EC8] = arrayfun(@(lambda) EE28(lambda),lambdaYY);  % lambda对应的E值

% EEE = [EC2, EC3, EC4, EC5, EC6, EC7, EC8]';

EEE = [EC2; EC3; EC4; EC5; EC6; EC7; EC8]';

EEE(1,:)'   % lambda(1)对应的E值

EEE(2,:)'   % lambda(2)对应的E值

EC2'   % E2值

EC3'   % E3值

function [e2, e3, e4, e5, e6, e7, e8] = EE28(lambda)

Gamma = 2.246/d*exp( (n1^2+n2^2)/8/n2^2 - (n1^2+n2^2)/(n1^2-n2^2)/n2^2*lambda^2/(pi*d)^2 );

Beta = 2*pi/lambda + Gamma*Gamma*lambda/4/pi;

c1 = (1-r1)^(1/2);

c2 = (1-r2)^(1/2);

c3 = (1-k1)^(1/2);

c5 = (1-k2)^(1/2);

c4 = (k1)^(1/2);

c6 = (k2)^(1/2);

kk1 = exp(Beta*L1*1j/2);

kk2 = exp(Beta*L2*1j);

e8 = (1j*c1*c2*c3*kk1*(1j*c6 + c2*c5^2*kk2 + c2*c6^2*kk2))/(1j + c2*c6*kk2 + c1*c2*c4*kk1^2*(1j*c6 + c2*c5^2*kk2 + c2*c6^2*kk2));

e3 = (c1*c3*(1j + c2*c6*kk2))/(1j + c2*c6*kk2 + c1*c2*c4*kk1^2*(1j*c6 + c2*c5^2*kk2 + c2*c6^2*kk2));

e7 = (1j*c1*c2*c3*c5*kk1)/(c1*c2*c4*kk1^2*(c2*c5^2*kk2 + c2*c6^2*kk2 + 1j*c6) + c2*c6*kk2 + 1j);

e4 = (1j*c1*(c1*c2*c3^2*kk1^2*(c2*c5^2*kk2 + c2*c6^2*kk2 + 1j*c6) + c1*c2*c4^2*kk1^2*(c2*c5^2*kk2 + c2*c6^2*kk2 + 1j*c6) + c4*(c2*c6*kk2 + 1j)))/(c1*c2*c4*kk1^2*(c2*c5^2*kk2 + c2*c6^2*kk2 + 1j*c6) + c2*c6*kk2 + 1j);

e2 = kk1*e8;

e5 = kk1*e3;

e6 = kk2*e7;

end

end