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的图。

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

谢谢。

回复此楼

» 猜你喜欢

读博已经有5人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有13人回复
博士申请都是内定的吗？已经有6人回复
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？已经有5人回复
博士读完未来一定会好吗已经有29人回复
投稿精细化工已经有4人回复
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！已经有4人回复
导师想让我从独立一作变成了共一第一已经有9人回复
心脉受损已经有5人回复
Springer期刊投稿求助已经有4人回复

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

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

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

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的回帖

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（ ...

那个简单的画图程序没写太严谨，那个没问题，问题是怎么在m文件里面加入diff的问题谢谢参与

赞一下

回复此楼

3楼2014-09-08 16:08:15

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