Singular matrix C in LSQ subproblem是什么原因 80C

这里使用的是scipy中的optimize优化函数

具体代码如下:
from scipy import optimize as opt
import numpy as np
def func(X,G):
sum=0.0
for i in np.arange(G.shape[0]):
sum=sum+X[3*i]
return -1*sum
def sigma(X,G,i):
sum=0.0
for j in np.arange(G.shape[0]):
sum=sum+G[i,j]*X[3*j+2]
return sum*X[3*i+2]
def Cons1(G):
cons=()
for i in np.arange(G.shape[0]):
cons=np.append(cons,{'type':'eq','fun':lambda X:X[3*i]-X[3*i+1]-sigma(X,G,i)})
return cons
def Cons2(G,I):
cons=()
for i in np.arange(G.shape[0]):
for j in np.arange(i,G.shape[0]):
cons=np.append(cons,{'type':'ineq','fun':lambda X:I[i,j]-(X[3*i+2]-X[3*j+2])*G[i,j]})
cons=np.append(cons,{'type':'ineq','fun':lambda X:(X[3*i+2]-X[3*j+2])*G[i,j]+I[i,j]})
return cons
def bou(G,P,vn,Pg):
bounds=np.mat(np.zeros([3*G.shape[0],2]))
for i in np.arange(G.shape[0]):
bounds[3*i,:]=Pg[i,:]
bounds[3*i+1,0]=0
bounds[3*i+1,1]=P[0,i]
bounds[3*i+2,0]=0.95*vn
bounds[3*i+2,1]=1.05*vn
return bounds
def x_start(G,P,vn,Pg):
x0=np.array(np.zeros([1,3*G.shape[0]]))
for i in np.arange(G.shape[0]):
x0[0,i]=float((Pg[i,0]+Pg[i,1]))/2
x0[0,i+1]=P[0,i]
x0[0,i+2]=vn
return x0
def opf(G,I,P,vn,Pg):
b=bou(G,P,vn,Pg)
cons=Cons1(G)
cons=np.append(cons,Cons2(G,I))
x=x_start(G,P,vn,Pg)
res=opt.minimize(func, args=(G),x0=x,bounds=b,constraints=cons,method='SLSQP',options={'disp': True})
return res.x
#以下为输入数据:
G=np.mat([[1/0.64,0,-1/0.64,0,0],
[0,1/0.8+1/0.8+1/0.5,-1/0.8,-1/0.8,-1/0.5],
[-1/0.64,-1/0.8,1/0.64+1/0.8+1/1.6,-1/1.6,0],
[0,-1/0.8,-1/1.6,1/1.6+1/0.8+1,-1],
[0,-1/0.5,0,-1,1/0.5+1]])
I=np.mat([[0,0,1000,0,0],
[0,0,1000,1000,1000],
[1000,1000,0,1000,0],
[0,1000,1000,0,1000],
[0,1000,0,1000,0]])
P=np.mat([0,0,1000,800,0.1])
vn=250
Pg=np.mat([[300,400],[1100,1200],[0,0],[0,0],[200,300]])
print(opf(G,I,P,vn,Pg))

1个回答

qq_38756935
baba_822 这个链接里面并没有说怎么解决,也是提出了一个和我类似的问题!
一年多之前 回复
Csdn user default icon
上传中...
上传图片
插入图片
抄袭、复制答案,以达到刷声望分或其他目的的行为,在CSDN问答是严格禁止的,一经发现立刻封号。是时候展现真正的技术了!
其他相关推荐
SVT:A Singular Value Thresholding Algorithm For Matrix Completion.pdf
这是一篇经典的矩阵填充算法文章,用到的是SVT,即奇异值阈值法。
Lsq一个MVC-CRUD例子
Lsq一个MVC-CRUD例子 希望可以帮助大家
matrix[0,:] 是什么意思
matrix[0,:]=line[0:3] 其中matrix[0 , :]代表什么?
singular control systems
广义系统的第一本著作,是研究广义系统必读的英文书籍。
矩阵链乘MATRIX C/C++
矩阵链乘实验,实现矩阵链乘算法,n=23,24,25,…,28,矩阵大小为p0p1…pn,随机产生正整数,范围5-120
C sharp matrix library
C# matrix library used for VS.NET C# matrix library used for VS.NET
C - Matrix(匈牙利 最大匹配数)
Give you a matrix(only contains 0 or 1),every time you can select a row or a column and delete all the '1' in this row or this column . Your task is to give out the minimum times of deleting all the '...
Meschach Matrix Computations in C
C语言的矩阵库——Meschach的说明书,Meschach Matrix Computations in C
乱序处理器中的LSQ简介
乱序处理器中访存指令的操作我一直看的不太懂,其中lsq在流水线中是怎么使用的我也不太清楚,今天看了篇论文,讲的是AMD的经典处理器K7,记录一下: LSU部件可以理解为存储器子系统的最高层,在该部件中包含Load Queue与Store Queue。其中Load Queue与Store Queue之间有着强烈的耦合关系,因此许多处理器系统将其合称为LSQ。在多数处理器的存储器子系统中,LSU是最顶...
是什么原因
编的程序在visual studio里全部运行通过,调试也没有问题。可做成安装程序安装后,一调试就出错误!生成安装程序都是按网上教程做的,也是生成成功的啊
<>是什么原因
刚学ASP,碰到个头大的问题:rnrn网站的第一个登录界面没什么问题,登录好后在做任何查询,都产生rnHTTP 500 - 内部服务器错误rnrn请问可能是什么原因,应该怎么调试
Singular Control Systems
1989年,Dm L Y系统地总结了众学者和他本人的研究成果,出版了{Singular Control Systems》一书 详细地介绍了广义系统的基本理论
是什么原因?
我设置了一个JOB按时间来执行一个存储过程,现在JOB也按时的运行但是总是报错!所以我的存储过程也没法工作。上面的情况到底是什么错误原因!存储过程看上去也没有呀!help me!!
Meschach: Matrix Computations in C
C语言矩阵运算库Meschach使用文档,由David E. Stewart和Zbigniew Leyk编写。 清晰,带目录,
singular value decomposition tutorial
singular value decomposition tutorial,适用于engineering的学生,在工程中有广泛应用
Cygwin完全安装以后的Singular CAS: 到底什么是Singular?
如下图所示,为什么rxvt,emacs都有存在于Singular菜单里面? 这个Singular是X11的概念呢,还是Cygwin实现的一种东东呢?rn[img=https://img-bbs.csdn.net/upload/201308/26/1377517168_458710.png][/img]
题目中选项C错的原因是什么
一道2011年的面试题,选出A、B、C、D中错误的选项,其中rn A. define N 10; int x[N];rn B. int N=10; int x[N]; rn C. int x[0..10];rn D. int x[];rnrnB和D知道是明显错的,但C说不清楚错误的原因,想知道0..10表示的到底是什么意思
C:/WINDOWS/winsg3a.dll加载出错是什么原因?
C:/WINDOWS/winsg3a.dll加载出错是什么原因?打开电脑后一段时间一出现,该怎么解决?
不能格式化C盘,是什么原因呀?
我想重装系统,可是不能格式化C盘呀。rn我用的是WIN98SE版本呀。rn
AS不能新建cpp和c文件是什么原因
这是我的AS右键打开的样子rn[img=https://img-bbs.csdn.net/upload/201705/18/1495117639_587582.png][/img]rn一般情况下应该是这样:rn[img=https://img-bbs.csdn.net/upload/201705/18/1495117742_494345.png][/img]rn而且cpp文件在项目的目录里被识别为Singleton文件rnrn我应该怎么设置我的AS来解决这个问题
Convolution Calder´on-Zygmund singular integral
Convolution Calder´on-Zygmund singular integralConvolution Calder´on-Zygmund singular integral
The Singular Value Decomposition, Applications and Beyond
The Singular Value Decomposition, Applications and Beyond 张志华的大牛作,适合有英文基础的
C语言处理Matrix的代码
一个C语言 Matrix 处理类代码
Singular Value Decomposition source code
powerful set of techniques for dealing with sets of equations or matrices. In many cases where Gaussian elimination and LU decomposition fail to give satisfactory results, this set of techniques, known as singular value decomposition, or SVD, will diagnose for you precisely what the problem is. In some cases, SVD will not only diagnose the problem, it will also solve it.
AN IMAGE WATERMARKING METHOD BASED ON THE SINGULAR
A robust image watermarking scheme in which a binary image is embedded in the singular values of selected DWT blocks in the horizontal and vertical sub-bands of a 1-level decomposition of a gray-scale image is proposed. The embedded blocks are selected by a secret key to enhance imperceptibility. A watermarked image that is perceptually indistinguishable from the original is obtained. The watermarking retrieval is non-blind and requires the use of parameters extracted during the watermarking process. The performance of the proposed algorithm is tested by comparing the retrieved watermark to the original watermark. Computer simulation results show that the algorithm is robust to common signal processing attacks such as Gaussian noise, cropping, and low pass filtering. It is also resistant to JPEG compression.
A multilinear singular value decomposition.pdf
We discuss a multilinear generalization of the singular value decomposition. There is a strong analogy between several properties of the matrix and the higher-order tensor decomposition; uniqueness, link with the matrix eigenvalue decomposition, first-order perturbation effects, etc., are analyzed. We investigate how tensor symmetries affect the decomposition and propose a multilinear generalization of the symmetric eigenvalue decomposition for pair-wise symmetric tensors
Matrix
Java写的仿Matrix的小程序<br>(骇客帝国中电脑屏幕上出现的神秘文字)
matrix
 用于矩阵计算  我不太会用
The Matrix
The Matrix
这个错误是什么原因
rn bmp:=Tbitmap.Create;rn mem:=cds.CreateBlobStream(cds.Fieldbyname('Graphic'),bmread);rn mem.Position:=0;rn[color=#FF0000] bmp.LoadFromStream(mem);[/color] 当运行到这里的时候 出现raised exception class EinvalidGraphic with message 'bitmap image is not valid'.process stoppedrn image1.Picture.Assign(bmp);rn bmp.Free;rn mem.Free ;rnrnrn如何解决
cmake编译错误是什么原因???
1>CMake is re-running because D:/MyGUI/Build/CMakeFiles/generate.stamp dependency file is missing.rn1>CMake Error: The source directory "D:/mygui" does not appear to contain CMakeLists.txt.rn1>Specify --help for usage, or press the help button on the CMake GUI.rn1>Project : error PRJ0019: 工具从"Checking Build System"
请问是什么原因?谢谢
访问http://localhost:8080/my/HelloWorldrn HelloWorld放在webapps/my/WEB-INF/classes下rnrnweb.xml里的内容rn rnrnrn rn contentTypern text/html;charset=Gb2312rn rn rn rn rn kkkrn HelloWorldrn rn rn kkkrn /HelloWorldrn rnrnrn显示错误:rnHTTP Status 404 - /my/HelloWorldrnrn--------------------------------------------------------------------------------rnrntype Status reportrnrnmessage /my/HelloWorldrnrndescription The requested resource (/my/HelloWorld) is not available.rnrnrn--------------------------------------------------------------------------------rnrnApache Tomcat/6.0.18rnrn但是可以访问html文件
提示它的原因是什么
[code=php][/code]rn上面的写法为啥当文件test.txt是0字节的时候提示这个呢rn[b][color=#FF0000]Warning: fread() [function.fread]: Length parameter must be greater than 0 in D:\APMServ5.2.6\www\htdocs\01\file.php on line 7[/color][/b]rnrn另外上面和下面这两种读取方式哪个更好点rn[code=php]",$con);rn echo $con;rn?>[/code]
接口错误,是什么原因?
internal interface IFormViewrn rn abstract event UpdateCommandEventHandler UpdateCommand;rnrn ViewController Controller get; rnrn System.Windows.Forms.Form Form get; rnrn IFormPanel TaskPanel get; rnrn rnrn public delegate void UpdateCommandEventHandler(object sender, EventArgsViewCommand e);rn rn rnrn错误:rnUpdateCommandEventHandler 接口不能声明类型,应怎么改?
错误提示是什么原因?
C:\Documents and Settings\俊\桌面\1\RPG.cpp(588) : error C2440: '=' : cannot convert from 'void *' to 'struct HBRUSH__ *'rn Conversion from 'void*' to pointer to non-'void' requires an explicit castrnError executing cl.exe.rnrn请问什么错误?怎么解决?rnrn
晕啊,是什么原因??????
晕啊,是什么原因??????rnINSERT INTO reg(sex) VALUES ('"&Option1(0).Caption &"')rn其中SEX为CHAR类型,不知道为什么,而Option1(0).Caption可以读出来值="男"rn但就是叉不进去哦!
这个异常是什么原因?
在测试一个BMP EJB时,遇到这样的问题:在得到一个引用以后,我先调用了它的一个set方法,接着调用了自己写的一个Update方法,其大致内容如下:rn tryrn rn String sql="update production set a=?,b=?where ID=? ";rn ps = connect.prepareStatement(sql);rn ps.setLong(1,this.a);rn ps.setString(2,this.b);rn ps.setLong(3,this.ID);rn int intResult=ps.executeUpdate();rn if(intResult!=1)rn throw new SQLException("更新记录出错!");rn rn rn catch(Exception ex)rn rn rn throw new EJBException("Exception in Update:"+ex);rn rn就发生了以下异常:rnRemoteException occurred in server thread; nested exception is: com.ibm.ejs.container.UncheckedException: ; nested exception is: javax.ejb.EJBException: Exception in ejbStore:com.ibm.websphere.ce.cm.StaleConnectionException: Connection is closed rn看得出来是和数据库的连接关闭了,但是我在调用这个方法之前只用了一个set方法,而且程序里面也没有显示的关闭连接的代码,为什么就关闭了连接啊?rn大虾们帮忙看看!
这种死锁的原因是什么
一个网站程序,访问比较频繁;主要是update更新数据的时候,发生“事务(进程 ID 75)与另一个进程已被死锁在 lock | communication buffer 资源上,且该事务已被选作死锁牺牲品。请重新运行该事务。”rnrn用网上找到的一个查找死锁原因的sql试了,结果如下:rnrn进程号SPID:55被进程号SPID:65阻塞(死锁),其当前进程执行的SQL语法如下rnEventType Parameters EventInfornPRC Event 0 sp_cursor;1
这个错误的原因是什么?
1 我的SQL7 在2000P下没有问题,可是在2000S下打开表(open table ->return all rows)的时候总是出现一个错误"无法定位程序输入点.....(一长串字符) 于mdt2fw.DLL 上" ,使用"SQL Server Query Analyzer"查询数据却没有问题, 这是什么原因?是不是sql7 不能在2000server下正常运行? rn2 以前针对sql7做的程序,如果改用SQL2000来做后台数据库,程序需不需要改动?rnrn盼高手指点。rn
相关热词 c# stream 复制 android c# c#监测窗口句柄 c# md5 引用 c# 判断tabtip 自己写个浏览器程序c# c# 字符串变成整数数组 c#语言编程写出一个方法 c# 转盘抽奖 c#选中treeview