| 查看: 4564 | 回复: 10 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
求助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; 跪求大神解惑。。。  |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有50人回复
-
filecode
已经有11人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助MATLAB编程
已经有8人回复
- 求助matlab写变限积分
已经有4人回复
- 求助Matlab
已经有0人回复
- 求助matlab图像
已经有0人回复
- 求助matlab求解三元超越方程
已经有4人回复
- 求助matlab三维图纠正程序
已经有2人回复
- 求助matlab有限差分程序
已经有1人回复
- 求助matlab写循环
已经有6人回复
- 求助MATLAB画图，二维曲线，总是出错，求助高手
已经有3人回复
- matlab求助
已经有4人回复
- Matlab求助
已经有6人回复
- 求助matlab问题
已经有5人回复
- 求助matlab程序
已经有3人回复
- 求助matlab问题
已经有5人回复
- 高金求助matlab解微分方程组
已经有12人回复
- 求助matlab程序
已经有4人回复
- 求助matlab积分的问题
已经有4人回复
- matlab求助
已经有1人回复
|
大神，这个我了解了，所以尝试把solve语句修改了一下，单独计算kx=solve(tan(t*(((sqrt(((w/c)^2)*(a1-((a4^2)/a1)))))-(kx^2)))==(a1*sqrt(((w/c)^2)*(a1-((a4^2)/a1))-(kx^2))*sqrt((kx^2)-((w/c)^2)))/((a1*((w/c)^2))-(kx^2)-(a4*kx*sqrt((kx^2)-((w/c)^2)))),kx);时，现出现警告： 警告: Cannot solve symbolically. Returning a numeric approximation instead. 这个时候能输出kx的数值解，所以我就比较疑惑为什么有数值解还不能赋值成功。。。 |
8楼2016-09-17 22:31:36
2楼2016-09-16 13:07:16
3楼2016-09-16 14:43:59
4楼2016-09-16 15:38:26
