代码搜索结果

%Neville插值法 function y=Neville(x1,y1,x) m=length(x); n=length(x1); for k=1:m for i=2:n for j=2:i y1(i,j)=((x(k)-x1(i))*y1(i-1,j-1)-(x(k)-x1(i-j+1))*y1(i,j-1))/(x1(i-j+1)-x1(i));

close all; clear all; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%具体实例%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% a=-1;b=1;c=-1;d=1; n=6;m=5;TOL=1e-10; ITMAX=500; %f=inline('x*exp(y)','x','y'); %ga=inline('0','x

x=1:8; y=[15.3,20.5,27.4,36.6,49.1,65.6,87.8,117.6]; y1=log(y); n=1;%表示用一次函数拟和 m=length(x); A=zeros(m,n+1); for i=1:2 A(:,i)=x.^(i-1)'; end B=A'*A; C=A'*y1'; re=B\C; a=exp(re(1)); b=r

a=0;b=1;m=10;T=0.5;N=50;af=1; f=inline('sin(pi*x)','x'); h=(b-a)/m; k=T/N; lmd=af^2*k/h^2; x=linspace(a,b,m+1); x=x(2:m); i=1:m-1; u=f(i.*h); for j=1:N t=j*k; u=trisys(-lmd*ones(m-2

function fe=fitfun() %连续函数的最佳逼近 %取基{1, x, ...} %fitfun.m %默认算例为课本：P60，例3.1 %原函数f(x)=x^(1/2)， 定义域 [1/4, 1] %结果：P(x) = 10/27 + 88/135x 平方误差=0.00010803 fs=input('\n输入原函数f(x)：[直接回车表

close all; clear all; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%具体实例%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% a=0;b=2;c=0;d=1; n=6;m=5;TOL=1e-10; ITMAX=200; f=inline('x*exp(y)','x','y'); ga=inline('0','x','y

%Newton插值法 function y=Newton(x1,y1,x); m=length(x); n=length(x1); for i=2:n for j=n:-1:i y1(j)=(y1(j)-y1(j-1))/(x1(j)-x1(j-i+1)); end end for j=1:m p(j)=y1(1); for k

%x=[-1.0,-0.5,0.0,0.5,1.0,1.5,2.0]; %y=[-4.447,-0.425,0.551,0.048,-0.447,0.549,4.552]; %x=-3:3; %y=[4,2,3,0,-1,-2,-5]; x=[0.1,0.2,0.3,0.4,0.5]; y=[5.1234,5.3057,5.5687,5.9378,6.4370]; n=2; m=le

%Langrange插值法 function y=Langrange(x1,y1,x); m=length(x); n=length(x1); for i=1:m y(i)=0; for j=1:n l=1; for k=1:n if (k~=j) l=(x(i)-x1(k)

% ambfn1.m - plots ambiguity function of a signal u_basic (row vector) % % The m-file returns a plot of quadrants 1 and 2 of the ambiguity function of a signal % The ambiguity function is defined as