- android - 多次调用 OnPrimaryClipChangedListener
- android - 无法更新 RecyclerView 中的 TextView 字段
- android.database.CursorIndexOutOfBoundsException : Index 0 requested, 光标大小为 0
- android - 使用 AppCompat 时,我们是否需要明确指定其 UI 组件(Spinner、EditText)颜色
我有一个 101x101
矩阵,我想以图形方式显示它。到目前为止,我使用了 matplotlib.pyplot
中的函数 matshow
,如下所示:
import numpy as np
import random
import matplotlib.pyplot as plt
A = np.zeros((101,101))
# do stuff to randomly change some values of A
plt.ion()
plt.matshow(A,cmap='PuBuGn')
plt.colorbar()
plt.show()
您应该将其视为物种之间的相互作用矩阵,正如您所见,只有三个物种存在强烈相互作用。这就是为什么考虑像我在 paper 中找到的下图这样的可视化看起来是合适的。 ,但我不知道这是否以及如何在 Python 中实现:
最佳答案
这应该可以解决问题:
import numpy as np
import matplotlib.pyplot as plt
def old_graph(A):
plt.matshow(A,cmap='PuBuGn')
plt.colorbar()
plt.title(r"abs$\left(\left[\mathbf{A}_{ij} \right ] \right )$ ; SIS=%d"%(sis,), va='bottom')
plt.show()
def new_graph(A, sis_list=np.zeros(0,int), symmetric=True, fig=None, pos=111):
#create and edit figure:
if fig is None:
fig = plt.figure()
ax = fig.add_subplot(pos, projection='polar')
ax.set_rgrids([1],[' '])
ax.set_rmax(1)
ax.set_thetagrids([])
ax.set_title(r"abs$\left(\left[\mathbf{A}_{ij} \right ] \right )$ ; SIS=%d"%(sis,), va='bottom')
colormap = plt.get_cmap('PuBuGn')
# make each species an angle value:
n_species = A.shape[0]
angles = np.linspace(0, 2*np.pi, n_species+1)
# the radius will always be r_max, and each line
# will always unite two points:
r = np.ones((2,))
# prepare list of lines to sort:
unordered_pairs_and_values = []
for index_line in xrange(n_species):
for index_column in xrange(index_line):
if symmetric:
value = A[index_line,index_column]
else: # not symmetric
value= abs(A[index_line,index_column]- A[index_column,index_line])
unordered_pairs_and_values.append([[angles[index_line],angles[index_column]], value])
# sort the lines (otherwise white lines would cover the 'important' ones):
ordered_pairs_and_values = sorted(unordered_pairs_and_values, key=lambda pair: pair[1])
# get the maximum value for scaling:
I_max = ordered_pairs_and_values[-1][1]
# plot every line in order:
for pair in ordered_pairs_and_values:
ax.plot(pair[0], r, color=colormap(pair[1]/I_max), linewidth=2, alpha=0.8)
# don't know how to add the colorbar:
#fig.colorbar(orientation='horizontal')
# mark the angles:
ax.plot(angles, np.ones(angles.shape), 'ko')
# mark the important angles (comment if you don't know which ones are these):
ax.plot(angles[sis_list], np.ones(sis_list.shape), 'ro')
fig.show()
if __name__ == '__main__':
n_species = 51
sis = 3 # strongly interacting species
sis_factor = 4.
A = np.zeros((n_species,n_species))
# do stuff to randomly change some values of A:
for index_line in xrange(n_species):
for index_column in xrange(index_line+1):
A[index_line,index_column] = np.random.random()
A[index_column,index_line] = A[index_line,index_column]
sis_list = np.random.randint(0,n_species,sis)
for species in sis_list:
A[species,:] *= sis_factor
A[:,species] *= sis_factor
for species2 in sis_list: # correct crossings
A[species,species2] /= sis_factor
# stuff to randomly change some values of A done
old_graph(A=A)
new_graph(A=A, sis_list=sis_list, symmetric=True)
#ax.set_rticks([])
)关于python - matshow() 的替代方案,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/40103855/
我应该编写一个函数来打印一组给定的三个数字中两个较大数字的平方和。 我对这种情况的处理相当笨拙。我没有编写返回一组 3 中最大的两个数字的函数,而是编写了函数,以便表达式减少到两个所需的数字。 # S
如果有人可以提供帮助,我将不胜感激。我一直在敲我的头一天试图让这个工作。我已经在互联网上搜索并重新阅读了手册,但我就是不明白。 guile << __EOF__ ( define heading-li
目前我正在处理一个方案问题,其中我们正在使用方案列表表示一个图。我们使用的第一个变体是表示为 的边列表图 '((x y) (y z) (x z)) 我们正在使用的图的第二个变体被称为 x 图,表示为
我正在尝试创建一个函数,该函数将两个函数作为参数并执行它们。 我尝试使用 cond ,但它只执行 action1 . (define seq-action (lambda (action1 act
我提前为我的原始英语道歉;我会尽量避免语法错误等。 两周前,我决定更新我对 Scheme(及其启示)的知识,同时实现我在手上获得的一些数学 Material ,特别是我注册的自动机理论和计算类(cla
Scheme中有没有函数支持分数的“div”操作? 意思是 - 11 格 2.75 = 4。 最佳答案 我认为你的问题的答案是:没有,但你可以定义它: #lang racket (define (di
我在scheme中实现合并排序,我必须通过定义两个辅助方法来实现:merge和split。 Merge 需要两个列表(已经按递增顺序)并将它们合并在一起。我这样做了如下: (define merge
尝试从终端加载方案文件。我创建了一个名为 test.scm 的文件,其中包含以下代码: (define (square x) (* x x)) (define (sum-of-squares x y)
我有以下代码: (define (howMany list) (if (null? list) 0 (+ 1 (howMany (cdr list))))) 如果我们执行以
我有点了解如何将基本函数(例如算术)转换为Scheme中的连续传递样式。 但如果函数涉及递归怎么办?例如, (define funname (lambda (arg0 arg1)
我正在尝试附加两个字符串列表 但我不知道如何在两个单词之间添加空格。 (define (string-concat lst1 lst2) (map string-append lst1
这个问题已经有答案了: How do I pass a list as a list of arguments in racket? (2 个回答) 已关闭 8 年前。 我有一个函数,它需要无限数量的
我对这段代码的工作方式感到困惑: (define m (list 1 2 3 '(5 8))) (let ((l (cdr m))) (set! l '(28 88))) ==>(1 2 3 (5 8
我正在为学校做一项计划作业,有一个问题涉及我们定义记录“类型”(作为列表实现)(代表音乐记录)。 我遇到的问题是我被要求创建一个过程来创建这些记录的列表,然后创建一个将记录添加到该列表的函数。这很简单
我有以下代码: (define (howMany list) (if (null? list) 0 (+ 1 (howMany (cdr list))))) 如果我们执行以
我正在尝试附加两个字符串列表 但我不知道如何在两个单词之间添加空格。 (define (string-concat lst1 lst2) (map string-append lst1
如何使用抽象列表函数(foldr、foldl、map 和 filter 编写函数),无需递归,消耗数字列表 (list a1 a2 a3 ...) 并产生交替和 a1 - a2 + a3 ...? 最
我试图找出在 Scheme 中发生的一些有趣的事情: (define last-pair (lambda (x) (if (null? (cdr x))
这个问题在这里已经有了答案: Count occurrence of element in a list in Scheme? (4 个答案) 关闭 8 年前。 我想实现一个函数来计算列表中元素出现
我正在尝试使用下面的代码获取方案中的导数。谁能告诉我哪里出错了?我已经尝试了一段时间了。 (define d3 (λ (e) (cond ((number? e) 0) ((e
我是一名优秀的程序员,十分优秀!