广告公司取名字大全免费查询_推广平台哪个好_百度免费建网站_互联网营销师有什么用

1. 输入值：如图
   
2. 生成的表格:如下

3.代码如下

import matplotlib.pyplot as plt
import numpy as np
####################################################
#给与的值，描绘图片
A_Low = 0.1153607
B_Low = 0.001453239
C_Low = 0.000508077
T_Low = 300.8246A_Mid = 0.1874887
B_Mid = 0.001586056
C_Mid = 0.0002653489
T_Mid = 300.8246A_High = -15.02923
B_High = 0.0003478767
C_High = 0.01089852
T_High = 300.9223##########################################################
def delta(P, A, B, C, T, coffd = 1):fTr50 =np.sqrt(T)fTterm= 100.24 + 0.2881*(T - 273.15)fPov  = Pa = B * fTr50b = C * fTtermc = A - fPov# print("****************P =" + str(P)+"****************"+"\n\r")# print("a =" + str(a)+"\n\r")# print("b =" + str(b)+"\n\r")# print("c =" + str(c)+"\n\r")return (b*b -4.0 *a *c),a,b,cdef solve_formula(P, A, B, C, T, coffd = 1):delta_result,a,b,c = delta(P, A, B, C, T, coffd = coffd)# print("delta_result:" + str(delta_result)+"\n\r")# print("a =" + str(a)+"\n\r")# print("b =" + str(b)+"\n\r")# print("c =" + str(c)+"\n\r")if delta_result<0 :print("0")return#ft = - C * (100.24 + 0.2881*(T - 273.15))# Q1经过测试为负值，所以舍弃，Q2_sqrt为Q-P的关系# Q1_sqrt = (- C * (100.24 + 0.2881*(T - 273.15)) / np.sqrt(coffd) - np.sqrt(delta_result))/(2 * B * np.sqrt(T)/coffd)#Q2_sqrt = (ft / np.sqrt(coffd) + np.sqrt(delta_result))/(2 * B * np.sqrt(T)/coffd)f =(-b + np.sqrt(delta_result))/(2 * a)f = f * coffd#return np.power(Q2_sqrt, 2)# print("f  =" + str(f )+"\n\r")return f * f#由于0到(1.6/100)*5 无解，所以从(1.6/100)*6
P1_Low = np.linspace((1.6/100)*6, 1.6, 100)  #logspace(a, b, n)生成一个数组，数据的第一个元素值为a，最后一个元素为b，n是总采样点数
P1_Mid = np.linspace(1.6, 15, 1000)
P1_High = np.linspace(15, 31, 1000)Q_Low = []  #一种有序、可变的数据结构，可以存储不同类型的元素
Q_Mid = []
Q_High = []for p in P1_Low:temp_low = solve_formula(p, A_Low, B_Low, C_Low, T_Low)#print("temp_low =" + str(temp_low)+"\n\r")Q_Low.append(temp_low)
for p in P1_Mid:temp_Mid= solve_formula(p, A_Mid, B_Mid, C_Mid, T_Mid)#print("temp_Mid =" + str(temp_Mid)+"\n\r")Q_Mid.append(temp_Mid)#Q_Mid.append(solve_formula(p, A_Mid, B_Mid, C_Mid, T_Mid))
for p in P1_High:temp_High= solve_formula(p, A_High, B_High, C_High, T_High)#print("temp_High =" + str(temp_High)+"\n\r")Q_High.append(temp_High)# Q_High.append(solve_formula(p, A_High, B_High, C_High, T_High))#P1 = float(13)  #输入P1的值
P1 = input("请输入:")
P1 = float(P1)
print("你输入的内容是:",P1)
if (P1>=(1.6/100)*6) & (P1 < 1.6):print(solve_formula(P1, A_Low, B_Low, C_Low, T_Low))
elif (P1 >=1.6) & (P1 < 15):print(solve_formula(P1, A_Mid, B_Mid, C_Mid, T_Mid))
elif (P1>=15) & (P1<=31) :print(solve_formula(P1, A_High, B_High, C_High, T_High))
else:print("输入P1值无效!由于0到(1.6/100)*5 无解,所以P1的范围为" + str((1.6/100)*6) + "<=P1<=" + str(31))#显示图表
plt.figure(0) 
plt.plot(P1_Low, Q_Low, '-', label='Low')
plt.plot(P1_Mid, Q_Mid, '-', label='Mid')
plt.plot(P1_High, Q_High, '-', label='High')#画图例
plt.legend()
#显示中文
plt.rcParams['font.sans-serif'] = ['SimHei']
#指定坐标轴名称
plt.xlabel('P1')
plt.ylabel('Qm')
#显示
plt.show()