حل مشکل کد در پایتون /محاسباتی

شرح پروژه

وقتی ضریب شکست که وابسته به طول موج هست وقتی داخل ماتریس قرار میگیره ارور میده

کد رو گذاشتم کپی کنید 

#s wave import numpy as np import math as m import matplotlib.pyplot as plt import cmath    N = 12 n0 = 1.0003#AIR npr = 1.5151#PRISM n1 = 2.20760#ZrO2 n2 = 1.45704#SiO2 ns = 1.0003#AIR l = 635  #nm  a = l / (4*n1) b = l / (4*n2) dp = l / (4*npr) ds = l / (4*ns) dd = 0.65  E0 = 4.14 wp = 1 aj = [0, 0.4, 1.75, 42.5, 5, 160] bj = [0, 0.059, 0.12, 0.37, 1.2, 1.3] wj = [0, 1.89, 2.03, 2.877, 2.275, 3.745] ew = [] ewr = [] w = np.linspace(300, 900, 1000) w1 = w[::-1] Ew = E0 + ((aj[1]*(wp2))/((wj[1]2)-((w1/(2.95e2))2)-(1j*(w1/(2.95e2))*bj[1])))+ ((aj[2]*(wp2))/((wj[2]2)-((w1/(2.95e2))2)-(1j*(w1/(2.95e2))*bj[2])))+ ((aj[3]*(wp2))/((wj[3]2)-((w1/(2.95e2))2)-(1j*(w1/(2.95e2))*bj[3]))) Ew1 = Ew + ((aj[4]*(wp2))/((wj[4]2)-((w1/(2.95e2))2)-(1j*(w1/(2.95e2))*bj[4])))+ ((aj[5]*(wp2))/((wj[5]2)-((w1/(2.95e2))2)-(1j*(w1/(2.95e2))*bj[5]))) e1 = np.real(Ew1) e2 = np.imag(Ew1) nd = ((np.sqrt((e12)+(e22)) + e1)/2)(1/2)  theta0 = (0.0*np.pi)/180  theta2 = np.arcsin((npr/n2)*(np.sin(theta0))) theta1 = np.arcsin((n2/n1)*(np.sin(theta2))) thetad = np.arcsin((n1/nd)*(np.sin(theta1)))#rad thetas = np.arcsin((n1/ns)*(np.sin(theta1)))  kp = (npr)*(2*np.pi)*(np.cos(theta0))*(1/w) k1 = (n1)*(2*np.pi)*(np.cos(theta1))*(1/w) k2 = (n2)*(2*np.pi)*(np.cos(theta2))*(1/w) kd = (nd)*(2*np.pi)*(np.cos(thetad))*(1/w) ks = (ns)*(2*np.pi)*(np.cos(thetas))*(1/w)  Dp = np.matrix([[1, 1], [(npr*np.cos(theta0)), (-npr*np.cos(theta0))]]) D1 = np.matrix([[1, 1], [(n1*np.cos(theta1)), (-n1*np.cos(theta1))]]) D2 = np.matrix([[1, 1], [(n2*np.cos(theta2)), (-n2*np.cos(theta2))]]) Dd = np.matrix([[1, 1], [(nd*np.cos(thetad)), (-nd*np.cos(thetad))]]) Ds = np.matrix([[1, 1], [(ns*np.cos(thetas)), (-ns*np.cos(thetas))]])  Dp_1 = (np.matrix([[(-npr*np.cos(theta0)), -1], [(-npr*np.cos(theta0)), 1]])) * (1/(-2*npr*np.cos(theta0))) D1_1 = (np.matrix([[(-n1*np.cos(theta1)), -1], [(-n1*np.cos(theta1)), 1]])) * (1/(-2*n1*np.cos(theta1))) D2_1 = (np.matrix([[(-n2*np.cos(theta2)), -1], [(-n2*np.cos(theta2)), 1]])) * (1/(-2*n2*np.cos(theta2))) Dd_1 = (np.matrix([[(-nd*np.cos(thetad)), -1], [(-nd*np.cos(thetad)), 1]])) * (1/(-2*nd*np.cos(thetad)))  Pp = np.matrix([[np.exp(1j*(kp*dp)), 0], [0, np.exp(-1j*(kp*dp))]]) P1 = np.matrix([[np.exp(1j*(k1*a)), 0], [0, np.exp(-1j*(k1*a))]]) P2 = np.matrix([[np.exp(1j*(k2*b)), 0], [0, np.exp(-1j*(k2*b))]]) Pd = np.matrix([[np.exp(1j*(kd*dd)), 0], [0, np.exp(-1j*(kd*dd))]])  DPd = np.dot(Dd, Pd) DPDd = np.dot(DPd, Dd_1) D1P1 = np.dot(D1, P1) D_12 = np.dot(D1_1, D2) P2D_2 = np.dot(P2, D2_1) M1 = np.dot(D1P1, D_12) M2 = np.dot(M1, P2D_2) M3 = (M2)(int(N/2)) M4 = np.dot(Dp_1, M3) M5 = np.dot(M4, DPDd) M = np.dot(M5, Ds) rn = (M[1, 0])/(M[0, 0]) R = (np.abs(rn))2 fig = plt.figure(figsize=(9,5)) ax  = fig.add_subplot(111) ax.set_ylabel(' Reflectance for TE waves(with defect) ') ax.set_xlabel(' \u03BB(nm) ') plt.style.use('ggplot') ax.plot(w, R, 'b')

مهارت ها و تخصص های مورد نیاز

درباره کارفرما

Sam19998

عضویت شش سال پیش

30 پروژه ثبت شده ،

0 پروژه در حال انجام ،

0 پروژه آماده دریافت پیشنهاد ،

نرخ پذیرش پیشنهاد 43%