is there way to convert this python code to matlab code? how can i convert python code to matlab? this is the code that I want to convert:
4 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
import numpy as np
from matplotlib import pyplot as plt
from scipy.misc import derivative as drv
from scipy.integrate import quad
#import sympy as sp
C1 = 3.74177e8; C2 = 1.43878e4; sig = 5.67e-8
T_static = 229.58 #changes with altitude, this is for altitude of 9 km
Tval = open("Tawvalue.txt")
Tvalue = Tval.read()
Tvalue = Tvalue.strip()
Tvalue = Tvalue.split('\n')
hval = open("hvalue.txt")
hvalue = hval.read()
hvalue = hvalue.strip()
hvalue = hvalue.split('\n')
T = []; H = []
for i in range(len(value)):
a = float(value[i])
b = float(hvalue[i])
H.append(b)
T.append(a)
def wall(T_w):
wallt = H[i]*T[i] -T_w*H[i] - e*sig*(T_w**4) + e*sig*(T_static**4)
return wallt
def power(lam):
return (e*C1)/((lam**5)*(np.e**(C2/(lam*c)) - 1))
c = T_static; e = 1; rad = quad(power,8,12)
Temp_value = []; optimum_emsvty = []
for i in range(len(H)):
e0=0.0; e1=1.0; delt = 1
while (abs(delt)>0.000001):
e = (e0 + e1)/2.0
c = 500
for j in range(5):
c = c - wall(c)/drv(wall,c)
powr = quad(power,8,12)
delt = powr[0] - rad[0]
if delt>0:
e1 = e
elif delt<0:
e0 = e
c = float(c)
optimum_emsvty.append(e)
Temp_value.append(c)
file = open("Optimum_Emissivity.txt","w")
for item in optimum_emsvty:
file.write("%s\n" %item)
file = open("Final_Wall_Temperature.txt","w")
for item in Temp_value:
file.write("%f\n" %item)
0 Kommentare
Antworten (1)
Steven Lord
am 5 Jun. 2024
Do you need to convert the code or do you just need to run the code? If the latter, see this section of the documentation.
0 Kommentare
Siehe auch
Kategorien
Mehr zu Call Python from MATLAB finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!