用matlab求解偏微分方程怎么求解1

温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!&&|&&
我喜欢不争取越众而出的刻意,我情愿站在原地,回头看看,等待水到渠成的相遇!
LOFTER精选
网易考拉推荐
用微信&&“扫一扫”
将文章分享到朋友圈。
用易信&&“扫一扫”
将文章分享到朋友圈。
阅读(6068)|
用微信&&“扫一扫”
将文章分享到朋友圈。
用易信&&“扫一扫”
将文章分享到朋友圈。
历史上的今天
在LOFTER的更多文章
loftPermalink:'',
id:'fks_081069',
blogTitle:'如何用MATLAB求一个函数的反函数',
blogAbstract:'&& syms n A L x&& f=L/(n^2*A)^(1/2)*atan((n^2*A)^(1/2)*x/(n^2-n^2*A*x^2-L^2)^(1/2)) &f =&L/(n^2*A)^(1/2)*atan((n^2*A)^(1/2)*x/(n^2-n^2*A*x^2-L^2)^(1/2))&&&& g=finverse(f)Warning: finverse(L/(n^2*A)^(1/2)*atan((n^2*A)^(1/2)*x/(n^2-n^2*A*x^2-L^2)^(1/2))) is not unique.& In sym.finverse at 43&g =',
blogTag:'',
blogUrl:'blog/static/',
isPublished:1,
istop:false,
modifyTime:0,
publishTime:0,
permalink:'blog/static/',
commentCount:0,
mainCommentCount:0,
recommendCount:0,
bsrk:-100,
publisherId:0,
recomBlogHome:false,
currentRecomBlog:false,
attachmentsFileIds:[],
groupInfo:{},
friendstatus:'none',
followstatus:'unFollow',
pubSucc:'',
visitorProvince:'',
visitorCity:'',
visitorNewUser:false,
postAddInfo:{},
mset:'000',
remindgoodnightblog:false,
isBlackVisitor:false,
isShowYodaoAd:false,
hostIntro:'我喜欢不争取越众而出的刻意,我情愿站在原地,回头看看,等待水到渠成的相遇!',
hmcon:'0',
selfRecomBlogCount:'0',
lofter_single:''
{list a as x}
{if x.moveFrom=='wap'}
{elseif x.moveFrom=='iphone'}
{elseif x.moveFrom=='android'}
{elseif x.moveFrom=='mobile'}
${a.selfIntro|escape}{if great260}${suplement}{/if}
{list a as x}
推荐过这篇日志的人:
{list a as x}
{if !!b&&b.length>0}
他们还推荐了:
{list b as y}
转载记录:
{list d as x}
{list a as x}
{list a as x}
{list a as x}
{list a as x}
{if x_index>4}{break}{/if}
${fn2(x.publishTime,'yyyy-MM-dd HH:mm:ss')}
{list a as x}
{if !!(blogDetail.preBlogPermalink)}
{if !!(blogDetail.nextBlogPermalink)}
{list a as x}
{if defined('newslist')&&newslist.length>0}
{list newslist as x}
{if x_index>7}{break}{/if}
{list a as x}
{var first_option =}
{list x.voteDetailList as voteToOption}
{if voteToOption==1}
{if first_option==false},{/if}&&“${b[voteToOption_index]}”&&
{if (x.role!="-1") },“我是${c[x.role]}”&&{/if}
&&&&&&&&${fn1(x.voteTime)}
{if x.userName==''}{/if}
网易公司版权所有&&
{list x.l as y}
{if defined('wl')}
{list wl as x}{/list}苹果/安卓/wp
请问老师们,如何使用matlab编程求出上述魔力值B,编写循环,学生不会求A的和的循环,谢谢啦。我自己搞的如下,但是不行*******************************%A为中间变量%Ti为第i个种子的生存时间, 即自种子发布起到现在所经过的时...间, 单位是周%T0为参数。T0 = 4%Si为第i个种子的大小,单位是GB5Ni为第i个种子当前的做种者数%N0为参数。N0 = 7%B为1小时中用户获得的做种魔力值点数5B0为参数,代表用户1小时获得魔力值的上限。B0 = 855L为参数。L = 400function b=moli(T0,Ti,N0,Ni,Si,B0,L)for i=1:Niai=ai+(1-10^(-T(i)/T0))*S(i)*(1+sqrt(2)*10^(-((N(i)-1)/(N0-1))endb=b0*2/pi*atan(ai/L)string=['您的魔力值为' num2str(b)]disp('string')***********************
载入中......
求救。自己顶。
细佬读过两年书,尘世中一介迷途小书童。
贴图被百度禁了
(24.05 KB)
22:42:27 上传
细佬读过两年书,尘世中一介迷途小书童。
看来没人解决了,只能把下一章看完继续。
细佬读过两年书,尘世中一介迷途小书童。
无限扩大经管职场人脉圈!每天抽选10位免费名额,现在就扫& 论坛VIP& 贵宾会员& 可免费加入
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
如有投资本站或合作意向,请联系(010-);
邮箱:service@pinggu.org
投诉或不良信息处理:(010-)
京ICP证090565号
论坛法律顾问:王进律师MATLAB帮助命令
本文所属图书&>&
本书结合高校数学课程教学和工程科学计算应用的需要,从实用角度出发,通过大量的算法实现,详尽、系统地介绍了MATLAB在高等数学问题求解中的应用。另外,为了帮助读者高效、直观地学习,作者对本书每章的重点内...&&
1.5& MATLAB帮助
MATLAB提供了强大的帮助,可以说学会使用MATLAB帮助系统,就意味着学会了自学MATLAB!用户通常在以下两种情况下需要使用MATLAB帮助系统。
在使用MATLAB的过程中遇到了一些问题,则可以查阅MATLAB帮助系统。
将MATLAB帮助系统作为学习资料来系统地学习MATLAB。
如果你具备一定的英语功底,那么建议你利用好MATLAB的帮助系统。MATLAB帮助系统主要提供了以下几种MATLAB帮助形式:帮助命令、帮助窗口、演示系统及远程帮助。
1.5.1& MATLAB帮助命令
MATLAB的所有命令、函数的M文件都有一个注释区。在该区中用纯文本形式简要地叙述了该函数的调用格式和输入、输出变量的含义。MATLAB主要提供了两个帮助命令:help和lookfor,这两个命令实现帮助的机制有所不同,下面分别予以介绍。
1. help命令
help命令通过在M文件的注释区找匹配的内容来显示帮助内容。例如,在命令行中执行如下语句可以得到函数help的所有文本内容。例如,已知magic.m函数的内容如图1-11所示。
然后在命令窗口中输入语句:
&& help magic
&MAGIC& Magic square.& % 函数文件中的在线帮助文档内容,即图1-11中矩形框中的内容
&&& MAGIC(N) is an N-by-N matrix constructed from the integers
&&& 1 through N^2 with equal row, column, and diagonal sums.
&&& Produces valid magic squares for all N & 0 except N = 2.
&&& Reference page in Help browser& % 以下内容提示进入帮助导航/界面
&&&&&& doc magic
有时我们不知道需要查询的函数的名称,这时可以根据该函数所属的主题来进行查询。可以在命令窗口中输入如下语句来查阅关于MATLAB系统中所有主题的帮助信息。
HELP topics:
Documents\MATLAB&&&&&&&&&& &- (No table of contents file)
matlab\general&&&&&&&&&&&& &- MATLAB通用命令
matlab\ops&&&&&&&&&&&&&&&&& &- 操作符和特殊字符
matlab\lang&&&&&&&&&&&&&&&& &- 语言结构
matlab\elmat&&&&&&&&&&&&&&& &- 初等矩阵和矩阵操作
matlab\randfun&&&&&&&&&&&&& &- 随机矩阵和随机流
matlab\elfun&&&&&&&&&&&&&&& &- 初等数学函数
matlab\specfun&&&&&&&&&&&&& &- 专用数学函数
matlab\matfun&&&&&&&&&&&&&& &- 矩阵函数-数值线性代数
matlab\datafun&&&&&&&&&&&&& &- 数据分析和傅里叶变换
matlab\polyfun&&&&&&&&&&&&& &- 插值与多项式
matlab\funfun&&&&&&&&&&&&&& &- 功能函数和微分方程求解器
&&&&&&&&&&&&&&&&&&&& &% 为节省篇幅,省略部分内容
rtw\accel&&&&&&&&&&&&&&&&&& &- (No table of contents file)
simulink\modeladvisor&&&& &- (No table of contents file)
stateflow\sfdemos&&&&&&&&& &- 状态流演示和采样
slvnv\rmidemos&&&&&&&&&&&&& &- (No table of contents file)
sam\studio&&&&&&&&&&&&&&&&&& &- (No table of contents file)
targets\autosar&&&&&&&&&&&& &- (No table of contents file)
symbolic\symbolic&&&&&&&&&& &- 符号数学工具箱
2. lookfor命令
lookfor命令用来在M函数文件中搜索包含某关键字的所有M函数文件名,例如,若要查询包含关键字jacobian的M文件,只需执行如下语句:
&& lookfor jacobian
numjac&&&&&&&&&&&&& &- Numerically compute the Jacobian dF/dY of function F(T,Y).
padess&&&&&&&&&&&& &- Sparse Pade approximation for use with Simulink/Jacobians
atamult&&&&&&&&&& &- Jacobian-matrix multiply
sfdnls&&&&&&&&&&& &- Sparse Jacobian via finite differences
gnns&&&&&&&&&&&&& &- Jacobian and error computation function for idgrey/idss
utEvalStateJacobian&- Gateway function for jacobian computation (dy/dx).
nlsf1&&&&&&&&&&&&& &- Nonlinear vector-valued function and Jacobian.
pdehypdf&&&&&&&&&& &- Jacobian for HYPERBOLIC.
pdeprbdf&&&&&&&&&& &- Jacobian for PARABOLIC.
csjac&&&&&&&&& &&- Compute Jacobian using complex-step differentiation.
说明:当知道某函数的函数名而不知道其用法时,help命令可以帮助用户准确了解此函数的用法。然而,当要查找一个不知其确切名称的函数名时,help命令就不能满足需要了。这种情况下,可以用lookfor命令来查询根据用户提供的关键字搜索到的相关函数。lookfor的查询机理为:它先对MATLAB搜索路径中的每个M文件的注释区的第一行进行扫描,一旦发现此行中含有所查询的字符串,则将该函数名及第一行注释全部显示在屏幕上。由此机理,用户可在自己的文件中加入在线注释。
您对本文章有什么意见或着疑问吗?请到您的关注和建议是我们前行的参考和动力&&
(window.slotbydup=window.slotbydup || []).push({
id: '2467141',
container: s,
size: '1000,90',
display: 'inlay-fix'
您的浏览器不支持嵌入式框架,或者当前配置为不显示嵌入式框架。
(window.slotbydup=window.slotbydup || []).push({
id: '2467142',
container: s,
size: '1000,90',
display: 'inlay-fix'
(window.slotbydup=window.slotbydup || []).push({
id: '2467143',
container: s,
size: '1000,90',
display: 'inlay-fix'
(window.slotbydup=window.slotbydup || []).push({
id: '2467148',
container: s,
size: '1000,90',
display: 'inlay-fix'查看: 56949|回复: 8|关注: 0
利用matlab函数怎么求一个函数的最大值呀
<h1 style="color:# 麦片财富积分
新手, 积分 5, 距离下一级还需 45 积分
现在有函数f(x)=4x^6-x-93,要求该函数在区间(0,100)上的最大值,怎么求呢?
<h1 style="color:# 麦片财富积分
关注者: 33
x=0:0.1:100;
y=4*x.^6-x-93;
[ymax,xmax]=max(y)复制代码
& && &&&1001
Matlab学无止境!
欢迎加入Matlab QQ交流群:
<h1 style="color:# 麦片财富积分
max(f(x))看看,找出f(x)的最大值。
论坛优秀回答者
关注者: 89
原帖由 cwjy 于
12:10 发表
x=0:0.1:100;
y=4*x.^6-x-93;
[ymax,xmax]=max(y)
& && &&&1001
这种情况下,未必是取大值,如果你把x再细分,最大值就又不一样了
这应该是函数求极值的问题
About me:/article-43-1.html
论坛优秀回答者
关注者: 89
原帖由 cwjy 于
12:10 发表
x=0:0.1:100;
y=4*x.^6-x-93;
[ymax,xmax]=max(y)
& && &&&1001
从函数表达式来看,其一阶导数为
y'=24*x^5 - 1
原函数在[0 0.5296]区间是单调递减,在[0.]上是单调递增的,最大值在x=100处
所以2楼的结果,从这点来看,值是对的
[ 本帖最后由 ljelly 于
12:31 编辑 ]
About me:/article-43-1.html
MATLAB 基础讨论板块优秀回答者
关注者: 99
原帖由 ljelly 于
12:27 PM 发表
从函数表达式来看,其一阶导数为
y'=24*x^5 - 1
原函数在[0 0.5296]区间是单调递减,在[0.]上是单调递增的,最大值在x=100处
所以2楼的结果,从这点来看,值是对的 ...
是的,最大值应该是在x=100处。
<h1 style="color:#5 麦片财富积分
关注者: 7
[x,fval,exitflag,output]=fminbnd('-(4*x^6-x-93)',0,100);
fval=-fval
fplot('4*x^6-x-93',[0,100])
<h1 style="color:# 麦片财富积分
关注者: 2
initpop=initializega(20,[0 100],'fitfit');
x=ga([0 100],'fitfit',[],initpop,[1e-25 1 0],'maxGenTerm',1000,'normGeomSelect',0.08,'heuristicXover',[20 3],'unifMutation',[4 0 0])
&&1.0e+012 *
& & 0.0000& & 4.0000
最大值为(1.0e+12)*4
function [sol,eval]=fitfit(sol,opt)
eval=4*x.^6-x-93;
<h1 style="color:# 麦片财富积分
Powered by}

我要回帖

更多关于 matlab求解偏微分方程 的文章

更多推荐

版权声明:文章内容来源于网络,版权归原作者所有,如有侵权请点击这里与我们联系,我们将及时删除。

点击添加站长微信