鱼C论坛

 找回密码
 立即注册
查看: 2976|回复: 3

怎么利用python把极半径和极角画出来(求助)

[复制链接]
发表于 2014-9-24 23:15:56 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能^_^

您需要 登录 才可以下载或查看,没有账号?立即注册

x
我是python3.4,我将极半径和极角存为了2个数组,现在我想将极半径和极角通过描点的方式画出一个图形,请问怎么搞啊?

各位大神帮帮忙,没有思路中。。。

将此极半径和极角画出图像并显示。。。thank you~~

我的源码和文件如下:

import math as m
f = open('c:\\circle.txt')
data1 = []
data2 = []
temp = []
data3 = []
value1 = []
value2 = []
temp1 = ''
temp2 = ''
temp3 = ''
#字符串切片,将XYZ储存到data1 data2 data3 并#去掉data3 中的\n
for i in f:
    (temp1,temp2,temp3) = i.split(' ')
    data1.append(temp1)
    data2.append(temp2)
    temp.append(temp3)

#去除\n
for each in temp:
    each = each[:16]
    data3.append(each)

#确定测量点数
len1 = len(data1)
len2 = len(data2)
#转换为极坐标,并判断处于哪个象限,math库中math.atan貌似只能计算(-90)----(90),这里应该可以优化????
#将极半径存到value1
#将极角存到value2

if len1 == len2 :
    for j in range(len1):
        r = m.sqrt(float(data1[j])**2 + float(data2[j])**2)
        value1.append(r)
        if float(data1[j]) > 0 and float(data2[j]) > 0:
            a = m.degrees(m.atan(float(data2[j])/float(data1[j])))
        elif float(data1[j]) < 0 and float(data2[j]) > 0 :
            a = 180 + m.degrees(m.atan(float(data2[j])/float(data1[j])))
        elif float(data1[j]) < 0 and float(data2[j]) < 0 :
            a = 180 + m.degrees(m.atan(float(data2[j])/float(data1[j])))
        else:
            a = 360 + m.degrees(m.atan(float(data2[j])/float(data1[j])))
        value2.append(a)
else:
    print('x number is diff from y number,pls check! ')
#验证
print(list(value1))
print(len(value1))
print(list(value2))
print(len(value2))

f.close()
input()

circle.zip

8.97 KB, 下载次数: 5

小甲鱼最新课程 -> https://ilovefishc.com
回复

使用道具 举报

发表于 2014-9-25 21:09:23 | 显示全部楼层
:titter:
小甲鱼最新课程 -> https://ilovefishc.com
回复

使用道具 举报

 楼主| 发表于 2014-9-25 22:45:07 | 显示全部楼层

笑有用吗???????你是哪个哦
小甲鱼最新课程 -> https://ilovefishc.com
回复 支持 反对

使用道具 举报

头像被屏蔽
发表于 2014-11-1 00:50:32 | 显示全部楼层
提示: 作者被禁止或删除 内容自动屏蔽
小甲鱼最新课程 -> https://ilovefishc.com
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

小黑屋|手机版|Archiver|鱼C工作室 ( 粤ICP备18085999号-1 | 粤公网安备 44051102000585号)

GMT+8, 2025-11-16 07:30

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表