我是想生成矩阵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

matlab里a1不能做变量,在matlab中将含有变量“w”的表达式存入矩阵元素,无法生成矩阵。哪里出问题了?...相关推荐

  1. matlab里面的cov函数,如何强制Matlab /八度cov函数使用可选参数

    我正在使用Matlab和Octave的cov(协方差)功能.实际上,我使用的是Octave,但最终它必须同时适用于这两种方法.此函数具有可选的第二个或第三个参数,以指示是否应使用N或N-1进行归一化. ...

  2. matlab里excel汉字怎么显示,[求助]matlab如何导入excel中的汉字,求助~~~

    第一种方法,你可以使用xlsread函数来读取excel中的数据 第二种方法,就是把字符转化为数字,使用函数str2num xlsread的使用 EXAMPLES: 1. Default operat ...

  3. matlab里如何输入圆周,如何在MATLAB中打印小圆圈?

    问题我正在尝试将小圆(直径= 0.02厘米等)实心圆圈打印成精确尺寸的pdf.我可以打印直径为2厘米,1厘米,0.5厘米的精确圆圈,但是当我试图低于~0.04厘米时,打印的图形不包括所有圆圈.让我告诉 ...

  4. 如何在matlab里输入复杂公式_在MATLAB作图中输出漂亮的公式

    函数 语法 效果 语法 效果 语法 效果 \sin\theta \cos\theta \tan\theta \arcsin\frac{L}{r} \arccos\frac{T}{r} \arctan\ ...

  5. matlab神经网络 时间序列,请问吧里有大神做过MATLAB时间序列神经网络(NARX)吗?...

    该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 请问吧里有大神做过MATLAB时间序列神经网络(NARX)吗?请教一下该神经网络的预测问题 我用网上的一个案例:知道2015年降雨我要预测该年水位.用往年 ...

  6. matlab narxnet,请问吧里有大神做过MATLAB时间序列神经网络(NARX)吗?

    该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 请问吧里有大神做过MATLAB时间序列神经网络(NARX)吗?请教一下该神经网络的预测问题 我用网上的一个案例:知道2015年降雨我要预测该年水位.用往年 ...

  7. Matlab里evalin和assignin的用法

    原文地址:Matlab里evalin和assignin的用法作者:了凡春秋 assignin与evalin     MATLAB通常的基本工作空间是'base'空间.MATLAB在程序运行过程中,将为 ...

  8. MATLAB里面的filter和filtfilt的C语言源代码

    MATLAB里面的filter和filtfilt的C语言源代码 嗯,算法非常简单,就是网上搜不到C代码实现.filter是个很万能的数字滤波器函数,只要有滤波器的差分方程系数,IIR呀FIR呀都能通过 ...

  9. matlab中ln4怎么表示,matlab里ln怎么表示

    matlab中图中颜色和形状表示_数学_自然科学_专业资料.颜色字符串有'c', 'm', 'y', 'r', 'g', 'b', 'w',和'k'.分别表示青,红紫,黄,红,绿,白...... 上下 ...

最新文章

  1. 仅展示近三天的动态设置_微信朋友圈设置三天可见了吗?
  2. MATLAB编写ode文件,MATLABODE45问题M文件为br/functiondq 爱问知识人
  3. LeetCode: 105. Construct Binary Tree from Preorder and Inorder Traversal
  4. 将tomcat的session信息通过memcached实现共享
  5. DataGridView使用技巧十:单元格表示值的自定义
  6. Thymeleaf3语法详解
  7. R语言第九讲 验证集法
  8. 大数据驱动乡村振兴共享共治机制研究
  9. 10.2 SQ02维护信息集
  10. 打开文件管理器_会声会影影片配置文件管理器
  11. Excel 2007中的新文件格式
  12. 容器技术Docker K8s 44 Serverless Kubernetes(ASK)详解-ASK弹性伸缩管理
  13. access统计行数 vb_VB代码行数统计工具
  14. 相机标定+极线矫正+elas生成视差图
  15. maxscale的causal_reads参数
  16. yum install 命令报错Could not retrieve mirrorlist
  17. 假若明天来临——《AI.未来》读后感3900字
  18. Struts2 框架项目新建教程(strut 2.5.20)(基于IDEA)
  19. cad延伸快捷键_史上最全CAD快捷键大全,熟练掌握画图快得飞起来(强烈建议收藏)...
  20. python+selenium自动化脚本例子

热门文章

  1. 几种机器学习平台的对比和选择
  2. linux下查看某个文件或目录占用磁盘空间的大小
  3. 【转】Linux 移动或重命名文件/目录-mv 的10个实用例子
  4. 新概念英语(1-37)Making a bookcase
  5. 结对编程1-基于GUI的四则运算生成器
  6. struts2基础----自定义拦截器
  7. 1.Weblogic通Eclipse调试配置(Weblogic同Eclipse调试配置技术)
  8. APUE 学习笔记(十) 高级I/O
  9. ZF环境要求及如何配置
  10. 大数据之-Hadoop3.x_MapReduce_序列化案例FlowDriver---大数据之hadoop3.x工作笔记0100