diag

Matlab(9)——矩阵变换

由 巧了我就是萌 提交于 2020-03-17 03:42:54
Matlab(9)——矩阵变换 文章目录 Matlab(9)——矩阵变换 一、对角阵 1.提取矩阵对角线上的元素 2.构造对角矩阵 3.应用 二、三角阵 1.上三角矩阵 2.下三角矩阵 三、矩阵的转置 四、矩阵的旋转 五、矩阵的翻转 六、矩阵的逆矩阵 一、对角阵 1.提取矩阵对角线上的元素 diag(A)：提取矩阵A对角线上的元素，形成一个列向量 diag(A,k)： 提取A 的第 k 条对角线上元素的列向量。k=0 表示主对角线，k>0 位于主对角线上方，k<0 位于主对角线下方。 2.构造对角矩阵 diag(v) ：以向量v为主对角线元素建立对角矩阵 D = diag(v,k) ：将向量 v 的元素放置在第 k 条对角线上。k=0 表示主对角线，k>0位于主对角线上方，k<0 位于主对角线下方。 3.应用 现有一n阶方阵A （1）要将A第一行元素乘r1，第二行元素乘r2，…，第n行元素乘以rn 可以建立对角矩阵：B=diag(r1,r2,…rn)，再B*A （2）要将A第一列元素乘c1，第二行元素乘c2，…，第n行元素乘以cn 可以建立对角矩阵：B=diag(c1,c2,…cn)，再A*B 二、三角阵 1.上三角矩阵 triu(A)：返回矩阵 A 的上三角部分。 triu(A,k)：返回位于 A 的第 k 条对角线上以及该对角线上方的元素。(k可以为负) 2.下三角矩阵

np.diag()

由 跟風遠走 提交于 2020-02-22 13:19:55
x = np . arange ( 9 ) . reshape ( ( 3 , 3 ) ) x array ( [ [ 0 , 1 , 2 ] , [ 3 , 4 , 5 ] , [ 6 , 7 , 8 ] ] ) np . diag ( x ) array ( [ 0 , 4 , 8 ] ) np . diag ( x , k = 1 ) array ( [ 1 , 5 ] ) np . diag ( x , k = - 1 ) array ( [ 3 , 7 ] ) np . diag ( np . diag ( x ) ) array ( [ [ 0 , 0 , 0 ] , [ 0 , 4 , 0 ] , [ 0 , 0 , 8 ] ] ) 来源： CSDN 作者： 侬开心就好 链接： https://blog.csdn.net/weixin_40717742/article/details/104441471

弹窗插件zDialog使用教程

由 ⅰ亾dé卋堺 提交于 2020-02-08 22:29:58
1.首先现在好zDialog然后复制项目中 2.配置zDialog解压以后images文件夹位置 images存放位置根据自己实际项目而定，zDialog.js中配置位置即可，如： var IMAGESPATH = '/images/'; //图片路径配置 3.弹出框html以及js代码编写： <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd"> <html xmlns="http://www.w3.org/1999/xhtml"> <head> <meta http-equiv="Content-Type" content="text/html; charset=UTF-8"> <meta name="Copyright" content="酷站代码 http://www.5icool.org/" /> <script type="text/javascript" src="/js/jquery.min.js"></script> <script type="text/javascript" src="/js/datepattern.js"></script> <script type=

MATLAB矩阵的翻转求逆旋转转置

由 浪子不回头ぞ 提交于 2020-01-10 03:46:52
对角阵 对角矩阵 只有对角线上有非零元素的矩阵 数量矩阵 对角线上的元素相等的对角矩阵 单位矩阵 对角线上都为1的对角矩阵 （1）提取矩阵的对角线上的元素 diag(A) 提取矩阵A主对角线元素，产生一个列向量。 diag(A,k) 提取第k条对角线上的元素 k=0是中间的那条对角线 （2）构造对角矩阵 diag(V) 产生一个主对角线为V的对角阵 diag(V,k) 向量v为第k条对角线 例子 先建立5阶矩阵A，然后将A的第一行元素乘以1，第二行乘以2，……，第五行乘以5. >> A=[4 5 6 7 8;9 5 3 1 6;4 6 2 3 6;5 3 5 4 6;5 6 4 6 4] A = 4 5 6 7 8 9 5 3 1 6 4 6 2 3 6 5 3 5 4 6 5 6 4 6 4 >> D=diag(1:5) D = 1 0 0 0 0 0 2 0 0 0 0 0 3 0 0 0 0 0 4 0 0 0 0 0 5 >> F=D*A F = 4 5 6 7 8 18 10 6 2 12 12 18 6 9 18 20 12 20 16 24 25 30 20 30 20 三角阵 上三角阵 矩阵的对角线一下的元素全为零的矩阵 （1）ttiu(A) triu(A,k) 下三角阵 …………………上……………………… (2)tril(A) tril(A,k) 矩阵的转置

数学建模系列：层次分析

由 为君一笑 提交于 2019-12-21 20:19:16
前言 主要是通过matlab求解微分方程的实例，先来看一个小小问题： \[ x^2+y+(x-2y)y^{'}=0 \] 可以通过 matlab 的 dsolve 函数求得通解,两行代码就能够解决问题，关于 dsolve 的具体用法，参考官方文档的介绍 clc,clear syms x %定义符号变量 y=dsolve('x^2+y+(x-2*y)*Dy==0','x') %运行结果 % y = % x/2 + ((4*x^3)/3 + x^2 + C1)^(1/2)/2 % x/2 - ((4*x^3)/3 + x^2 + C1)^(1/2)/2 % C1是常数 问题一 求解常微分方程组的通解 代码： clc,clear syms f(x) g(x) [f1,g1]=dsolve('D2f+3*g==sin(x)','Dg+Df==cos(x)'); f1=simplify(f1) %化简 g1=simplify(g1) % 运行结果 % f1 = % C2 - sin(x)/3 + t*cos(x) + (3^(1/2)*C3*exp(3^(1/2)*t))/3 - (3^(1/2)*C4*exp(-3^(1/2)*t))/3 % g1 = % % sin(x)/3 - (3^(1/2)*C3*exp(3^(1/2)*t))/3 + (3^(1/2)*C4*exp(-3^(1

数学建模系列：层次分析

由 梦想与她 提交于 2019-12-06 04:16:54
前言 主要是通过matlab求解微分方程的实例，先来看一个小小问题： \[ x^2+y+(x-2y)y^{'}=0 \] 可以通过 matlab 的 dsolve 函数求得通解,两行代码就能够解决问题，关于 dsolve 的具体用法，参考官方文档的介绍 clc,clear syms x %定义符号变量 y=dsolve('x^2+y+(x-2*y)*Dy==0','x') %运行结果 % y = % x/2 + ((4*x^3)/3 + x^2 + C1)^(1/2)/2 % x/2 - ((4*x^3)/3 + x^2 + C1)^(1/2)/2 % C1是常数 问题一 求解常微分方程组的通解 代码： clc,clear syms f(x) g(x) [f1,g1]=dsolve('D2f+3*g==sin(x)','Dg+Df==cos(x)'); f1=simplify(f1) %化简 g1=simplify(g1) % 运行结果 % f1 = % C2 - sin(x)/3 + t*cos(x) + (3^(1/2)*C3*exp(3^(1/2)*t))/3 - (3^(1/2)*C4*exp(-3^(1/2)*t))/3 % g1 = % % sin(x)/3 - (3^(1/2)*C3*exp(3^(1/2)*t))/3 + (3^(1/2)*C4*exp(-3^(1

convert python ndarray to theano tensor type variable

由 匿名 (未验证) 提交于 2019-12-03 00:58:01
可以将文章内容翻译成中文,广告屏蔽插件可能会导致该功能失效(如失效，请关闭广告屏蔽插件后再试): 问题: I have ndarray like : diag = [] diag.append(np.diag([1,1,0])) diag.append(np.diag([0,1,1])) diag [array([[1, 0, 0], [0, 1, 0], [0, 0, 0]]), array([[0, 0, 0], [0, 1, 0], [0, 0, 1]])] How can I convert it into Theano tensor variable of type float 64, matrix ? As I need to perform dot operation like Theano.dot(diag, X) where X is shared variable of type float 64, matrix. 回答1: Just create a SharedVariable like this diag_ = theano.shared(np.array(diag).astype("float64")) theano.dot(diag_, X) http://deeplearning.net/software/theano/library/compile

oracle12c的日志查看

由 匿名 (未验证) 提交于 2019-12-03 00:40:02
查看GI日志：切换到grid用户 查看DB日志：切换到oracle的目录下 执行[oracle@swnode1 ~]$ adrci [oracle@swnode1 ~ ]$ adrci ADRCI: Release 12.2 . 0.1 . 0 - Production on Fri Jul 6 16 : 05 : 45 2018 Copyright (c) 1982 , 2017 , Oracle and/ or its affiliates. All rights reserved. ADR base = " /u01/app/oracle " adrci > show homes ADR Homes: diag /clients/user_oracle/ host_2908371209_107 diag /rdbms/orcl/gnnt1 alert日志就在上述位置 [oracle@swnode1 ~]$ tail -f /u01/app/oracle/diag/rdbms/orcl/gnnt1/alert/log.xml < txt > Please see the VKTM trace file for more details: </ txt > </ msg > < msg time =‘2018-07-05T09:12:34.163+08:00‘ org_id =

ios预览文件

由 本秂侑毒 提交于 2019-12-02 18:50:43
var diag = new Dialog(); //diag.Width = window.screen.width; diag.Height = window.screen.height; diag.URL = url; 嵌套显示会出现内容显示不全：父及div需要设置-webkit-overflow-scrolling: touch;overflow-y: scroll;height:100rem才行没高度不行包裹才行 来源： https://www.cnblogs.com/otsf/p/11758366.html

python算法对音频信号处理Sonification ：Gauss-Seidel迭代算法

由 ↘锁芯ラ 提交于 2019-12-01 07:57:49
原文链接： http://tecdat.cn/?p=7620 可以将44.1kHz单通道.wav文件中的一秒读取到长度为44100的数组（称为b）中。给定矩阵A，我们寻求系统Ax = b的解。通过Gauss-Seidel的迭代，向量如果我们将b记录的录音，则将一些白噪声作为我们的初始猜测，并在每次交替中写出Ax，我们会观察到b中高音调的音符首先变得可听，而同时白噪声的音调分解。 最初的12秒.wav文件的音频（白噪声） initialAx.wav 初始Ax，残差和残差FFT的图： ​ ​ ​ 经过一轮迭代，高音 转化gauss_seidel_out000000.wav 在光谱中可以看到一些结构： ​ ​ ​ 第二次 迭代 ： gauss_seidel_out000001.wav ​ ​ ​ 第三次 迭代 ： gauss_seidel_out000002.wav ​ ​ ​ 第四次 迭代 ： gauss_seidel_out000003.wav ​ ​ ​ 这一切都在python中完成。将.wav文件加载到数组中，在scipy中还不错。为了避免缓存问题，必须使用稀疏矩阵类，因为12秒的.wav文件需要一个大小为12 * 44100的数组。这是我使用的TridiagonalMatrix类代码片段： from numpy import * #a tridiagonal matrix
订阅 diag