| 查看: 1121 | 回复: 4 | ||
susulucky金虫 (小有名气)
|
[求助]
matlab方程求教
|
|
函数部分： function epsiul_m =epsiul( Lambda ) syms epsiul_m w %%%%%%% parameters of gold epsiul_L=5.9673; w_D=2113.6*2*pi*1e12; r_D=15.92*2*pi*1e12; omigal_L=650.07*2*pi*1e12; T_L=104.86*2*pi*1e12; deltepsiul=1.09; c=3*1E8; [epsiul_m,w] = solve('epsiul_m=epsiul_L-w_D^2/(w*(w+j*r_D))-deltepsiul*omigal_L^2/((w^2-omigal_L^2)+j*T_L*w)','w=2*pi*c/(Lambda*(sqrt(real(epsiul_m)/epsiul_0)))'); end epsiul_0=(1/(36*pi))*1e-9; 程序部分： Lambda=1.55e-6; epsiul_m =epsiul( Lambda); n_gold=sqrt(epsiul_m/epsiul_0); 为什么运行后最后出来的结果还是包含epsiul_m的表达式呢，怎么修改直接出来一个数值，而不包含epsiul_m和其它已知数的表达式，，比如最后结果连pi都未带入，本来结果应该就是一个数值呢，而不是一个表达式呢。我不知道我哪里出错了。 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : epsiul.m
- 附件 2 : n.m
2013-06-17 09:46:11, 562 bytes
2013-06-17 09:46:16, 129 bytes
» 猜你喜欢
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有5人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有7人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有7人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab解偏微分方程求助
已经有12人回复
- 求Matlab解非线性方程代码
已经有12人回复
- 用matlab如何解方程？
已经有6人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- matlab求解非线性方程组，错误提示怎么解决
已经有5人回复
- 各位朋友，如何用matlab编写这个方程，并求解？
已经有15人回复
- 求助 MATLAB解方程组-fslove
已经有7人回复
- matlab求解方程中的参数
已经有21人回复
- matlab拟合方程参数时初值的选择
已经有15人回复
- matlab解微分方程
已经有10人回复
- 求matlab 解方程的问题
已经有5人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlab求解非线性方程组
已经有16人回复
- matlab解微分方程组
已经有15人回复
- 求教origin拟合简单的渗流方程（又称逾渗方程），做复合材料的应该会遇到吧
已经有20人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- MATLAB解方程出现??? Subscripted assignment dimension mismatch.
已经有3人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】求教 多元指数方程如何解 有软件否？
已经有6人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 求助一个差分方程的MATLAB实现
已经有3人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
- 【求助】新手求教 关于MATLAB 解超越方程 和画图
已经有15人回复
susulucky
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1089.7
- 散金: 1028
- 红花: 3
- 帖子: 261
- 在线: 91.6小时
- 虫号: 1043976
- 注册: 2010-06-19
- 性别: MM
- 专业: 光学
2楼2013-06-20 15:25:22
兵封年代
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1786.8
- 红花: 5
- 帖子: 92
- 在线: 52.5小时
- 虫号: 1884381
- 注册: 2012-07-09
- 性别: GG
- 专业: 能源化工
3楼2013-06-20 15:51:08
susulucky
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1089.7
- 散金: 1028
- 红花: 3
- 帖子: 261
- 在线: 91.6小时
- 虫号: 1043976
- 注册: 2010-06-19
- 性别: MM
- 专业: 光学
|
syms epsiulm w epsiul_L=5.9673; w_D=2113.6*2*pi*1e12; r_D=15.92*2*pi*1e12; omigal_L=650.07*2*pi*1e12; T_L=104.86*2*pi*1e12; deltepsiul=1.09; c=3*1E8; epsiul_0=(1/(36*pi))*1e-9; Lambda=1.55e-6; [epsiulm,w]= solve('epsiulm-epsiul_L-w_D^2/(w*(w+i*r_D))-deltepsiul*omigal_L^2/((w^2-omigal_L^2)+i*T_L*w)=0','(w*Lambda/(2*pi*c))^2-real(epsiulm)/epsiul_0=0'); n_gold=sqrt(epsiulm/epsiul_0); 不用函数其实应该是这样子的。但是出来的结果总是还包含epsiulm。本来结果应该直接就是一个数值的 |
4楼2013-06-21 11:15:56
兵封年代
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1786.8
- 红花: 5
- 帖子: 92
- 在线: 52.5小时
- 虫号: 1884381
- 注册: 2012-07-09
- 性别: GG
- 专业: 能源化工
5楼2013-06-21 15:55:58
10