| 查看: 573 | 回复: 0 | ||
[求助]
在matlab中将含有变量“w”的表达式存入矩阵元素，无法生成矩阵。哪里出问题了？
|
|
我是想生成矩阵tg，然后det（tg）=0，求解w的值。但是一直无法生成矩阵。求大神指教 clc; tic; epssys=1.0e-6; syms w a1=w*0.003; a2=w*1.376e-004; kt2=w*3.315e-004; b2 = kt2; r=0.01; bh2 = @(nu,z)besselj(nu,z)-1i*bessely(nu,z); tg=zeros(5,5); for n2=1:5 for n3=1:5 n=n2-3; n1=n3-3; A11=n*besselj(n,a1*r)/r-a1*besselj(n+1,a1*r); B11=n*bh2(n,a1*r)/r-a1*bh2(n+1,a1*r); C21=1i*n*besselj(n,b2*r); A121=0; A122=(2*n*lml1*a1+2*lml1*a1)*besselj(n+1,a1*r)/r; A123=lml1*a1^2*besselj(n+2,a1*r); A12=A121-A122+A123; B121=0; B122=(2*n*lml1*a1+2*lml1*a1)*bh2(n+1,a1*r)/r; B123=lml1*a1^2*bh2(n+2,a1*r); B12=B121-B122+B123; C221=(2*lmn2*1i*n^2-2*lmn2*1i*n)*besselj(n,b2*r)/(kt2*r^2); C222=2*lmn2*1i*n*b2*besselj(n+1,b2*r)/(kt2*r); C22=C221-C222; C131=2*lmn2*(1i*n^2+n^2)*besselj(n,a2*r)/r^2; C132=2*lmn2*1i*n*a2*besselj(n+1,a2*r)/r; C13=C131-C132; C341=b2^2*lmn2*n*besselj(n,b2*r)/(kt2*r); C342=b2^3*lmn2*besselj(n+1,b2*r)/kt2; C34=C341-C342; F5=C22*C13*C34; F10=C21*C13*C34; Tn1=A11*(F1-F5)-A12*(F7-F10); Tn2=B12*(F7-F10)-B11*(F1-F5); Tn=vpa(Tn1/Tn2); tg_1=Tn; tg(n2,n3)=str2double(vpa(tg_1)); end; end; tg |
回复此楼» 猜你喜欢
材料调剂
已经有6人回复
-
化工求调剂
已经有7人回复
-
321求调剂
已经有10人回复
-
材料专硕322分
已经有13人回复
-
材料295
已经有3人回复
-
11408 359分 求调剂
已经有3人回复
-
320分，材料与化工专业，求调剂
已经有24人回复
-
286求调剂
已经有7人回复
-
化学调剂
已经有10人回复
-
357求调剂
已经有3人回复
找到一些相关的精华帖子，希望有用哦~
matlab里如何实现用向量赋值数组元素并且时数组实现运算？
已经有10人回复
- matlab循环生成的字符，如何随着逐步放入字符矩阵中？
已经有6人回复
- matlab 将 数值或者矩阵copy到粘贴板上用什么命令
已经有4人回复
- 一个随机矩阵线性相关的概率
已经有5人回复
- MATLAB子矩阵的问题
已经有11人回复
- 我用matlab求矩阵的特性向量矩阵，怎么出错了？
已经有13人回复
- Matlab 7.0 基础教程 清华大学
已经有29人回复
- matlab中fmincon的fun是矩阵运算怎么办？
已经有5人回复
- matlab求解正交矩阵。
已经有10人回复
- matlab简单矩阵运算问题请教
已经有7人回复
- 循环计算出的矩阵的保存
已经有10人回复
- 请教matlab用数组给参数矩阵赋值的问题
已经有11人回复
- 请教MATLAB中的LMI求解问题，急！！！
已经有7人回复
- matlab求逆矩阵是不是有问题？
已经有13人回复
- matlab循环计算转为矩阵 30金币求助
已经有7人回复
- 如何用Matlab画矩阵的散点图？
已经有3人回复
- 菜鸟求助~matlab怎么样让一个1xn的向量中的每个元素是个矩阵？
已经有11人回复
- matlab中系数是矩阵的积分
已经有5人回复
- matlab产生随机二维矩阵，这些随机数只能是-1或者1？
已经有9人回复
- 再次求助MATLAB的值构成矩阵问题。。。
已经有6人回复
- 【求助】关于matlab的高维数组
已经有6人回复
- 【求助】matlab曲面拟合表达式
已经有10人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
