这是一篇测试文章。

1213

from sympy import *
from math import radians, pi
import random
import numpy as np
import matplotlib.pyplot as plt

r = 100
n = 2
position = [i * 40 for i in range(2)]
rou = []
theta = []
i = 0
k = 0
for thetaj_index in range(len(position)):
    thetaj = radians(position[thetaj_index])
    for thetai_index in range(len(position)):
        thetai = radians(position[thetai_index])
        i += 1
        print('外面', i)
        if thetai >= thetaj:  # 此处要求了thetaj要大于thetai
            continue
        for alpha0jd in range(10, 90, 10):
            for alpha0id in range(10, 90, 10):
                # alpha0j = radians(alpha0jd)
                # alpha0i = radians(alpha0id)
                deltatheta = 0.001
                alpha0j= radians(30) + random.uniform(-deltatheta,deltatheta)
                alpha0i= radians(10) + random.uniform(-deltatheta,deltatheta)
                b = sin(alpha0j - thetaj) * sin(alpha0i) - sin(alpha0i - thetai) * sin(alpha0j)
                a = cos(alpha0i - thetai) * sin(alpha0j) - cos(alpha0j - thetaj) * sin(alpha0i)
                if a == 0:
                    continue
                thetak = pi + atan(b / a)
                k += 1
                print('里面', k)
                if thetak <= thetaj and ((thetak - thetai) > pi or (thetak - thetaj) > pi):
                    continue  # 限制thetak必须要大于thetaj
                theta.append(thetak * (180 / pi))
                rouk = r * sin(thetak + alpha0i - thetai) / sin(alpha0i)
                rou.append(rouk)
                print('thetai',(180 / pi) * thetai,'thetaj', (180 / pi) * thetaj, 'alpha0i',(180 / pi) * alpha0i,'alpha0j', (180 / pi) * alpha0j,
                      'rouk',(180 / pi) * rouk,'thetak', (180 / pi) * thetak)
solver = []
solvet = []
theta_circle = np.linspace(0, 2 * np.pi, 1000)
r = 100  # 圆形半径
x = r * np.cos(theta_circle)
y = r * np.sin(theta_circle)
plt.plot(x, y, color='black')
for i in range(len(rou)):
    print(rou[i], theta[i])
    solver.append(rou[i] * cos(theta[i]))
    solvet.append(rou[i] * sin(theta[i]))
plt.scatter(solver, solvet, s=70,color='red',label='计算值')
plt.axis('equal')
roustandard = [100 for _ in range(9)]
thetastandard = [radians(40*k) for k in range(9)]
plt.scatter(0,0,color='blue',s=100)
for i in range(len(roustandard)):
    if i ==0 or i==1:
        plt.scatter(roustandard[i] * cos(thetastandard[i]), roustandard[i] * sin(thetastandard[i]), color='blue', s=100,
                    label='标准值')
        plt.text(roustandard[i] * cos(thetastandard[i])+10, roustandard[i] * sin(thetastandard[i]),i+1,fontsize = 10)
        continue
    elif i == 4:
        plt.scatter(roustandard[i] * cos(thetastandard[i]), roustandard[i] * sin(thetastandard[i]), color='green', s=100,
                    label='标准值')
        plt.text(roustandard[i] * cos(thetastandard[i]) + 10, roustandard[i] * sin(thetastandard[i]), i + 1, fontsize=10)
        continue
    plt.scatter(roustandard[i]*cos(thetastandard[i]),roustandard[i]*sin(thetastandard[i]),color='black',s=100,label='标准值')
    plt.text(roustandard[i] * cos(thetastandard[i]) + 10, roustandard[i] * sin(thetastandard[i]), i + 1, fontsize=10)
plt.text(5,5,0, fontsize=10)
plt.savefig('pic1.jpg',dpi=300)
plt.show()