matlab中使用二分法对线性方程组求根

functiontestclearclcf0=@(x)x-cos(x);f2fen(f0,0,pi/2)%下面是二分法的函数部分functionx2=f2fen(f,a,b)x1=[a;b];p1=f

用矩阵表示需要拟合的数据,把大括号改为中括号即可>> X=[1,2,3,4,5,6,7,8,9];Y=[9,8,7,6,5,4,3,2,1];[P,S]=polyfit(X,Y,

输出是对的.因为他找到的是B(2:n,2)中的第一个元素.你可以好好看看find函数的作用.你如果写下面代码应该就对了：A=B(2:n,2);A(q=find(B(2:n,2)==max(abs((B

xn长a,hn长b,卷积结果yn长度为a+b-1如果结果的序列yn还要参与后续的运算,例如和xn相加,那么长度就不一致了,需要截断一般的xn或者hn非零的序列很短,两端的值为零或者趋于零,例如高斯函数

f=inline('x^2-x-2');>>[c,err,yc]=bisect(f,0,3,0.01)c=2.0010err=0.0059yc=0.0029-----------%使用二分法求解上面超

functiontesta=1;b=2;tol=1e-4;N=10000;k=0;fa=f(a);fork=1:Np=(a+b)/2;fp=f(p);if(fp==0||(b-a)/2再问：粘贴了您的

对(x,y)绘图,以x为横坐标,y为纵坐标,拆线图（默认）

这是源代码:在matlab中保存为:bisection.mfunction rtn=bisection(fx,xa,xb,n,delta)% 二分法解方程% fx是由方程

在Matlab里输入：x=[20,18,16,14,12,10,8,6,4,2,0,-2,-4,-6,-8,-10,-12,-14,-16,-18,-20];y=[-0.98106,-0.88756,

可以,只要在二分法外面,加一个循环就可以了.比如,以i为变量从1到1000,只要f(x(i-1))*f(xi)

程序如下：clear,clc;a=0;%a=input('inputa:');b=1;%b=input('inputb:');err=10^-5;y1=a*exp(a)

那调用两次不就ok了.区间你是用a、b定义的嘛.

你x10个值,y11个值,而且591.0也有误吧r=corrcoef(x,y);%r就是相关系数R=r^2;k=polyfit(x,y,1);scatter(x,y,'.');holdonx1=200

建议你在第三行插入一个end,以完成函数的定义而且函数的定义只有在m文件里面才能做到.（这很重要!）还有,好像没有与while相应的end对了,你的函数在[-pi/2,pi/2]里面只有一个解0,所以

就是把s与r排列到一起就相当于在r矩阵后面加了一个s矩阵扩大了r的列数

functiontestclearclc%实验方程：3*x.^2+x+2*exp(x)=0%原函数f=@(x)3*x.^2+x-2*exp(x);%导函数df=@(x)6*x+1-2*exp(x);%

A=[123-4-2;-3-4-12135;2100-310;4149-137];%A为系数矩阵fori=2:4forj=i:5fork=1:(i-1)A(i,j)=A(i,j)-A(i,k)*A(k

二、一元线性回归2．1．命令polyfit最小二乘多项式拟合[p,S]=polyfit（x,y,m）多项式y=a1xm+a2xm-1+…+amx+am+1其中x=（x1,x2,…,xm）x1…xm为（

这是由于matlab里面有个关于默认自变量规定问题一般情况下默认自变量为表达式中距离x最近的字母你可以用这种方法来求解[x0]=solve('2600*(z0+15)=(z0-24985)*x1')[

B=REGRESS(Y,X)