k=1.4

import math
import numpy as np
import pandas as pd
#Ma1到p之比
def Ma2p(Ma1,k):
    return 2*k*Ma1**2/(k+1)-(k-1)/(k+1)
def Ma2rho(Ma1,k):
    return (k+1)*Ma1**2/(2+(k-1)*Ma1**2)
def Ma2T(Ma1,k):
    return 1/Ma1**2*(2/(k+1))**2*(k*Ma1**2-(k-1)/2)*(1+(k-1)/2*Ma1**2)
def Ma2Ma(Ma1,k):
    return math.sqrt((Ma1**2+2/(k-1))/(2*k*Ma1**2/(k-1)-1))
def Ma2sigma(Ma1,k):
    #Ma到总压比
    rho=Ma2rho(Ma1,k)
    p=Ma2p(Ma1,k)
    return rho**(k/(k-1))/p**(1/(k-1))
N=302
N1=6
dMa1=0.01
Ma1min=1
B=np.zeros([N,N1])
k=1.4
for i in range(0,N):
    Ma1=Ma1min+i*dMa1
    B[i,0]=Ma1
    B[i,1]=Ma2Ma(Ma1,k)
    B[i,2]=Ma2p(Ma1,k)
    B[i,3]=Ma2rho(Ma1,k)
    B[i,4]=Ma2T(Ma1,k)
    B[i,5]=Ma2sigma(Ma1,k)
df=pd.DataFrame(B)
# 将DataFrame写入Excel文件
df.to_excel('output.xlsx', index=False, header=False)

绘制相关图像的代码

import numpy as np
import matplotlib.pyplot as plt
import math
# 设置字体为支持中文的字体，例如：
plt.rcParams['font.sans-serif'] = ['SimHei']  # 使用黑体
plt.rcParams['axes.unicode_minus'] = False  # 用来正常显示负号
def Ma2p(Ma1,k):
    return 2*k*Ma1**2/(k+1)-(k-1)/(k+1)
def Ma2rho(Ma1,k):
    return (k+1)*Ma1**2/(2+(k-1)*Ma1**2)
def Ma2T(Ma1,k):
    return 1/Ma1**2*(2/(k+1))**2*(k*Ma1**2-(k-1)/2)*(1+(k-1)/2*Ma1**2)
def Ma2Ma(Ma1,k):
    return math.sqrt((Ma1**2+2/(k-1))/(2*k*Ma1**2/(k-1)-1))
def Ma2sigma(Ma1,k):
    #Ma到总压比
    rho=Ma2rho(Ma1,k)
    p=Ma2p(Ma1,k)
    return rho**(k/(k-1))/p**(1/(k-1))

Ma1_values = np.linspace(1, 3, 100)  
k=1.4
Ma2=[Ma2Ma(Ma1,k) for Ma1 in Ma1_values]
Ma2T_values = [Ma2T(Ma1,k) for Ma1 in Ma1_values]
p=[Ma2p(Ma1,k) for Ma1 in Ma1_values]
rho=[Ma2rho(Ma1,k) for Ma1 in Ma1_values]
sigma=[Ma2sigma(Ma1,k) for Ma1 in Ma1_values]

# 绘制图像
plt.figure()
plt.plot(Ma1_values,Ma2,label='Ma2')
plt.plot(Ma1_values, Ma2T_values, label='温度比')
plt.plot(Ma1_values,p,label='压力比')
plt.plot(Ma1_values,rho,label='密度比(速度比)')
plt.plot(Ma1_values,sigma,label='总压比')

plt.xlabel('Ma1')
plt.ylabel('正激波参数')
plt.title('正激波参数与Ma1的关系')
plt.legend()
plt.grid(True)
plt.show()

图像示例

相关表格