清空筛选条件
清空筛选条件

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 次查看（过去 30 天）
Aditi
Aditi 2024-6-5
回答： Steven Lord 2024-6-5
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)

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

回答（1 个）

Steven Lord
Steven Lord 2024-6-5
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.

类别

在 Help Center 和 File Exchange 中查找有关 Call Python from MATLAB 的更多信息

标签

产品


版本

R2023b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by