24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

qmrain

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1016.5
散金: 190
帖子: 213
在线: 175.4小时
虫号: 823835
注册: 2009-08-07
性别: GG
专业: 凝聚态物性 II ：电子结构

[求助] 关于matlab求导和积分的两个问题已有1人参与

1、我写了一个m文件，比如是F=abc(x,y)，然后写了一个画图程序来画图F随着x或者y的变化。
下面我想在m文件中加上一步，即求F关于x或者y的导数，然后再利用画图程序调用m文件来画图比如
x=0:1：10；k=size(x,2)
for i=1:k
F=abc(x,y);
end
plot（x,F）

是不是直接可以在m文件最后一行加上 F=diff（F）；然后用画图程序来画图？这好像涉及到对x赋值过了，怎么再求导？

此外，好像还涉及到x和F最后矩阵对等的问题，总是出现 A（I）=B 有错误。

2、积分的问题：同上，F=abc(x,y)，这么一个m文件写好，我需要在程序中间引入另外一个变量z；但是这个z可以在m文件中最后一步积分积掉。

那么如何在m文件中引入z，最后顺利积掉，而不会出错，也不用改变F=abc（x,y）然后还可以画F关于x，y的图。

当然，上述好像都可以进行符号计算，然后再微分，积分，然后根据结果来画图。可是问题是有时候符号计算涉及到的量很多，最后形式非常复杂，不易得到。

谢谢。

回复此楼

» 猜你喜欢

急招9月入学博士，要有4级、最晚7月硕士毕业。精密电机驱控课题；学位材料已经有6人回复
有没有快的中文核心比较快录用的，纳米材料光催化已经有4人回复
本人42，博士刚毕业，现在找不到工作，怎么办？:( 已经有21人回复
河北省自然基金已经有6人回复
博士申请已经有5人回复
有人投过CCC中国控制会议吗？已经有3人回复
3,4-二羟基苯乙酮如何纯化？已经有5人回复
国基评审已经有10人回复
2026-博士申请已经有4人回复
考研调剂已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2014-09-07 21:42:41

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

qmrain

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1016.5
散金: 190
帖子: 213
在线: 175.4小时
虫号: 823835
注册: 2009-08-07
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

2楼: Originally posted by herowolf at 2014-09-08 00:23:14
x=0:1：10；k=size(x,2)
for i=1:k
F=abc(x,y);
end
plot（x,F）

你这样写有点问题。

你想画出x和F的变化的话，你得给定一个y值。
然后， F得是和x一样size的向量。

可以写成， F(i)。

F=diff（ ...

如果程序里面某一步涉及到求max或者min，怎么利用syms？

赞一下

回复此楼

4楼2014-09-08 16:16:41

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 6 个回答

herowolf

木虫 (正式写手)

应助: 25 (小学生)
金币: 3766.3
红花: 8
帖子: 731
在线: 381.1小时
虫号: 3283264
注册: 2014-06-19
性别: GG
专业: 流体力学

【答案】应助回帖

★
感谢参与，应助指数 +1
qmrain: 金币+1 2014-09-08 16:09:10

x=0:1：10；k=size(x,2)
for i=1:k
F=abc(x,y);
end
plot（x,F）

你这样写有点问题。

你想画出x和F的变化的话，你得给定一个y值。
然后， F得是和x一样size的向量。

可以写成， F(i)。

F=diff（F）这个好像也不行，好好查一查 diff命令。求导好像得在symbolic tool box 弄。
积分和求导我是一般在syms 里面弄的。

赞一下

回复此楼

2楼2014-09-08 00:23:14

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖