24小时热门版块排行榜    

查看: 4205  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

duguwuhao

新虫 (小有名气)

[求助] 求助matlab for循环中的赋值问题。 已有2人参与

dm=0.001*2*pi*c;
    w=-0.5*2*pi*c:dm:0.5*2*pi*c;
    Kx=zeros(size(w));
    Ay=zeros(size(w));
for m=1:length(w);
    syms ay kx ky;
    y=(0.0382*8*(pi^3)*(c^4)*sqrt(ainf-(wp.^2)./(w.^2)))./((wp^2).*w);
    axx=ainf*((1-(wp*(w+1i*y)*wp))./(w.*(((w+1i*y).^2)-((wc)^2))));   %无单位
    axy=ainf*(((wp^2)*wc*1i)./(w.*(((w+y*1i).^2)-((wc)^2))));   %无单位
    azz=ainf.*((1-(wp^2))./(w.*(w+1i*y)));   %无单位
    a1=real(axx);
    a2=imag(axx);
    a3=real(axy);
    a4=imag(axy);
    a5=real(azz);
    a6=imag(azz);
    f1=ay-sqrt((kx.^2)-((w./c).^2));
    f2=kx-((sqrt(((a4.*ay./2).^2))+(a1.*(w./c).^2)-((a1.*ay.*ky)./(tan(ky.*t))))-(a4.*ay./2));
    f3=ky-(sqrt((((w./c).^2).*(a1-((a4.^2)./a1)))-kx.^2));
    [kx,ay]=solve('f1','f2','f3');
    Kx(m)=kx;
    Ay(m)=Ay;
end


运行提示:在赋值 A( = B 中,A 和 B 中的元素数目必须相同。
出错行数:    Kx(m)=kx;
    Ay(m)=Ay;

跪求大神解惑。。。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

duguwuhao

新虫 (小有名气)

修改了一下,将w改成w(m),还是出错貌似,是我solve的问题。没有得出数值解

发自小木虫Android客户端
3楼2016-09-16 14:43:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

duguwuhao

新虫 (小有名气)

那个笑脸不是我卖萌啊。。。。我也不知道怎么出现的。。。。。。
2楼2016-09-16 13:07:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

duguwuhao

新虫 (小有名气)

我检查了下,程序,貌似是slove语句的问题,单独拿出来,程序出现警报:警告: Cannot solve symbolically. Returning a numeric approximation instead. 跪求大神解决
4楼2016-09-16 15:38:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)

【答案】应助回帖

感谢参与,应助指数 +1
https://www.mathworks.com/help/symbolic/solve.html

Return Numeric Solutions
Try solving the following equation. The symbolic solver cannot find an exact symbolic solution for this equation,
and therefore issues a warning before calling the numeric solver. Because the equation is not polynomial, an
attempt to find all possible solutions can take a long time. The numeric solver does not try to find all numeric
solutions for this equation. Instead, it returns only the first solution it finds.
5楼2016-09-16 21:56:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见