哪位高手帮我看下我的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);
请各位高手指点其中的错误及改正方法。

你的公式应该写错了,Suu和Sww似乎完全一样
=====================================
你那个Hw应该是自己构建的一个函数吧?
二次型之后求积分不是所有函数都能求出解析解的
这不一定是程序出问题,有可能是算法本身就不对
温馨提示:内容为网友见解,仅供参考
第1个回答  2010-12-29
你的公式没错么? 挨行运行一下的话你的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))]
额的神啊……本回答被网友采纳
第2个回答  2010-12-29
最好把乘除全换成点成的形式。“.*”"./"

大侠们帮忙看看,MATLAB中编的程序哪里出错了
程序哪里出错了?程序主要存在的错误有,1、e应该为轮换搜索的次数 2、用solve()函数求解的t, s 为全部解,其解既有实数,又有虚数,所以在执行 C=X+t*e1 、C=X+s*e2 命令时出现了错误(软件不知道要用你哪个值去计算)3、因而声明变量时,应变量t, s 为实数解,syms t s real 4、...

哪位大神帮忙看一下matlab 程序哪里出错了 为什么总:显示未定义函数或变...
fprintf(fid,'%12.8f\\r',strain_rate);fclose(fid);这个中的strain_rate,应该是strain_rate(m,n)这种形式吧,因为在循环外,m,n应该是具体的数值。我不太懂MATLAB程序,但是按计算机程序来看的话,应该是这个问题。

matlab 打开出错
一、系统的运行库比较旧,而软件需要的是更新版本的运行库;或者系统里根本就没有软件需要的运行库; 二、系统的运行库损坏的说; 三、软件需要的运行库与系统的语言版本不一致。如图零、查毒杀毒,清理系统插件; 一、如果是IE出现这个错误,你应该检查插件了(例如用360),以及尝试打开IE菜单“工具”、“internet选项”...

matlab 优化设计部分,程序如下,大神看看哪里出错了,我是新手。
exitflag的值是-2,表示没有找到可行解。所谓优化,是要在满足所有约束条件的解的集合里面找到目标函数取极值的解,而现在的问题是,求解器根本就找不到能够满足全部约束的解,而lb、ub也是一种约束,不满足是很正常的。比较大的可能是,题主所给的约束条件本身是存在冲突的,最好认真检查一下约束条件...

MATLAB程序出错了,怎么样修改?
题主提供的代码没有错误,但运行时出错,其问题在于,题主把自定义函数与主程序的运行顺序搞反了。方法1、分别将y1=Q(r)、y2=R(q)和y3=S(r)自定义函数保存为Q.m、R.m和S.m文件,再将主程序保存为main.m文件,然后在命令窗口中,执行 >>main %回车 方法2、建立main.m文件,按下列顺序...

请高手帮我看一下我的matlab程序错在哪里了? max(min{f1,f2,f3,f4,f...
(1)要满足整数解的要求,7个变量应该有7个方程的约束,你刚加了两个方程不满足整数解是很正常的;(2)怎么知道最小的最大值就一定为-7呢?(另,我这里试过,加入第4项约束时,最小最大值的确还是-7)帮你改写了一下非线性约束条件:function [C,ceq] = ghun(x)C = []; ceq = ...

matlab lsqcurvefit函数问题, 运行不出结果,一直出错,高手帮帮忙吧,谢...
disp('a=',num2str(a));%显示结果 你的这个语句错了,程序无法运行的。另外不知道你的a和t是什么东西,无法进行下一步的验证了。你把这个改好之后运行下试试看吧。不行把a和t贴出来看看

Matlab7运行不了……
应该是CPU的问题,你可以用下面的方法试试:1.先确定你安装程序的位置,假设在C盘。C:\\MATLAB7\\bin\\win32文件夹中有个atlas_Athlon.dll的文件。2.右击我的电脑选择"属性",在"高级"中点"环境变量",在"系统变量"中点"新建",在其中输入:变量名:BLAS_VERSION 变量址:C:\\MATLAB7\\bin\\win32\\...

我的matlab程序问什么出现这种情况:??? Error using ==> imread_百度...
会给出提示:IMREAD Read image from graphics file。A = IMREAD(FILENAME,FMT)就是读取图片的功能,所以你得有图片。在上面给出的D:\\Matlab7.0\\ORL\\ORL\\ORL\\s110\\位置整一个1.pgm就ok了,这个错误属于调用的错误,跟程序无关。加一个1.pgm应该就没问题了。你试一下看看 不可以的话再联系。

求大神看看matlab函数哪里出错了
首先,数组1\/a应该写成1.\/a;其次,因为有1.\/a和1.\/(1-a),所以,a=0:0.001:1不对,可以改成a=0.1:0.001:0.9;具体代码如下:clear clc a=0.1:0.001:0.9;m=4*pi*10e-7;I=1;t1=m*I\/(2*pi);t2=(1-a).\/a.*sqrt(a.^2+(1-a).^2);t3=a.\/(1-a).*sqrt(a....

相似回答