| 查看: 534 | 回复: 4 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
zhanghcui木虫 (小有名气)
|
[交流]
【求助】 matlab 符号转化double
|
||
|
各位高人， 怎样把符号表达式转变为double表达式 下面是我的主程序和出错提示： syms x y z1 z2 m=3 n=3 h=0.0003; w=0.0005; a=0.0001 e=h/w %z1=@pp ff=sym(zeros(m,n)) for i=1:m %feval(z1,x,y,i) z1=x.^i+y.^i for j=1:n % z2=@pp %feval(z2,x,y,j) z2=x.^j+y.^j z1x=diff(z1,x) z2x=diff(z2,x) z1y=diff(z1,y) z2y=diff(z2,y) ff(i,j)=e^2*z1x*z2x+z1y*z2y ff(i,j)=char(ff(i,j)) ff(i,j)=strrep(ff(i,j),'*','.*') ff(i,j)=strrep(ff(i,j),'/','./') ff(i,j)=strrep(ff(i,j),'\','.\') ff(i,j)=strrep(ff(i,j),'^','.^') ss=intout(ff(i,j),0,1,0,1) end end 出错： Warning: Inputs must be character arrays or cell arrays of strings. > In zzz at 23 ??? Error using ==> strrep Conversion to double from sym is not possible. Error in ==> zzz at 23 ff(i,j)=strrep(ff(i,j),'*','.*') [ Last edited by luoqiquan on 2009-12-19 at 18:14 ] |
回复此楼» 猜你喜欢
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组招收2026硕士研究生1名
已经有0人回复
-
山西大同大学物理学硕士研究生招收调剂生
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有241人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
一志愿物理所凝聚态物理352求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
070200求调剂，一志愿某211，288分
已经有17人回复
zhanghcui
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2116.1
- 散金: 140
- 红花: 2
- 帖子: 208
- 在线: 229.4小时
- 虫号: 659790
- 注册: 2008-11-22
- 专业: 传递过程
5楼2008-12-08 09:04:02
mingdong
荣誉版主 (知名作家)
- 应助: 0 (幼儿园)
- 贵宾: 4.548
- 金币: 42102.8
- 散金: 72
- 红花: 8
- 帖子: 5922
- 在线: 456.4小时
- 虫号: 633286
- 注册: 2008-10-22
- 性别: GG
- 专业: 能源化工
- 管辖: 计算模拟
★
xuefei06(金币+1,VIP+0):感谢参与并提意见
xuefei06(金币+1,VIP+0):感谢参与并提意见
|
你好 数值仿真版开有matlab专贴，请参考： http://muchong.com/bbs/viewthread.php?tid=759829&fpage=1 |
2楼2008-12-04 16:50:56
zhanghcui
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2116.1
- 散金: 140
- 红花: 2
- 帖子: 208
- 在线: 229.4小时
- 虫号: 659790
- 注册: 2008-11-22
- 专业: 传递过程
3楼2008-12-05 13:06:56
★ ★ ★ ★ ★
mingdong(金币+2,VIP+0):感谢您的帮助！
xuefei06(金币+3,VIP+0):感谢你给出这么详细的答案，希望对楼主有所帮助
mingdong(金币+2,VIP+0):感谢您的帮助！
xuefei06(金币+3,VIP+0):感谢你给出这么详细的答案，希望对楼主有所帮助
|
syms x y z1 z2 m=3 n=3 h=0.0003; w=0.0005; a=0.0001 e=h/w %z1=@pp ff=sym(zeros(m,n)) for i=1:m %feval(z1,x,y,i) z1=x.^i+y.^i for j=1:n % z2=@pp %feval(z2,x,y,j) z2=x.^j+y.^j z1x=diff(z1,x) z2x=diff(z2,x) z1y=diff(z1,y) z2y=diff(z2,y) ff(i,j)=e^2*z1x*z2x+z1y*z2y ff(i,j)=char(ff(i,j)) ff(i,j)=strrep(ff(i,j),'*','.*') ff(i,j)=strrep(ff(i,j),'/','./') ff(i,j)=strrep(ff(i,j),'\','.\') ff(i,j)=strrep(ff(i,j),'^','.^') ss=intout(ff(i,j),0,1,0,1) end end 你去掉后面的几个转换 syms x y z1 z2 m=3 n=3 h=0.0003; w=0.0005; a=0.0001 e=h/w %z1=@pp ff=sym(zeros(m,n)) for i=1:m %feval(z1,x,y,i) z1=x.^i+y.^i for j=1:n % z2=@pp %feval(z2,x,y,j) z2=x.^j+y.^j z1x=diff(z1,x) z2x=diff(z2,x) z1y=diff(z1,y) z2y=diff(z2,y) ff(i,j)=e^2*z1x*z2x+z1y*z2y ff(i,j)=char(ff(i,j)) end end 输出结果是 ff = [ 34/25, 18/25*x+2*y, 27/25*x^2+3*y^2] [ 18/25*x+2*y, 36/25*x^2+4*y^2, 54/25*x^3+6*y^3] [ 27/25*x^2+3*y^2, 54/25*x^3+6*y^3, 81/25*x^4+9*y^4] 不知道是不是你相要的， 你要是想吧*换成别的运算符号 你再做一个替代循环吧 另外intout有这个函数吗？你是不是写错了 |
4楼2008-12-05 15:16:12
10