引子
Matlab中有一个函数叫做find,可以很方便地寻找数组内特定元素的下标,即:Find indices and values of nonzero elements。
这个函数非常有用。比如,我们想计算图1中点Q(x0, y0)抛物线的最短距离。一个可以实施的方法是:计算出抛物线上所有点到Q点的距离,找到最小值,用find函数找到最小值对应的下标,即M点横坐标和纵坐标对应的元素的下标,M点到Q点的距离就是最短距离。
首先给出Matlab使用find函数实现的代码:
a = linspace(-5,5,1000); b = a .^2; x0 = 4; y0 = 4; dis = sqrt((a - x0).^2 + (b - y0).^2); mm = find (dis == min(dis)); a0 = a(mm); b0 = b(mm); disMin = sqrt((a0 - x0).^2 + (b0 - y0).^2); plot(a, b); hold on; scatter(x0, y0, 'k*'); scatter(a0, b0, 'k*'); xx = [a0, x0]; yy = [b0, y0]; plot(xx, yy);
NumPy中的where函数
Syntax: np.where(conditions, [x,y])
具体实现代码如下:
import numpy as np import math import matplotlib.pyplot as plt a = np.linspace(-5, 5, 10000) b = a * a x0 = 4 y0 =4 num = np.linspace(0, len(a) - 1, len(a)) dis = np.linspace(0, 0, len(a)) for k in num: k = int(k) dis[k] = dis[k] + math.sqrt((a[k] -x0) **2 + (b[k] - y0) **2) disMin = min(dis) disMinIndex = np.where(dis == disMin) disMin0 = math.sqrt((a[disMinIndex] - x0) **2 + (b[disMinIndex] - y0) **2) print('The mininum distance:',disMin) print('The mininum distance:',disMin0) print(type(dis)) a0 = a[disMinIndex] b0 = b[disMinIndex] fig = plt.figure(figsize = (6,6), dpi = 200) ax1 = plt.subplot(1,1,1) line11 = ax1.scatter(a,b,s = 1) line12 = ax1.scatter(x0, y0, s = 100, marker = '*', color = 'darkorange') line13 = ax1.scatter(a0, b0, s = 100, marker = '*', color = 'darkorange') line14 = ax1.plot([x0,a0],[y0,b0], color = 'darkorange') line15 = ax1.text(4.2,4,'Q(x0,y0)') line16 = ax1.text(0.6,5, 'M(a0,b0)') line18 = plt.xlim(-5,5) line17 = plt.ylim(0,25) plt.savefig('C:/Users/BRIAR/Desktop/index.png') plt.show()
The mininum distance: 1.943317035
The mininum distance: 1.9433170350024023
class ‘numpy.ndarray'
List中的index函数
Syntax: List.index(aimElement)
注意:此处需将NumPy数组转换成List格式的数据。
具体实现代码如下:
import numpy as np import math import matplotlib.pyplot as plt a = np.linspace(-5, 5, 10000) b = a * a x0 = 4 y0 =4 num = np.linspace(0, len(a) - 1, len(a)) dis = np.linspace(0, 0, len(a)) for k in num: k = int(k) dis[k] = dis[k] + math.sqrt((a[k] -x0) **2 + (b[k] - y0) **2) disMin = min(dis) disList = dis.tolist() disMinIndex = disList.index(disMin) disMin0 = math.sqrt((a[disMinIndex] - x0) **2 + (b[disMinIndex] - y0) **2) print('The mininum distance:',disMin) print('The mininum distance:',disMin0) print(type(disList)) a0 = a[disMinIndex] b0 = b[disMinIndex] fig = plt.figure(figsize = (6,6), dpi = 200) ax1 = plt.subplot(1,1,1) line11 = ax1.scatter(a,b,s = 1) line12 = ax1.scatter(x0, y0, s = 100, marker = '*', color = 'darkorange') line13 = ax1.scatter(a0, b0, s = 100, marker = '*', color = 'darkorange') line14 = ax1.plot([x0,a0],[y0,b0], color = 'darkorange') line15 = ax1.text(4.2,4,'Q(x0,y0)') line16 = ax1.text(0.6,5, 'M(a0,b0)') line18 = plt.xlim(-5,5) line17 = plt.ylim(0,25) plt.savefig('C:/Users/BRIAR/Desktop/index.png') plt.show()
The mininum distance: 1.943317035
The mininum distance: 1.9433170350024023
class ‘list'
以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持。
免责声明:本站文章均来自网站采集或用户投稿,网站不提供任何软件下载或自行开发的软件! 如有用户或公司发现本站内容信息存在侵权行为,请邮件告知! 858582#qq.com
稳了!魔兽国服回归的3条重磅消息!官宣时间再确认!
昨天有一位朋友在大神群里分享,自己亚服账号被封号之后居然弹出了国服的封号信息对话框。
这里面让他访问的是一个国服的战网网址,com.cn和后面的zh都非常明白地表明这就是国服战网。
而他在复制这个网址并且进行登录之后,确实是网易的网址,也就是我们熟悉的停服之后国服发布的暴雪游戏产品运营到期开放退款的说明。这是一件比较奇怪的事情,因为以前都没有出现这样的情况,现在突然提示跳转到国服战网的网址,是不是说明了简体中文客户端已经开始进行更新了呢?
更新日志
- 雨林唱片《赏》新曲+精选集SACD版[ISO][2.3G]
- 罗大佑与OK男女合唱团.1995-再会吧!素兰【音乐工厂】【WAV+CUE】
- 草蜢.1993-宝贝对不起(国)【宝丽金】【WAV+CUE】
- 杨培安.2009-抒·情(EP)【擎天娱乐】【WAV+CUE】
- 周慧敏《EndlessDream》[WAV+CUE]
- 彭芳《纯色角3》2007[WAV+CUE]
- 江志丰2008-今生为你[豪记][WAV+CUE]
- 罗大佑1994《恋曲2000》音乐工厂[WAV+CUE][1G]
- 群星《一首歌一个故事》赵英俊某些作品重唱企划[FLAC分轨][1G]
- 群星《网易云英文歌曲播放量TOP100》[MP3][1G]
- 方大同.2024-梦想家TheDreamer【赋音乐】【FLAC分轨】
- 李慧珍.2007-爱死了【华谊兄弟】【WAV+CUE】
- 王大文.2019-国际太空站【环球】【FLAC分轨】
- 群星《2022超好听的十倍音质网络歌曲(163)》U盘音乐[WAV分轨][1.1G]
- 童丽《啼笑姻缘》头版限量编号24K金碟[低速原抓WAV+CUE][1.1G]