matlab中funcl,求大神帮忙解答如何在Matlab function中使用fsolve函数

matlab function中的代码如下：

function [CT_,alpha_] = fcn(P,L_,V0)

coder.extrinsic('fsolve');

coder.extrinsic('fun');

global m n v;

m=P;

n=L_;

v=V0;

a=zeros(1,2);

a=fsolve('fun',[4,1.5],optimset('Display','off'));  ****错误位置在这里

CT_=a(2);

alpha_=a(1);

fun函数代码如下：(基本就是一个插值运算)

function [ e ] = fun( x )

global m n v;

a=[-10 -2 6 12 18 24 28];

Ct=[0 0.97 1.93 2.90 3.85];

Cx=[0.25 0.325 0.5 0.65 0.875 1.080 1.14;-0.12 0.25 0.6 0.93 1.24 1.34 1.41;-0.47 0.125 0.64 1.125 1.52 1.64 1.72;-0.8 -0.41 0.525 1.12 1.68 1.97 2.08;-1.125 -0.425 0.325 0.975 1.75 2.35 2.6];

Cxe=[0.009 0.007 0.01 0.0180 0.0795 0.0075 0.0080;-0.0060 0.0075 0.012 0.016 0.024 0.015 0.013;-0.0075 0.0125 0.015 0.022 0.0265 0.015 0.011;-0.0080 0.015 0.016 0.022 0.0265 0.0215 0.0155;-0.0040 0.0160 0.0155 0.0200 0.0250 0.0250 0.0215];

Cxsp=[0.001 0.001 0.0005 0 -0.00033 -0.00033 0;0.00133 0.00116 0.00017 -0.00033 -0.00067 -0.00017 0.00083;0.0015 0.00133 -0.00008 -0.00050 -0.00083 0 0.0015;0.00167 0.0020 0.0020 0.0010 0.00017 0 0;0.00167 0.00167 0.00067 -0.00017 -0.00067 0 0];

Cxf=[0.005 0.005 0.005 0.005 0.005 0.005 0.005;0.01750 0.01250 0.01150 0.01250 0.01250 0.01250 0.01250;0.02750 0.02 0.01750 0.02 0.02 0.02 0.01750;0.02750 0.02500 0.02500 0.02500 0.03 0.02750 0.02750;0.02750 0.03 0.03250 0.03250 0.03750 0.03500 0.03500];

Cz=[-0.1 1 2.1 2.3 2.2 2 1.61;1.9 3.15 4.1 4.5 4.49 3.6 3.15;3 4.15 5.32 5.8 6.1 5.45 4.5;3.71 4.9 6.15 6.85 7.35 7.1 6.25;4.2 5.55 6.9 7.7 8.4 8.8 8.15];

Cze=[0.06 0.0344 0.0396 0.0286 0.04 -0.005 0.005;0.0290 0.0460 0.0520 0.0440 0.0430 -0.0120 -0.0100;0.0350 0.0570 0.0570 0.0500 0.0440 0.0100 -0.0100;0.0490 0.0650 0.0570 0.0540 0.0355 0.0260 0.0050;0.0710 0.0700 0.0640 0.0480 0.0446 0.0410 0.0120];

Czsp=[0.00333 -0.00667 -0.00667 -0.00667 -0.00333 -0.00167 0;0.00133 -0.00917 -0.00950 -0.00883 -0.00530 -0.00300 -0.00133;0.00033 -0.01160 -0.01 -0.01160 -0.00750 -0.00750 -0.003;0.0005 -0.01 -0.01 -0.01160 -0.00916 -0.00583 0.00083;0.00133 0.00116 -0.00667 -0.00867 -0.00925 -0.01030 0];

Czf=[0.005 0.0035 0.005 0.005 0.01 0.01 0.015;0.01750 0.01250 0.01250 0.01150 0.01250 0.01 0.015;0.03 0.02250 0.02000 0.01750 0.01500 0.01 0.015;0.03500 0.02250 0.02500 0.02 0.01750 0.01250 0.01250;0.04250 0.03250 0.02750 0.02250 0.02 0.01500 0.01];

rou=1.205;

S=55.7;

ELE=-1.4;

Cx0=interp2(a,Ct,Cx,x(1),x(2),'spline');

Cxe0=interp2(a,Ct,Cxe,x(1),x(2),'spline');

Cxsp0=interp2(a,Ct,Cxsp,x(1),x(2),'spline');

Cxf0=interp2(a,Ct,Cxf,x(1),x(2),'spline');

CD=(Cx0+Cxf0*5+Cxe0*ELE+Cxsp0*10);

Cz0=interp2(a,Ct,Cz,x(1),x(2),'spline');

Cze0=interp2(a,Ct,Cze,x(1),x(2),'spline');

Czsp0=interp2(a,Ct,Czsp,x(1),x(2),'spline');

Czf0=interp2(a,Ct,Czf,x(1),x(2),'spline');

CL=(Cz0+Czf0*5+Cze0*ELE+Czsp0*10);

X=-CD*0.5*rou*v^2*S;

L=CL*0.5*rou*v^2*S;

e(1)=m-X;

e(2)=n-L;

错误提示：Call to MATLAB function aborted: 在赋值 A(I)=B 中，B 和 I 中的元素数目必须相同。

全局变量data store memory也都设置完成了。不知道问题出在哪里，按理说fsolve输出的就是一个二维向量，请大神帮忙解答万分感谢

捕获.JPG (34.08 KB, 下载次数: 2)

2017-10-27 20:22 上传