MATLAB解方程与函数极值_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
文档贡献者贡献于
评价文档：
16页免费21页免费52页免费28页免费28页免费 28页免费31页免费27页免费28页免费28页免费
喜欢此文档的还喜欢8页免费2页免费41页免费7页免费6页免费
MATLAB解方程与函数极值|
把文档贴到Blog、BBS或个人站等：
普通尺寸(450*500pix)
较大尺寸(630*500pix)
大小：134.50KB
登录百度文库，专享文档复制特权，财富值每天免费拿！
你可能喜欢MATLAB软件与基础数学实验实验,和,数学,数学实验,和数学软件
扫扫二维码，随身浏览文档
手机或平板扫扫即可继续访问
MATLAB软件与基础数学实验
举报该文档为侵权文档。
举报该文档含有违规或不良信息。
反馈该文档无法正常浏览。
举报该文档为重复文档。
推荐理由：
将文档分享至：
分享完整地址
文档地址：
粘贴到BBS或博客
flash地址：
支持嵌入FLASH地址的网站使用
html代码：
&embed src='/DocinViewer-4.swf' width='100%' height='600' type=application/x-shockwave-flash ALLOWFULLSCREEN='true' ALLOWSCRIPTACCESS='always'&&/embed&
450px*300px480px*400px650px*490px
支持嵌入HTML代码的网站使用
您的内容已经提交成功
您所提交的内容需要审核后才能发布，请您等待！
3秒自动关闭窗口4:45:12【 转载互联网】 作者： &&|&责编：李强
&&& &为了解决用户可能碰到关于"一段画牛顿迭代收敛域的matlab程序，求解释"相关的问题，突袭网经过收集整理为用户提供相关的解决办法，请注意，解决办法仅供参考，不代表本网同意其意见,如有任何问题请与本网联系。"一段画牛顿迭代收敛域的matlab程序，求解释"相关的详细问题如下:d=6;tol=1e-5;maxIter=100;r=-2:0.01:2;
%实部虚部的范围[x y]=meshgrid(r);
%产生实部虚部二维网格Z=x+1i*y;
%Z对应网格的虚平面%% Define fuctionf=@(x,d) (x.^d)-1;fprime=@(x,d) d*(x.^(d-1));%% Perform Newton iterationsfor k=1:maxI
Z=Z-(f(Z,d)./fprime(Z,d));end%% Find d roots of unity, and the
maskrenderMat=0;for j=1:d
root=exp(2*pi*1i/d)^j;
% the jth root
Mj=abs(Z-root);
Z中每点都这个根的距离
% Each root gets a unique number in [1,d]
mask=(Mj&=tol)*j;
%Mj&=tol返回满足误差的逻辑矩阵
%满足误差部分为1*j，不满足部分为0
renderMat=renderMat+
%加起来之后renderMat中收敛于第j个根的区域数据都是j
%那么收敛于第j个根的区域都是同一种颜色endcolormap(hsv(d+1));
% Set the color mapimagesc(r,r,renderMat) % Render the fractalxlabel('Re(Z)');ylabel('Im(Z)');h=set(h,'ytick',(2*(0:d)+1)*d/(d+1)/2);str=arrayfun(@(x)num2str(x,'%.2f'),exp(2*pi*1i/d).^(1:d),'uniformoutput',false);set(h,'yticklabel',[{'未收敛'},str]);这是大神给的解答，基本看懂了，但是这个算法还是没看明白是怎样确定哪些区域是趋向于某个根的？root=exp(2*pi*1i/d)^j; 这一步什么意思，为什么冒出个exp（2pi*1i/d)^j什么意思?===========突袭网收集的解决方案如下===========
解决方案1：1的开d次方有d各根，第一个根是exp(2*pi*1i/d)其余的根分辨是这个根的 2，3，4，....，d倍所以root=exp(2*pi*1i/d)^j，当j取1时就是第一个根，j取2时就是第二个根，......
================可能对您有帮助================
问：这是迭代函数x=x-(H*sqrt(1+2*x/H)-x*arch(H/x+1)-b)/(1/sqrt(1+2*x/H)-...答：x=x-(H*sqrt(1+2*x/H)-x*arch(H/x+1)-b)/(1/sqrt(1+2*x/H)-arch(H/x+1))-x*((H/x^2+(H^2/x^3+H/x^2)/sqrt(H^2/x^2+2*H/x))/(H/x+1+sqrt(H^2/x^2+2*H/x))) arch后面少了一个)===========================================问：（1）牛顿迭代法算法步骤： 给定初始值 ， 为根的容许误差， 为 的容许...答：syms x f=x^x-10; df=diff(f,x); eps=1e-6; x0=10; cnt=0; MAXCNT=200; %最大循环次数 while cnt===========================================问：（1）牛顿迭代法算法步骤： 给定初始值 ， 为根的容许误差， 为 的容许...答：% 用牛顿下山法求解方程 function [x,k]=myfun_newton(f,x0,emg) % f表示非线形方程 % x0迭代初值，此种方法是局部收敛，初值要选择恰当 % emg是精度指标 % k,u分别表示迭代次数和下山因子 % d1表示非线形方程f在x0处的导数值 [f1,d1]=feval(f,x...===========================================问：求12-3x+2cosx=0在【2,4】内的解，Matlab程序如下： %清除所有内...答：建立.m文件直接运行试过没？===========================================问：要一个MATLAB用牛顿迭代求解非线性方程的程序程序，一定要能迭代的，最...答：给你一个完整版：% 牛顿法求解非线性方程组 function main() f = @(x)log(x+sin(x)); % 测试函数 df = @(x)(1+cos(x))/(x+sin(x)); % 导数函数 x0 = 0.1; % 迭代初值 x = TestNewton(f, df, x0) % 牛顿法求解 function x = Test...===========================================问：function [x]=Newtondd(f,x0,x1,e) j=1; flag=0; for n=10^j for i=0:n ...答：我也认为是 if abs(x_a-x_b)& 这一句可能一直被执行，所以没有给出x值而报错。而else 没有一次命中。待会儿再仔细看看===========================================问：function [x]=Newtondd(f,x0,x1,e) j=1; flag=0; for n=10^j for i=0:n ...答：例如求解方程组:x1^2+x2^2=4;x1^2-x2^2=1 可以编写程序如下： x0=1;y0=1;n=10;tol=1e-6; x(1)=x0;y(1)=y0; i=1;u=[1 1];k(1)=1; while(norm(u)&tol*norm([x(i),y(i)]')) A=2*[x(i),y(i);x(i),-y(i)]; b=[4-x(i)^2-y(i)^2,1-x(i)^2+y(i)^2...===========================================问：用5种颜色画出x^5-1=0的收敛域. 不怎么会matlab,老师布置的作业就硬着头...答：首先，你把a附错值了，应该是a=[1,0,0,0,0,-1] 其次，你把牛顿收敛的精度定得太高了，我i3-2100的电脑跑了2分多钟才把图给跑出来。 程序有没结束运行，看命令窗口的的引导符有没显示就可以了。估计是你的电脑太慢，等不到图出现。 PS.在matlab中...===========================================问：小弟的matlab学的不好,数值计算中的两个题关于牛顿迭代法的，急,谢谢!!!...答：源程序： %%牛顿迭代法求解非线性方程 function x=nanewton(fname,dfname,x0,e,N) %%fname,dfname分别表示f(x)及其导数的函数句柄或者内嵌函数，N为迭代次数上限 if nargin===========================================2012MATLAB培训2012
扫扫二维码，随身浏览文档
手机或平板扫扫即可继续访问
2012MATLAB培训
举报该文档为侵权文档。
举报该文档含有违规或不良信息。
反馈该文档无法正常浏览。
举报该文档为重复文档。
推荐理由：
将文档分享至：
分享完整地址
文档地址：
粘贴到BBS或博客
flash地址：
支持嵌入FLASH地址的网站使用
html代码：
&embed src='/DocinViewer-4.swf' width='100%' height='600' type=application/x-shockwave-flash ALLOWFULLSCREEN='true' ALLOWSCRIPTACCESS='always'&&/embed&
450px*300px480px*400px650px*490px
支持嵌入HTML代码的网站使用
您的内容已经提交成功
您所提交的内容需要审核后才能发布，请您等待！
3秒自动关闭窗口matlab中用函数句柄定义了一个函数func1=@(x,y)(x-10)*x+xy+y^2-8,如果求解方程func1=0，用solve或fsolve_百度知道
matlab中用函数句柄定义了一个函数func1=@(x,y)(x-10)*x+xy+y^2-8,如果求解方程func1=0，用solve或fsolve
solve(eq)中的eq要求是字符串或syms，函数句柄无法用。我的目的是给定x=[0:0.01:2]，如果求出对应的y值，希望有简介高效的办法，谢谢！
我有更好的答案
dy=zeros(2,1); dy(1)=y(2); dy(2)=2*x*y(2)/(1+x^2); ode45中要取函数句柄， [x,y]=ode(@sdre,[0 10],[1 3]); 我已
没说完嘛，y(2)是什么意思？报错啊
好专业的问题，深奥得很
其他类似问题
等待您来回答
下载知道APP
随时随地咨询
出门在外也不愁