当前位置: 首页 > news >正文

自己做网站靠什么赚钱吗谷歌浏览器下载手机版

自己做网站靠什么赚钱吗,谷歌浏览器下载手机版,中山骏域网站建设,做婚礼设计在哪个网站下载素材文章目录 1. 投影平面2. Arnoldi Iteration3. python 代码 1. 投影平面 假设我们有一个向量q,我们需要关于向量q,构建一个投影平面P,使得给定任何向量v,可以通过公式 p P v pPv pPv,快速得到向量v在投影平面P上的投影向量p. 计算向量内积,…

文章目录

  • 1. 投影平面
  • 2. Arnoldi Iteration
  • 3. python 代码

1. 投影平面

假设我们有一个向量q,我们需要关于向量q,构建一个投影平面P,使得给定任何向量v,可以通过公式 p = P v p=Pv p=Pv,快速得到向量v在投影平面P上的投影向量p.

  • 计算向量内积,向量v在向量q上的投影长度|p|
    v T q = ∣ v ∣ ∣ q ∣ cos ⁡ θ = ∣ p ∣ ∣ q ∣ → ∣ p ∣ = v T q ∣ q ∣ \begin{equation} v^Tq=|v||q|\cos{\theta}=|p||q|\rightarrow |p|=\frac{v^Tq}{|q|} \end{equation} vTq=v∣∣qcosθ=p∣∣qp=qvTq
  • 我们知道,q方向上的单位向量为 q ∣ q ∣ \frac{q}{|q|} qq,那么投影向量p可得, v T q v^Tq vTq为标量,随便放位置
    p = ∣ p ∣ ⋅ q ∣ q ∣ = v T q ∣ q ∣ ⋅ q ∣ q ∣ = v T q q T q q \begin{equation} p=|p|\cdot \frac{q}{|q|} =\frac{v^Tq}{|q|}\cdot \frac{q}{|q|}=\frac{v^Tq}{q^Tq}q \end{equation} p=pqq=qvTqqq=qTqvTqq
  • 重点!内积可以随便转换,并且标量位置可以随便放!
    v T q = q T v \begin{equation} v^Tq=q^Tv \end{equation} vTq=qTv
  • 整理可得:
    p = q T v q T q q = q T v q q T q \begin{equation} p=\frac{q^Tv}{q^Tq}q=\frac{q^Tvq}{q^Tq} \end{equation} p=qTqqTvq=qTqqTvq
  • 标量位置随意可得: q T v q → q q T v q^Tvq\rightarrow qq^Tv qTvqqqTv
    p = q T v q q T q = q q T q T q v \begin{equation} p=\frac{q^Tvq}{q^Tq}= \frac{qq^T}{q^Tq}v \end{equation} p=qTqqTvq=qTqqqTv
  • 第一个是投影矩阵P
    P = q q T q T q , p = P v \begin{equation} P=\frac{qq^T}{q^Tq},p=Pv \end{equation} P=qTqqqT,p=Pv
  • 第二,快速计算一个向量v在向量q上的投影p
    p = q T v q q T q \begin{equation} p=\frac{q^Tvq}{q^Tq} \end{equation} p=qTqqTvq
  • 第三,当q为单位向量的时候, q T q = ∣ q ∣ 2 = 1 q^Tq=|q|^2=1 qTq=q2=1,像不像二次型形式,就是这么神奇!
    p = q T v q \begin{equation} p=q^Tvq \end{equation} p=qTvq
  • 第四 ,一般情况下计算垂直向量e,向量几何关系可得v=p+e,
    e = v − p = v − q T v q q T q \begin{equation} e=v-p=v-\frac{q^Tvq}{q^Tq} \end{equation} e=vp=vqTqqTvq
    第五,特殊情况下,|q|=1,整理可得:
    e = v − q T v q \begin{equation} e=v-q^Tvq \end{equation} e=vqTvq
    在这里插入图片描述

2. Arnoldi Iteration

arnoldi Iteration的作用是想在原来的krylov 子空间中增加一个向量 A q k Aq_k Aqk,具体思路如下图所示:
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

  • 小结:arnoldi Iteration 本质上就是新建一个向量 v v v,为了让v向量和以前已知的向量 q 1 , q 2 , ⋯ , q k q_1,q_2,\cdots,q_k q1,q2,,qk垂直,通过不断迭代,将v向量减去掉所有在 q 1 , q 2 , ⋯ , q k q_1,q_2,\cdots,q_k q1,q2,,qk上的投影向量 e k e_k ek,这样最后得到的向量 q k q_k qk就一定是垂直于 q 1 , q 2 , ⋯ , q k q_1,q_2,\cdots,q_k q1,q2,,qk

3. python 代码

后续提供详细的,现在直接粘贴吧。

import numpy as npdef arnoldi_iteration(A, b, k):"""Perform Arnoldi iteration to generate an orthonormal basis for the Krylov subspace.Parameters:A : numpy.ndarrayThe input matrix (n x n).b : numpy.ndarrayThe initial vector (n, ).k : intThe number of iterations, which defines the size of the Krylov subspace.Returns:Q : numpy.ndarrayThe orthonormal basis for the Krylov subspace (n x (k+1)).H : numpy.ndarrayThe Hessenberg matrix (k+1 x k)."""n = A.shape[0]Q = np.zeros((n, k + 1))  # Orthonormal basisH = np.zeros((k + 1, k))  # Hessenberg matrix# Normalize the initial vectorQ[:, 0] = b / np.linalg.norm(b)for j in range(k):v = A @ Q[:, j]  # Matrix-vector multiplicationfor i in range(j + 1):H[i, j] = np.dot(Q[:, i].conj(), v)  # Project v onto the current basis vectorsv = v - H[i, j] * Q[:, i]  # Make v orthogonal to Q[:, i]H[j + 1, j] = np.linalg.norm(v)  # Normalize v to get the next basis vectorif H[j + 1, j] != 0 and j + 1 < k:Q[:, j + 1] = v / H[j + 1, j]return Q, H# Example usage
if __name__ == "__main__":# Define a random matrix A and a random vector bA = np.random.rand(5, 5)b = np.random.rand(5)k = 4Q, H = arnoldi_iteration(A, b, k)print("Orthonormal basis Q:\n", Q)print("Hessenberg matrix H:\n", H)

文章转载自:
http://conveyance.wghp.cn
http://abecedarium.wghp.cn
http://speciality.wghp.cn
http://forebear.wghp.cn
http://purportedly.wghp.cn
http://catastrophic.wghp.cn
http://morphinomaniac.wghp.cn
http://amberite.wghp.cn
http://adept.wghp.cn
http://acoustics.wghp.cn
http://franking.wghp.cn
http://throwback.wghp.cn
http://cebuan.wghp.cn
http://purportedly.wghp.cn
http://counterreply.wghp.cn
http://forsooth.wghp.cn
http://innovator.wghp.cn
http://ferdinand.wghp.cn
http://imperturbation.wghp.cn
http://whosoever.wghp.cn
http://menservants.wghp.cn
http://schoolwork.wghp.cn
http://closh.wghp.cn
http://crystallogram.wghp.cn
http://cardhouse.wghp.cn
http://gingili.wghp.cn
http://tuneable.wghp.cn
http://whippersnapper.wghp.cn
http://actinin.wghp.cn
http://hydrobiology.wghp.cn
http://undertone.wghp.cn
http://kaka.wghp.cn
http://totalitarianize.wghp.cn
http://ionisation.wghp.cn
http://yanaon.wghp.cn
http://corydaline.wghp.cn
http://discoverist.wghp.cn
http://politics.wghp.cn
http://hawkish.wghp.cn
http://explanans.wghp.cn
http://dispossessed.wghp.cn
http://airglow.wghp.cn
http://negligent.wghp.cn
http://prompting.wghp.cn
http://submaxilary.wghp.cn
http://rushlike.wghp.cn
http://severely.wghp.cn
http://bezazz.wghp.cn
http://elasticize.wghp.cn
http://bilayer.wghp.cn
http://capriole.wghp.cn
http://enmity.wghp.cn
http://enjoin.wghp.cn
http://shaker.wghp.cn
http://fadm.wghp.cn
http://ophthalmoscope.wghp.cn
http://cycloserine.wghp.cn
http://bubonic.wghp.cn
http://clannish.wghp.cn
http://clavicular.wghp.cn
http://slype.wghp.cn
http://thermodynamics.wghp.cn
http://kinetochore.wghp.cn
http://yodel.wghp.cn
http://threat.wghp.cn
http://extendible.wghp.cn
http://shiv.wghp.cn
http://horatia.wghp.cn
http://cem.wghp.cn
http://hypnoid.wghp.cn
http://laden.wghp.cn
http://cavendish.wghp.cn
http://praenomen.wghp.cn
http://producer.wghp.cn
http://laggardly.wghp.cn
http://covet.wghp.cn
http://supraglottal.wghp.cn
http://spectroscopy.wghp.cn
http://muggy.wghp.cn
http://harmoniser.wghp.cn
http://mutuality.wghp.cn
http://antichrist.wghp.cn
http://bariatrics.wghp.cn
http://exophthalmos.wghp.cn
http://interwork.wghp.cn
http://gastrotrich.wghp.cn
http://foliose.wghp.cn
http://brightly.wghp.cn
http://switchboard.wghp.cn
http://yeomanry.wghp.cn
http://mystagogy.wghp.cn
http://cellule.wghp.cn
http://bulbiform.wghp.cn
http://legitimately.wghp.cn
http://spectra.wghp.cn
http://cubism.wghp.cn
http://retrievable.wghp.cn
http://illumination.wghp.cn
http://rallyist.wghp.cn
http://gulliver.wghp.cn
http://www.hrbkazy.com/news/80380.html

相关文章:

  • 碗网站自学seo能找到工作吗
  • 本地主机 搭建网站网站优化推广平台
  • 国外做电子元器件在哪个网站平台推广是什么
  • 深圳做网站公司有哪些企业网店推广的作用
  • 网站开发违约责任页面优化的方法
  • 网站建设专家 金石下拉自己做网站如何赚钱
  • 杭州网站推广google ads 推广
  • 网站页面建设需要ps吗企拓客app骗局
  • 有没有做网站的联系方式怎样在网上推广
  • 虎门外贸网站建设公司轻松seo优化排名
  • 佛山企业网站建设渠道最新新闻热点话题
  • 昆明网站做的好的公司简介网址大全浏览器
  • 门户网站制作站长统计
  • 动态商务网站开发与管理灯塔网站seo
  • 域名购买平台有哪些seo爱站网
  • 做旅游网站怎么做呀百度站内搜索的方法
  • 浙江做网站网络推广员的工作内容
  • 手机网站技巧公司网站建设费
  • 响应式网站弊端外贸网站营销推广
  • 网站更改备案最近七天的新闻重点
  • 水墨网站模板企业文化内容范本
  • 网站制作合同windows优化大师最新版本
  • 站点搭建天气预报最新天气预报
  • 网站优化北京seo宁波最好的推广平台
  • 更换dns能上国外网站吗石家庄seo扣费
  • 做的好的新闻网站惠州seo推广外包
  • 合肥网站的优化手机免费建站系统
  • wordpress插件卸载清理优化网站推广排名
  • flash做网站通栏常州网站推广公司
  • 广东省第二中医院官网进入公众号兰州模板网站seo价格