이 질문을 팔로우합니다.
- 팔로우하는 게시물 피드에서 업데이트를 확인할 수 있습니다.
- 정보 수신 기본 설정에 따라 이메일을 받을 수 있습니다.
inner matrix dimensions must agree error in the function?
조회 수: 1 (최근 30일)
이전 댓글 표시
function [f] = objfun(W,VCVMatx)
%function for varience of the portfolio with swtN stocks
f =W*VCVMatx*W';
this is my function I am trying to run in the main program. However, when run in the main program, it keeps giving me :
Error using * Inner matrix dimensions must agree.
W is a 10(N)X1 vector, while VCVMatx is an NXN vector?
Does anyone see anything wrong, considering I multiplied it by a transposed W'
I tried running fmincon, but in the main program it says the function matrix dimensions do not agree.
f = (W'*VCVMatx)*W;
i=1;
while i<=rows(ERGrid);
ExpectedReturn = ERInc(i);
[W] = fmincon('objfun',x0,swtOther_A,swtOther_b,swtOther_Aeq,swtOther_beq,swtMinWt1,swtMaxWt1,[],options);
ERGrid(2:swtN+1,i) = W;
i=i+1;
end;
W' is a 1X10 vector while VCVMatx is a 10X10 matrix and W is 10X1 vector
fmincon cannot run as a result.
답변 (1개)
Geoff Hayes
2014년 5월 30일
The first matrix product (in the equation) is the problem. Since W is Nx1 and VCVMatx is NxN, then W*VCVMtx will fail due to the incompatible dimensions of the two matrices. If you are just expecting a single value, then the equation should be re-written as
f = W'*VCVMtx*W;
댓글 수: 29
Geoff Hayes
2014년 5월 30일
Type the commands that Star Strider indicated just prior to evaluating f to ensure that the dimensions of the matrices are compatible.
Geoff Hayes
2014년 5월 30일
편집: Geoff Hayes
2014년 5월 30일
Try the following to get the number of rows and columns of VCVMatx
VCVMatx_size = size(VCVMatx)
Sameer
2014년 5월 30일
i did that, and it says
Inner matrix dimensions must agree.
Error in objfun (line 8) f = (W'*VCVMatx)*W;
so it is basically giving me the same error as before.
Star Strider
2014년 5월 30일
I’ve never seen the size function return ‘undefined’ as output.
Go back and see how you created VCVMtx. If the code creating it is not very long, post it or the relevant part of it.
Sameer
2014년 5월 30일
swtN = cols(Data); VCVMatx = cov(Data)*252;
data is imported from a .csv file and has 10 (swtN) stocks.
these equation are in the main program
Geoff Hayes
2014년 5월 30일
@Sameer - you said a couple of comments ago that W is 1x10 (which is different from the original statement of 10x1). If that is the case, and VCVMtx is still 10x10 (what did size(VCVMtx) return?), then your original line of code (not the one I suggested above) should be valid
f = W*VCV*W';
Sameer
2014년 5월 30일
when i type that same equation in the main program ( f = W'*VCVMtx*W;) it seems to work.
but it does not work when it is in the function.
Geoff Hayes
2014년 5월 30일
Strange that it doesn't work in the function. Could you attach the function (m file) to a comment, using the paper clip button?
Geoff Hayes
2014년 5월 30일
편집: Geoff Hayes
2014년 5월 30일
The code is
function [f] = objfun(W,VCVMatx)
%function for varience of the portfolio with swtN stocks
global VCVMatx;
W_size = size(W);
VCVMatx_size = size(VCVMatx);
f = (W'*VCVMatx)*W;
end
Why is there a global variable with the same name as the input parameter? If this global variable has not been set then it is probably empty and so the error is thrown. Please comment out the line
% global VXVMatx;
and try again so that the second input (to the function) is used.
Sameer
2014년 5월 30일
i took out the globa variable in the function and now it is saying:
Error using objfun (line 7) Not enough input arguments.
Error in fmincon (line 564) initVals.f = feval(funfcn{3},X,varargin{:});
Error in FEHW7 (line 183) [W] = fmincon('objfun',x0,swtOther_A,swtOther_b,swtOther_Aeq,swtOther_beq,swtMinWt1,swtMaxWt1,[],options);
Caused by: Failure in initial user-supplied objective function evaluation. FMINCON cannot continue.
Geoff Hayes
2014년 5월 30일
Okay - I don't have the fmincon function, but reading the documentation at http://www.mathworks.com/help/optim/ug/fmincon.html indicates that the function handle passed into fmincon must be for a function that accepts a vector and returns a scalar. Is that why you have tried to get around the limitation of only one input being passed - by using a global variable?
If that is the case, then I suppose you have to put that global variable back and remove the second input to your objfun. Then, assuming that the input vector is 1x10 and that the global variable VCVMtx is 10x10, then your f should be changed to
f = (W*VCVMatx)*W';
and the global variable VCVMatx must be set prior to the evaluation of fmincon.
Geoff Hayes
2014년 5월 30일
I think the problem that you were originally experiencing was the the global variable VCVMatx was not defined to be a 10x10 matrix, it was just empty. And so if you multiply a non-empty matrix ( W ) with an empty matrix ( VCVMatx ) the message Error using * Inner matrix dimensions must agree is thrown. So it didn't matter which combination of f we were using (transpose of W on left or right) it would still fail.
At some point, the global variable must have been set. And now, with the correct order of the matrix multiplications for f, it works as expected.
But you will have to remember to set the global variable VCVMatx...
Sameer
2014년 5월 30일
would you happen to know how to place two linear equality constraints in one fmincon function? thanks.
Geoff Hayes
2014년 5월 30일
No, unfortunately I have never used that function (nor have the toolbox for it).
참고 항목
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!오류 발생
페이지가 변경되었기 때문에 동작을 완료할 수 없습니다. 업데이트된 상태를 보려면 페이지를 다시 불러오십시오.
웹사이트 선택
번역된 콘텐츠를 보고 지역별 이벤트와 혜택을 살펴보려면 웹사이트를 선택하십시오. 현재 계신 지역에 따라 다음 웹사이트를 권장합니다:
또한 다음 목록에서 웹사이트를 선택하실 수도 있습니다.
사이트 성능 최적화 방법
최고의 사이트 성능을 위해 중국 사이트(중국어 또는 영어)를 선택하십시오. 현재 계신 지역에서는 다른 국가의 MathWorks 사이트 방문이 최적화되지 않았습니다.
미주
- América Latina (Español)
- Canada (English)
- United States (English)
유럽
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom (English)
아시아 태평양
- Australia (English)
- India (English)
- New Zealand (English)
- 中国
- 日本Japanese (日本語)
- 한국Korean (한국어)