Question

哪位高手帮我看下我的matlab程序哪里出错了，运行不了。谢谢啊

clc
clear
syms w;
m=44000;I=2879000;
Sh=0.01;Sa=0.01;Wh=0.2741;Wa=0.5884;Ro=1.225;B=32.5;H1=1;H2=1;H3=1;H4=1;A1=1;A2=1; A3=1;A4=1;
Cl=-0.09118;Cld=-7.2464;Cm=0.01324;Cmd=0.9352;Cd=0.09855;U=10;
v10=10;
Suu=4*(1200*w/(2*pi*v10))^2/(1+(1200*w/(2*pi*v10))^2)^(4/3)*0.4*v10^2*2*pi/w;
Sww=4*(1200*w/(2*pi*v10))^2/(1+(1200*w/(2*pi*v10))^2)^(4/3)*0.4*v10^2*2*pi/w;
M=[m 0;0 I];
C=[2*Sh*m*Wh-1/2*Ro*B^2*Wh*H1, -1/2*Ro*B^3*Wa*H2; -1/2*Ro*B^3*Wh*A1, 2*Sa*I*Wa-1/2*Ro*B^4*Wa*A2];
K=[m*Wh^2-1/2*Ro*B^2*Wh^2*H4, -1/2*Ro*B^3*Wa^2*H3; -1/2*Ro*B^3*Wh^2*A4, I*Wa^2-1/2*Ro*B^4*Wa^2*A3];
Hw=(-w^2*M+1i*w*C+K)^(-1);
Sll=Ro^2*U^2*B^2*(Cl*Suu+1/4*(Cld+Cd)^2*Sww);
Smm=Ro^2*U^2*B^2*(Cm*Suu+1/4*(Cmd)^2*Sww);
Sxx=[Sll,0;0,Smm];
Syy=((-w^2*M-1i*w*C+K)^(-1))*Sxx*Hw.';
sigma2=quad(Syy,0,10)
我的问题可以简化为：
clc
clear
w=sym('w');
Suu=2*w;
M=[2 0;0 3];
C=[5 4;2 3];
K=[2 7;8 2];
Hw=(-w^2*M+1i*w*C+K)^(-1);
Sll=2*Suu;
Smm=3*Suu;
Sxx=[Sll 0;0 Smm];
Syy=((-w^2*M-1i*w*C+K)^(-1))*Sxx*Hw.';
sigma2=int(Syy,w,0,10)

错误提示为：
??? Error using ==> mupadinit at 35
Unable to initialize kernel

Error in ==> sym.sym>sym.sym/symr at 252
[S,err] = mupadmex(' ',9,x);

Error in ==> sym.sym>sym.sym/symnumeric at 152
S{k} = symr(double(x(k)));

Error in ==> sym.sym>sym.sym at 116
S = cell2sym(S,symnumeric(x,a));

Error in ==> sym.mtimes at 8
A = sym(A);
请各位高手指点其中的错误及改正方法。

Answer 1

你的公式应该写错了，Suu和Sww似乎完全一样
=====================================
你那个Hw应该是自己构建的一个函数吧？
二次型之后求积分不是所有函数都能求出解析解的
这不一定是程序出问题，有可能是算法本身就不对

Answer 2

你的公式没错么？ 挨行运行一下的话你的Hw数值
Hw =

[ (- 8528159124494628956369266752693348662244073731737375473664000000000*w^6 - 12380977601206541953095044273483367978653714212275236110336000000*i*w^5 + 4990549318479859735471438691035268596961214207190123930208698368000*w^4 + 1581583182235756273456744400172724198686993569996940483850403840*i*w^3 - 915035855514724894555245930860394249228407493962922558225324703744*w^2 + 329006630691047429361819236694767586407845016991079156519272448*i*w + 43808099672010008806976756722999221496354659557459819917370982400)/(375239001477763674080247737118507341138739244196444520841216000000000000*w^8 - 247149661908037673843896476632703881577987112502599940365581549568000000*w^6 + 56388101964802810245692940986570903206844046628286842224039446511106624*w^4 - 4890643570056508395408011388631047431988746555285293030560197295475404*w^2 + 142026825638946551113325961208191554980794269544895836191903537130625), (36647341086075437200740895120542654827602651253148278063104000000*i*w^5 + 26196975811216848524897190482618026410262066107872772645453824000*w^4 - 15088210209520337371007274356140299525463720832811368796513107968*i*w^3 - 7624350345421494062135747611985225054449215275504358252875874304*w^2 + 919322937828028650102038748047901485552319025176899413486862336*i*w + 419513884193898829640533702502313973057231821723760681864396800)/(375239001477763674080247737118507341138739244196444520841216000000000000*w^8 - 247149661908037673843896476632703881577987112502599940365581549568000000*w^6 + 56388101964802810245692940986570903206844046628286842224039446511106624*w^4 - 4890643570056508395408011388631047431988746555285293030560197295475404*w^2 + 142026825638946551113325961208191554980794269544895836191903537130625)]
[ -(706398690289271679249343325137984029997400064*w^2 + 53072305784372120418215736574720973360594944)/(15526694757753819378582479437299240646666616832000000*i*w^5 - 6219058297241482253868244536613906263253914222592000*w^4 - 5775521007771350276231086331324279751647195593787744*i*w^3 + 1584156151845743907977946092831159317065462451001872*w^2 + 342798149565372139155713618175797295308685074879818*i*w - 82797887274642026682748885692892558052061807502775), ((706398690289271679249343325137984029997400064*w^2 + 53072305784372120418215736574720973360594944)*(- 19619932779746576182128811180032000*i*w^3 + 5349339825418572605549915921902532*w^2 + 1444580213050153743728242643499089*i*w - 398099443569567870705236651955232))/((2569863464629712977509710345249856*w^2 + 193075923687036703698015977623876)*(15526694757753819378582479437299240646666616832000000*i*w^5 - 6219058297241482253868244536613906263253914222592000*w^4 - 5775521007771350276231086331324279751647195593787744*i*w^3 + 1584156151845743907977946092831159317065462451001872*w^2 + 342798149565372139155713618175797295308685074879818*i*w - 82797887274642026682748885692892558052061807502775))]
额的神啊……

Answer 3

最好把乘除全换成点成的形式。“.*”"./"