78 lines
2.7 KiB
Python
78 lines
2.7 KiB
Python
import matplotlib.pyplot as plt
|
||
|
||
def decimal_to_8bit_binary(decimal_number):
|
||
if decimal_number < 0 or decimal_number > 255:
|
||
raise ValueError("十进制数必须在0到255之间。")
|
||
return f'{decimal_number:08b}' # 使用8位表示
|
||
|
||
def read_numbers_from_file(filename):
|
||
with open(filename, 'r') as file:
|
||
binary_numbers = [int(line.strip()) for line in file.readlines() if line.strip().isdigit()]
|
||
return binary_numbers
|
||
|
||
def to_signed_decimal(binary_str):
|
||
# 将二进制字符串转换为有符号十进制数
|
||
if binary_str[0] == '1':
|
||
# 如果最高位是1，表示是负数
|
||
return int(binary_str, 2) - (1 << len(binary_str))
|
||
else:
|
||
return int(binary_str, 2)
|
||
|
||
def main():
|
||
filename = "E:\\02.JobFile\\01.项目资料\\03.风电机组三合一采集器项目开发\\01.机舱采集器\\01.软件\\others\\ADS1278_manualclock\\ADS1278.txt"
|
||
|
||
# 从文件读取十进制数
|
||
decimal_numbers = read_numbers_from_file(filename)
|
||
|
||
# 列表以存储结果
|
||
binary_numbers = []
|
||
|
||
# 每三个十进制数组合成一个24位的二进制数
|
||
for i in range(0, len(decimal_numbers), 3):
|
||
chunk = decimal_numbers[i:i+3]
|
||
if len(chunk) == 3: # 确保是三个数
|
||
# 转换每个十进制数为8位二进制数，并连接
|
||
binary_number = ''.join(decimal_to_8bit_binary(num) for num in chunk)
|
||
signed_decimal = to_signed_decimal(binary_number)
|
||
scaled_value = signed_decimal / 8388607 * 2.5 * 4
|
||
binary_numbers.append(scaled_value)
|
||
else:
|
||
print(f"剩余的数字不足3个, 忽略: {chunk}")
|
||
|
||
# 计算平均数
|
||
average = sum(binary_numbers) / len(binary_numbers)
|
||
|
||
# 计算最大值和最小值
|
||
maximum = max(binary_numbers)
|
||
minimum = min(binary_numbers)
|
||
|
||
# 计算中位数
|
||
sorted_numbers = sorted(binary_numbers)
|
||
n = len(sorted_numbers)
|
||
if n % 2 == 0: # 偶数个数
|
||
median = (sorted_numbers[n // 2 - 1] + sorted_numbers[n // 2]) / 2
|
||
else: # 奇数个数
|
||
median = sorted_numbers[n // 2]
|
||
|
||
# 输出结果
|
||
for value in binary_numbers:
|
||
print(value)
|
||
|
||
print(f"平均值: {average}")
|
||
print(f"最大值: {maximum}")
|
||
print(f"最小值: {minimum}")
|
||
print(f"中位数: {median}")
|
||
|
||
|
||
# 绘制时域图
|
||
plt.figure(figsize=(10, 5))
|
||
plt.plot(binary_numbers, label='Scaled Values', color='b')
|
||
plt.title('Time Domain Signal')
|
||
plt.xlabel('Sample Index')
|
||
plt.ylabel('Amplitude')
|
||
plt.grid(True)
|
||
plt.legend()
|
||
plt.show()
|
||
|
||
if __name__ == "__main__":
|
||
main() |