24小时热门版块排行榜    

查看: 4562  |  回复: 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

新虫 (小有名气)

引用回帖:
6楼: Originally posted by guokeqin at 2016-09-17 06:57:58
solve函数,不能计算点乘,建议你采用solve函数计算出来函数关系式,然后手动输入方程再操作

我把点乘去掉显示:
在赋值 A( = B 中,A 和 B 中的元素数目必须相同。
我重新编了一下solve的公式:
kx=solve(tan(t*(((sqrt(((w(m)/c)^2)*(a1-((a4^2)/a1)))))-(kx^2)))==(a1*sqrt(((w(m)/c)^2)*(a1-((a4^2)/a1))-(kx^2))*sqrt((kx^2)-((w(m)/c)^2)))/((a1*((w(m)/c)^2))-(kx^2)-(a4*kx*sqrt((kx^2)-((w(m)/c)^2)))),kx);
但是在kx赋值到Kx的时候出现错误,另外,solve语句单独拿出来进行计算的时候会出现警告:
警告: Cannot solve symbolically. Returning a numeric approximation instead. 这个时候我能输出kx。很疑惑。
7楼2016-09-17 22:28:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

duguwuhao

新虫 (小有名气)

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

duguwuhao

新虫 (小有名气)

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

发自小木虫Android客户端
3楼2016-09-16 14:43:59
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +24 医学老男孩 2026-08-13 50/2500 2026-08-15 06:38 by 学员8dgXkO
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 04:45 by 4wMiSEwB6436
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-15 00:16 by 4wMiSEwB6436
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 HFw0lei2R37i 2026-08-14 8/400 2026-08-15 00:16 by 4wMiSEwB6436
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +5 Ldrop2023 2026-08-13 5/250 2026-08-14 18:48 by ssxclkj
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +6 majunge000 2026-08-11 8/400 2026-08-14 14:44 by 医学老男孩
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
信息提示
请填处理意见