版块导航
正在加载中...
客户端APP下载
论文辅导
申博辅导
登录
注册
帖子
帖子
用户
本版
应《网络安全法》要求,自2017年10月1日起,未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用,请尽快对帐号进行手机号验证,感谢您的理解与支持!
24小时热门版块排行榜
>
论坛更新日志
(586)
>
虫友互识
(52)
>
休闲灌水
(32)
>
导师招生
(21)
>
公派出国
(17)
>
硕博家园
(16)
>
考博
(14)
>
论文投稿
(14)
>
教师之家
(12)
>
考研
(11)
>
博后之家
(10)
>
基金申请
(10)
>
论文道贺祈福
(6)
>
招聘信息布告栏
(4)
>
留学生活
(4)
>
功能材料
(2)
小木虫论坛-学术科研互动平台
»
计算模拟区
»
计算模拟
»
matlab符号变量的一个问题
5
1/1
返回列表
查看: 995 | 回复: 4
只看楼主
@他人
存档
新回复提醒
(忽略)
收藏
在APP中查看
匿名
用户注销
(小有名气)
应助: 0
(幼儿园)
金币: 58.7
散金: 66
红花: 3
帖子: 141
在线: 95.5小时
虫号: 0
注册: 2010-09-21
性别: GG
专业: 固体与表面物理化学
本帖仅楼主可见
» 猜你喜欢
求国际会议网站
已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。
已经有0人回复
物理学I论文润色/翻译怎么收费?
已经有228人回复
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有19人回复
反铁磁体中的磁性切换:两种不同的机制已成功可视化
已经有0人回复
26申博推荐:南京航空航天大学国际前沿院光学方向招收博士生!
已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
» 本主题相关价值贴推荐,对您同样有帮助:
【活动or资源】vasp晶体结构优化
已经有142人回复
【求助】求稳定性的问题~~急~~
已经有8人回复
【整理】今天整理硬盘上的资料,偶然发现的好玩的小东西就直接发上来吧
已经有105人回复
【求助】请教一个关于解方程的问题
已经有8人回复
【求助】 求MAtlab求解程序!!!
已经有20人回复
【求助】因不能追加赏币重新发帖:急求Tchebichef矩不变量或Krawtchouk矩不变量代码
已经有19人回复
【求助】求问一个寻找最优解问题,是泛函问题么
已经有22人回复
1楼
2011-11-10 11:17:38
已阅
同方向广播
申请计算强帖
回复此楼
编辑
查看我的主页
dbb627
荣誉版主
(著名写手)
专家经验: +4
计算强帖: 12
应助: 289
(大学生)
贵宾: 0.589
金币: 24640.4
散金: 551
红花: 61
沙发: 1
帖子: 1246
在线: 1794.8小时
虫号: 149791
注册: 2005-12-29
性别: GG
专业: 污染控制化学
管辖:
计算模拟
【答案】应助回帖
★ ★
cenwanglai(金币+2): ~~ 2011-11-10 12:40:16
ywxing(金币+5): thanks 2011-11-10 12:49:14
>> syms x y a b
>> d=x+y+a+b
xdai=[1 3 5 7];
ydai=[2 3 5 6];
adai=3;
bdai=4;
d =
a + b + x + y
>> d=subs(d)
d =
a + b + x + y
>> d=subs(d,{'x','y','a','b'},{xdai, ydai,a,b})
d =
[ a + b + 3, a + b + 6, a + b + 10, a + b + 13]
>> d=subs(d,{'x','y','a','b'},{xdai, ydai,adai,bdai})
d =
10 13 17 20
赞
一下
(1人)
回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
2楼
2011-11-10 11:33:29
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
匿名
用户注销
(小有名气)
应助: 0
(幼儿园)
金币: 58.7
散金: 66
红花: 3
帖子: 141
在线: 95.5小时
虫号: 0
注册: 2010-09-21
性别: GG
专业: 固体与表面物理化学
本帖仅楼主可见
赞
一下
3楼
2011-11-10 12:51:44
已阅
申请计算强帖
回复此楼
编辑
查看我的主页
匿名
用户注销
(小有名气)
应助: 0
(幼儿园)
金币: 58.7
散金: 66
红花: 3
帖子: 141
在线: 95.5小时
虫号: 0
注册: 2010-09-21
性别: GG
专业: 固体与表面物理化学
本帖仅楼主可见
赞
一下
4楼
2011-11-10 13:36:10
已阅
申请计算强帖
回复此楼
编辑
查看我的主页
dbb627
荣誉版主
(著名写手)
专家经验: +4
计算强帖: 12
应助: 289
(大学生)
贵宾: 0.589
金币: 24640.4
散金: 551
红花: 61
沙发: 1
帖子: 1246
在线: 1794.8小时
虫号: 149791
注册: 2005-12-29
性别: GG
专业: 污染控制化学
管辖:
计算模拟
你的定义方式不对,要么所有的Gdr_1,Gdr_2相关的放到funTM中定义,或者定义Gdr_1,Gdr_2为全局变量,比较简单的事按照我下面的方式定义
匿名函数
CODE:
syms f_1 f_2 r_1 r_2
G=(f_1*f_2+r_2*f_2^2)/(r_1*f_1^2+2*f_1*f_2+r_2*f_2^2);
Gdr_1=f_1^2*(f_1*f_2 + r_2*f_2^2)/(2*f_1*f_2 + f_1^2*r_1 + f_2^2*r_2)^2;
Gdr_2=f_2^2*(f_1*f_2 + r_2*f_1^2)/(2*f_1*f_2 + f_1^2*r_1 + f_2^2*r_2)^2;
%定义符号表达式
f_11=[0.925 0.846 0.763 0.674 0.579 0.479 0.256 0.133];
f_22=[0.075 0.154 0.237 0.326 0.421 0.521 0.744 0.867];
F_22=[0.441 0.699 0.753 0.828 0.864 0.900 0.968 0.983];
r_11=0.09;r_22=10;
%准备给符号表达式赋值
Gdr_1=subs(Gdr_1,{'f_1','f_2','r_1','r_2'},{f_11,f_22,r_11,r_22});
Gdr_2=subs(Gdr_2,{'f_1','f_2','r_1','r_2'},{f_11,f_22,r_11,r_22});
G=subs(G,{'f_1','f_2','r_1','r_2'},{f_11,f_22,r_11,r_22});
d=F_22-G;
myfun=@(x,xdata)xdata(1,:)*x(1)+xdata(2,:)*x(2);
xdata=[Gdr_1;Gdr_2];
x0=[0.1 10];
x = lsqcurvefit(@(x,xdata) myfun(x,xdata),x0,xdata,d)
结果
Local minimum found.
Optimization completed because the size of the gradient is less than
the default value of the function tolerance.
x =
-0.0514 0.0967
赞
一下
回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
5楼
2011-11-10 14:27:53
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
相关版块跳转
第一性原理
量子化学
计算模拟
分子模拟
仿真模拟
程序语言
我要订阅楼主
ywxing
的主题更新
5
1/1
返回列表
如果回帖内容含有宣传信息,请如实选中。否则帐号将被全论坛禁言
普通表情
龙
兔
虎
猫
百度网盘
|
360云盘
|
千易网盘
|
华为网盘
在新窗口页面中打开自己喜欢的网盘网站,将文件上传后,然后将下载链接复制到帖子内容中就可以了。
信息提示
关闭
请填处理意见
关闭
确定