matlab里a1不能做变量,在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
matlab里a1不能做变量,在matlab中将含有变量“w”的表达式存入矩阵元素,无法生成矩阵。哪里出问题了?...相关推荐
- matlab里面的cov函数,如何强制Matlab /八度cov函数使用可选参数
我正在使用Matlab和Octave的cov(协方差)功能.实际上,我使用的是Octave,但最终它必须同时适用于这两种方法.此函数具有可选的第二个或第三个参数,以指示是否应使用N或N-1进行归一化. ...
- matlab里excel汉字怎么显示,[求助]matlab如何导入excel中的汉字,求助~~~
第一种方法,你可以使用xlsread函数来读取excel中的数据 第二种方法,就是把字符转化为数字,使用函数str2num xlsread的使用 EXAMPLES: 1. Default operat ...
- matlab里如何输入圆周,如何在MATLAB中打印小圆圈?
问题我正在尝试将小圆(直径= 0.02厘米等)实心圆圈打印成精确尺寸的pdf.我可以打印直径为2厘米,1厘米,0.5厘米的精确圆圈,但是当我试图低于~0.04厘米时,打印的图形不包括所有圆圈.让我告诉 ...
- 如何在matlab里输入复杂公式_在MATLAB作图中输出漂亮的公式
函数 语法 效果 语法 效果 语法 效果 \sin\theta \cos\theta \tan\theta \arcsin\frac{L}{r} \arccos\frac{T}{r} \arctan\ ...
- matlab神经网络 时间序列,请问吧里有大神做过MATLAB时间序列神经网络(NARX)吗?...
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 请问吧里有大神做过MATLAB时间序列神经网络(NARX)吗?请教一下该神经网络的预测问题 我用网上的一个案例:知道2015年降雨我要预测该年水位.用往年 ...
- matlab narxnet,请问吧里有大神做过MATLAB时间序列神经网络(NARX)吗?
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 请问吧里有大神做过MATLAB时间序列神经网络(NARX)吗?请教一下该神经网络的预测问题 我用网上的一个案例:知道2015年降雨我要预测该年水位.用往年 ...
- Matlab里evalin和assignin的用法
原文地址:Matlab里evalin和assignin的用法作者:了凡春秋 assignin与evalin MATLAB通常的基本工作空间是'base'空间.MATLAB在程序运行过程中,将为 ...
- MATLAB里面的filter和filtfilt的C语言源代码
MATLAB里面的filter和filtfilt的C语言源代码 嗯,算法非常简单,就是网上搜不到C代码实现.filter是个很万能的数字滤波器函数,只要有滤波器的差分方程系数,IIR呀FIR呀都能通过 ...
- matlab中ln4怎么表示,matlab里ln怎么表示
matlab中图中颜色和形状表示_数学_自然科学_专业资料.颜色字符串有'c', 'm', 'y', 'r', 'g', 'b', 'w',和'k'.分别表示青,红紫,黄,红,绿,白...... 上下 ...
最新文章
- 仅展示近三天的动态设置_微信朋友圈设置三天可见了吗?
- MATLAB编写ode文件,MATLABODE45问题M文件为br/functiondq 爱问知识人
- LeetCode: 105. Construct Binary Tree from Preorder and Inorder Traversal
- 将tomcat的session信息通过memcached实现共享
- DataGridView使用技巧十:单元格表示值的自定义
- Thymeleaf3语法详解
- R语言第九讲 验证集法
- 大数据驱动乡村振兴共享共治机制研究
- 10.2 SQ02维护信息集
- 打开文件管理器_会声会影影片配置文件管理器
- Excel 2007中的新文件格式
- 容器技术Docker K8s 44 Serverless Kubernetes(ASK)详解-ASK弹性伸缩管理
- access统计行数 vb_VB代码行数统计工具
- 相机标定+极线矫正+elas生成视差图
- maxscale的causal_reads参数
- yum install 命令报错Could not retrieve mirrorlist
- 假若明天来临——《AI.未来》读后感3900字
- Struts2 框架项目新建教程(strut 2.5.20)(基于IDEA)
- cad延伸快捷键_史上最全CAD快捷键大全,熟练掌握画图快得飞起来(强烈建议收藏)...
- python+selenium自动化脚本例子
热门文章
- 几种机器学习平台的对比和选择
- linux下查看某个文件或目录占用磁盘空间的大小
- 【转】Linux 移动或重命名文件/目录-mv 的10个实用例子
- 新概念英语(1-37)Making a bookcase
- 结对编程1-基于GUI的四则运算生成器
- struts2基础----自定义拦截器
- 1.Weblogic通Eclipse调试配置(Weblogic同Eclipse调试配置技术)
- APUE 学习笔记(十) 高级I/O
- ZF环境要求及如何配置
- 大数据之-Hadoop3.x_MapReduce_序列化案例FlowDriver---大数据之hadoop3.x工作笔记0100