MATLAB m文件如何求w(i)的积分？已知w(i+1)=(1+T*(a.^2*k1+b.^2*k2/I*u))*w(i)...._百度知道
MATLAB m文件如何求w(i)的积分？已知w(i+1)=(1+T*(a.^2*k1+b.^2*k2/I*u))*w(i)....
y=[];v(i+1)=(1+T*(k1+k2)&#47,yr(1; n=length(t);V=30; I=4140;yr(1)=0; %期望轨迹的开始位置p=[],' %初始线速度xr=[]; %#车辆的横摆角速度; x=[];m*u(i))*v(i)+(T*(a*k1-b*k2)/2+v(i)*u(i)*sin(se(i)));for i=1:300);3).1:k&#39.44;).1*(u(i)*v(i)*cos(se(i))-(v(i).4; % #车辆横摆角度
x(i+1)=x(i)+0;yr=[]; % #车辆y位移
seo(i+1)=seo(i)+atan(2*(a+b)*(x(i)-xr(i))/m*u-u*T)*w(i)-T*k1*seo(i)/xr(1)=0;m; %#车辆的横摆角度se(1)=pi/%#车辆前轮转角vc=30;
f=trapz(T;r' plot(xr(1; ser=[];r=0.26;I*u;f=[].3; a=1.1;vc=[],w(i));xe=[];seo=[].^2);3.1*vc*cos(pi/ k1=80000; se=[]; plot(x(1.^2*k2/
yr(i+1)=yr(i)+0;v(1)=0;% #车辆v的计算w(i+1)=(1+T*(a:100.097;y(1)=0;x(1)=0;
xr(i+1)=xr(i)+0:300).1*vc*sin(pi/ T=0; w(1)=0;L;I*u))*w(i)-T*a*k1*seo(i)&#47，初始角速度v=[];b=1;I+T*(a*k1-b*k2)*v(i)/;2.1*(cos(se(i))*v(i),'3);ye=[].1*f.^2/ %实际轨迹的开始位置w=[];2);
se(i+1)=se(i)+0;k2=80000;% #车辆w的计算 end figure:n
L=V*T:300):300).^2*sin(se(i)))&#47,y(1; %#采样周期u=[]; % #车辆x位移
y(i+1)=y(i)+0t=0:0;
ser(i)=pi&#47.^2*k1+b
求w(i)的积分;I*u))*w(i)-T*a*k1*seo(i)&#47已知w(i+1)=(1+T*(a.^2*k1+b.^2*k2&#47，咋样求;I+T*(a*k1-b*k2)*v(i)/I*u
我有更好的答案
2*[t+C1+ln|sint+cost|+C2]=1/(sint+cost)dt=∫(cost+sint)/sintdt+t=∫d(sint)/(sint+cost)dt=∫1/(sint+cost)dt=∫d(sint+cost)/(sint+cost)dt则A+B=∫cost/(sint+cost)dt+∫sint&#47解;(sint+cost)=ln|sint+cost|+C2联立解得A=1&#47：令A=∫1/tant)+1]dt=∫(cott+1)dt=∫cost/(sint+cost)dt+∫sint/2*(t+ln|sint+cost|)+C副产品B=∫sint/(sint+cost)dt=∫dt=t+C1A-B=∫cost/(sint+cost)dtB=∫sint/(cott+1)dt=1/2*(t-ln|sint+cost|)+C如果是∫[(1/(tant+1)dt=∫cost/(sint+cost)dt=∫(cost-sint)&#47
有谁知道吗
其他类似问题
为您推荐：
等待您来回答
下载知道APP
随时随地咨询
出门在外也不愁matlab中M文件for循环求解题目意思是求 y(i)=(x(i+1)-x(i))/x(i)中的每个yfunction y(i)=rk(x)x=[3.9 5.3 7.2 9.6 12.9 17.1 13.2 31.4 38.6 50.2 62.9 76.0 92.0 106.5 123.2 131.7 150.7 179.3 204.5 226.5 251.4]for i=1:21y(i)=(x(i+1)-x(i))/x(i);endy运行结果出现这个 该怎么改 Undefined function or variable 'x'.
函数：function y=digui(x)n=length(x);for ii=1:n-1
y(ii)=(x(ii+1)-x(ii))/x(ii);end主程序：x=[3.9 5.3 7.2 9.6 12.9 17.1 13.2 31.4 38.6 50.2 62.9 76.0 92.0 106.5 123.2 131.7 150.7 179.3];y=digui(x)
步奏后面 能 加个说明吗？
n=length(x);这个是什么意思
一般情况下，定义函数时，没有必要赋值，当然有时也需要：
function y=digui(x)
定义函数digui.m
n=length(x);
% 向量x的长度为n
for ii=1:n-1
y(ii)=(x(ii+1)-x(ii))/x(ii);
x=[3.9 5.3 7.2 9.6 12.9 17.1 13.2 31.4 38.6 50.2 62.9 76.0 92.0 106.5 123.2 131.7 150.7 179.3];
y=digui(x)
为您推荐：
其他类似问题
扫描下载二维码matlab遗传算法程序出错 function result=sga(n,a,b,pc,pm,e)for i=1:50if (b-a)/e>2^(i)m=i+1;elsei=i+1;endendpopusize=n;chromlength=m;j=1;popu=round(rand(popusize,chromlength));while j
newpopu(i,:)=[popu(i,1:cpoint) popu(i+1,不能这样写,你是不是想把一行分成两行写?那么第一行末尾要加“...”的.
为您推荐：
其他类似问题
扫描下载二维码计算机二级考试题_百度知道
计算机二级考试题
}（9） 有以下程序段;n&quot,b,4.dat”; printf(&quot,20 C）11;
y2=fun( 【19】 ),21
C）2、填空题(每空2分,6:case ′I′：CCDAB
41-50;t=x[p1];
while(a：18;(k*(k+1));#include &lt,a;
i+=m+1.dat&}else if(x[m]}t=x[p0];n&
C）y=x:v1+=1; printf(&quot,p;!=b)
while(a&,2;j） { printf（&,i,*（z++））,b); p0=m.1f,&c），即将t所指字符串复制两次到s所指内存空间中, k=2;&#92,++*x,b=14,int n){ int p0。若运行时输入命令行,(37)以下程序的输出结果是#include void prt （int *x;a［p］）;3）{i+=2; p=&a［0］［0］;} main（） { int a=10;
程序的运行结果是A）CDEF B）ABEF
C）ABCD D）CDAB（43）有以下程序
#include &lt,u,12345;fun(while((_【16】_______);%d,&x),x[1][0]);m=i+x+y、选择题（（1）-（10）每题2分;=3,1.age).c=101;回车&}（13）以下程序统计从终端输入的字符中大写字母的个数,%2,3,&a,v1=4，否则可能导致数据丢失（50）有以下程序
#include &lt,int y)
{static int m=0.h&
f=add,int n） { int i,v2=7
D）v0=13,b=13,c=%d\ &#47,5
B）5,&b;%d&quot,2,42;%d\;n&%d&
i1=f(2),20
12;10,6),&c）;}
char a= ‘A’,m;
B）12, bbxxyyC)%lf&quot,22,3
B）2，若运行时从键盘输入,bbxxyy
D)} } main（） { int i=0，不会改变对应实参的值B）函数可以返回地址值C）可以给指针变量赋一个整数作为地址值D）当在程序的开头包含文件stdio,a+b);,7,ps&
printf(&func1（i）;n\} } void func2（int i） { printf（&
ps=&s[0],8;
{ prinf(&} main() {int a=1，正确的是 x=-1； %d%d&n\ /
main(int argc,c);
程序执行后的输出结果是
C）b D）a二;i++) printf(&quot: case ′e′，合法的一组C语言数值常量是A) 0x8A
C);k=N*N*5.10
则以下函数调用语句错误的是
A）k=(*f)(a;En&quot,y2,b=10;
int *p1=&a,c=%d＼n&quot,sin(x)
4; x=y=1,c=3,21 D）11; printf(&quot,k,b);i&;n&quot,使下面程序段的功能为完全相同s=1.0;10;%d&
printf(&}x,c=19 （19)下列程序段的输出结果是 int a=1234; B）x;,char *t){ char *p=t;i++){ n=n[i-1]*2+1;for(u=0：BDDBA
21-30;0;.h&quot,&b,c=19.0;
fun2(‘E’,&I++) scanf(&quot:#define N 100 则以下叙述中正确的是A)宏定义行中定义了标识符N的值为整数100B)在编译程序对C源程序进行预处理时用100替换标识符NC) 对C源程序进行编译时用100替换标识符ND)在运行时用100替换标识符N（49）以下叙述中正确的是
A）C语言中的文件是流式文件;
strcpy2(str1 , B）scanf(“%d”; } } while（c;n&%d %d＼n&quot,3; printf（&quot,j;;printf(&quot，从键盘输入 ADescriptor（表示回车）,i2); while（;;
fp=fopen(“abc,char *argv[])
{int i=1：proc。（6） 以下程序的输出结果是【6】;,b=3,++m; char st［］=&quot. 程序调试
6,%d& D）printf(“%d\A）-1
D）32(48) 若程序中有宏定义行.0&#47,b,8;for(u=0,2,123;
scanf(&quot。#include &lt:%d\***a=%d; } A）3
D）随机数 (36) 有以下程序#includevoid fun(int *a;
while(b&gt；后，主函数中将调用fun函数计算;
for(i=0;,int b){return+b},-2L;
y1=fun( 【18】 );%c&quot,s),1:case ′A′，则能正确执行的语句是A）scanf(&quot.177
D) 028（15）以下叙述中正确的是 A）构成C程序的基本单位是函数 B）可以在一个函数中定义另一个函数 C）main（）函数必须放在其他函数之前 D）所有被调用的函数一定要在调用之前进行定义 （16）若在定义语句,123;for(i=0;n&} while（;}（12）下列程序中的函数strcpy2()实现字符串两次复制, ‘F’),bbxxyy(39) 有以下程序 #include main（） { FILE *w&（46）若有以下定义和语句
union data
{}（14）以下程序中, j=10，调用strcpy2后;;%d& p=j;3;;{int t,*s++
13;}}（7）以下程序运行后的输出结果是 【7】 ;
printf(“%d %d\ D）scanf(“%d”,把最大的数与最后一个数对换，123;回车&gt:case ′E′;efgh&);;double fun(double x)
return(x*x-2*x+6); s［0］=k;stdio,5（29)有以下程序
fun(int x)
{int p：int a,其它依次类推;math.54321,7,i+' fclose（fp），判断i& printf(& for（i=j;
D）k=f(a.k&%d＼n& sum=fun((a++？,c);
k=fun(j;string：int a=25.队列是一种特殊的线性表;%d＼n&quot,aa),i2;#'.s--;,共40分)（1）在面向对象方法中,j; float b=123;I&v0=%d：1-10;26,z; } } A）4
D）6 （34）以下函数返回a所指数组中最大值所在的下标值 fun（int *a; }}（11）以下程序中函数f的功能是在数组x的n个数(假定n个数互不相同)中找出最大最小数; prt （&a,20,fun(7)),%d;}printf(&＼n&;i&lt.x[0]
12;%d\age);%d&#92,x,a;}程序的运行结果是A）732
D）874(23) 设有以下程序段 int x=0;%d&quot，a,将其中最小的数与第一个数对换; n/ 【13】______=t; printf(&
程序的运行结果是A）12
7（44）设有以下语句
typedef struct TT
{char c,st［i］）,v2=8
C）v0=7,m); fclose（fp），它的根结点数目是A）有且只有1
B）1或多于1
D）至少2（3）在面向对象方法中; float f,&quot,b); x[n-1]=t; void func2（int i）、填空题答案; x[p0]=x[n-1]!='for( y&
{FILE *fp,i,当执行p++,8;while(n,11
C）11,&p).共60分））（1）下列选项中不属于软件生命周期开发阶段任务的是
A）详细设计
B）概要设计
C）软件维护
D）软件测试（2）树是结点的集合; printf（&quot,(*f)(),v1=4;C)scanf(&quot,a;jif (a}main(){int c[10]={1;n&
while(i& fscanf（fp，类之间共享属性和操作的机制称为 【1】 ;;}（8）有以下程序, char b){printi(“%b%c”,1fp,*p=&c，则输出结果为A）2 0 0
D）2 4 6（41）有以下程序
int add(*fun函数的功能是将a所指数组元素从大到小排序*/\sex)); int v0=1;，若t所指字符串为efgh,b;;
if(x==0║x= =1) return(3),sum),1;
C）k=*f(a; fp=fopen（&
printf(&%d&void strcpy2(char *s; printf（&quot,n); } 在下划线处应填入的内容是 A）i=p
B）a［p］=a［i］
D）p=i （35）下列程序执行后的输出结果是 main（） { int a［3］［3］;%d＼n&
char sex,v1=%d. 3 7 5
B)0;之后.#include #include main(){ int num[26]={0},b=14,k, s［2］;;n”,%d&#92,str2)，接着执行以下选项中的语句,9;*给a赋值 b=20.4,3;=n
B)scanf(&}main(){char ss[10]=”acc”,k);j共执行的次数是 main（） { int i=0;
fseek(fp,7; i++）p［i］=i+1,&(ps&
k=fun(j:case ′U′,7:&while(_【9】_______){ s=s+1. 1
9;printf(“%s; if（i&lt,21,y2);scanf(&quot,str1),&quot,v1=4,s）;;
则以下语句正确的是A）x=10;%d& } s+=i!=0）s+=++x,4. 54321
7; printf(&quot,v2=%d＼n& _【10】_______,s); do{ switch（c=getchar（）） { case ′a′。例如、选择题答案,2; printf(& } A）11;
D) scanf(&quot，（11）-（50）每题1分; for （;%d,v2=0; char c,&n）;s=%f&#92,aa[10]=”bbxxyy”;%3d%3d%3d\;abcd&quot,v1,10);n&quot,41,0。main(){
int m=011;）,SEEK_END);3）{i+=2,i;i=j=x[0]，11&lt,b=13;
while( 【15】 =*p++).5
D）1234;k=1;);d1;
若运行时输入,2,b, int*z） { printf（&,”wb+”),&quot,t;}CIN,num[1]中统计字母B的个数，则以下语句的输出结果是
s=32;i++)Printf(&n”;%s\
}s[5]。请填空;
则下面叙述中正确的是
A）可以用TT定义结构体变量
B）TT是struct类型的变量C）可以用CIN定义结构体变量
D）CIN是struct TT类型的变量（45）有以下结构体说明,str2[]=&quot，因此只能顺序存取数据
B）打开一个已存在的文件并进行了写操作后：CABCD
BBBDB二;b) a-=b.5，实现信息隐蔽是依靠A）对象的继承
B）对象的多态
C）对象的封装 D）对象的分类（4）程序设计语言的基本成分是数据成分;for (i=0;,a;s=%f\ifor(j=i+1,22 (38) 有以下程序，则下面程序的运行结果是 #include main（） { char c,9:v0+=1,&quot, int*y;n&quot,p).h&gt。#include &Enter a,s）; A）是死循环
B）循环执行两次
C）循环执行一次
D）有语法错误 （28）有以下程序 fun(int x,a).5、运算成分,42;n&quot，下列说法错误的是
A）p向高地址移了一个字节
B）p抽高地址移了一个存储单元C）p向高地址移了两个字节
D）p与a+1等价（31）有以下程序
{int a=1;n”,11
D）11;,41; } 程序运行后的输出结果是 A）2;%d\ printf（&quot,&b),b;,a=5;%d&u++) printf(&
B）k=add(a. (x+8),c=19
C）### a=25，函数fun的功能是计算x2-2x+6!=0)t++,b;%3d%3d&#92，&k,y;!&; }
printf( “%c%c”：y1=(x+8)2-2(x+8)+6y2=sin2(x)-2sin(x)+6请填空,c).用#号结束输入,b),1（18) 若已定义; if（i&gt,5;;n&
则以下scanf函数调用语句中错误引用结构体变量成员的是
A）scanf(“%s”,z),c),7,5!x;
while(*s++=*t++);while((*t++=*s++).#include #define M 5#define N M+Mmain(){%c&i++，则程序输出结果是 【8】 ;}main(){ int a[10];printf(&quot,n)。 （4）线性表的存储结构主要分为顺序存储结构和链式存储结构,b);,v1=0，合并形成一个新的字符串,a; if（i&（17) 有以下程序 main() {int x!=0){ d=n%10,21.h时;%d&m{ if(x[m]&gt.（5） 【5】 的任务是诊断和改正程序中的错误，才能读到第1个数据D）当对文件的读（写）操作完成之后;i&n&quot. 元组
{int j=1,i1; }
程序运行后输出结果是 A）30
D）出错（21）若运行以下程序时:& printf（& prt （&a,k;k++) s=s+1，变量k中的值为 int k=3; p0=p1=0,x);4;}一,4,v1=4,v2）,3;2
D）n（n+1）&#47.#include main(){ int i.继承
2,9,c++), *p.#include viod f(int x[];
char ch[]=”abcd”,y1;
for(I=0;n”,b=%d：1;
执行后的输出结果是A）1
D）4（32）下列程序执行后的输出结果是 void func1（int i）,1234,9;y--)if(y%3= =0) printf(&quot,41 B）11;;d1;2（6）下列叙述中,++y，原有文件中的全部数据必定被覆盖
C）在一个程序中当对文件进行了写操作后;,n[4]={1},int n)&#47,x[0][1];,1;iif（a［i］&gt,c);，其中 【3】 反映了模块内各成分之间的联系；程序输出的结果是
A）***a=25，可以给指针变量赋NULL（27）以下程序段的描述,b=2,%d。（3）耦合和内聚是评价模块独立性的两个主要标准;n&a) b-=a;.dat&u&lt: case ′u′：#include void fun(char *t;; k=s［1］*10;i++)fwriter&ch[];%d&#92：DCBCA
31-40.0/请填空,c),5,a,3;n&,1.c=getchar()
i2=f(1+1)。 （2）一个关系表的行称为 【2】 ,*p=&a;）;u&lt,k);
void fun1( ){ a= ‘C’l b= ‘D’;}main(){ char str1[100]=&quot,且变量已正确定义和赋值for(s=1,0};
#define N 5
#define M N+1
#define f(x) (x*M)
{int i1.h&gt,5; printf(& return（p）.h&gt,d),i,请填空、控制成分和 A）语句成分
B）变量成分
C）对象成分
D）传输成分（5）冒泡排序在最坏情况下的比较次数是 A）n/ double c=12345,30};}
printf(“%d&#92.1f&
printf(&quot,1,num), 31
12;),t,%s\,n=0;以下三目运算符（.5 （20）有以下程序 main() {int a=0;,2,&s[0],k=1;for(i=1;%d;=10,b;
}main(){ A）不定值
D）20 （26）以下叙述中错误的是A）改变函数形参的值.#include main(){ int n=12345;r&
fread(&t;}
执行后的输出结果是A）7
D）0（30）在16位编译系统上; void func1（int i） { printf（&quot.h&gt：）所构成语句的执行后a&u++) scanf(&quot,4;n”;2
C）n（n-1）/ a=10;; int i=20,y2=%lf&#92,n）; fp=fopen（&)if(isupper(c)) num[c-‘A’]+= 【17】_______,fp);* 输出计算结果 *&#47,v2=12 (22) 有以下程序#includemain(){int y=9，可以作为软件开发工作的基础和依据 C）便于用户,j）;-;}
{ do {x=x*x?printf(&quot!=′＼n′）,*ps,s[0];
p=x-fun(x=2),x[i]);
printf(&quot,i,v0,s); int a[4]:case ′O′,*p=&c,b++b:printf(&quot.顺序 5;=25&&b--&lt,%d＼n&}程序运行的结果是A）1,y1，右有定义int a[ ]={10； fprintf（ }（15） 以下程序的输出结果是_【20】______;
printf(&quot.5,2; for（i=1，不需要操作系统的支持B）数据库技术的根本目标是要解决数据的共享问题C）数据库管理系统就是数据库系统D）以上三种说法都不对（10）下列叙述中正确的是 A）每个C程序文件中都必须有一个main()函数 B）在C程序的函数中不能定义另一个函数
C）C程序可以由一个或多个函数组成
D）在C程序中main()函数的位置是固定的（11）将E-R图转换到关系模式时;,6.name);10,int y){return(x+y);A' return m,a［1］［2］）;;=2&&c,b=14,&a);
printf(“%c\%c;C）scanf(“%c”, s=0: case ′o′,4,8,n,c=19
D）### a=26;func2（i）;
该程序生成的可执行文件名为;fun(c+4;;（41）有以下程序
#include &x[p1]= _【12】_______.456：246&lt：CACDC
11-20,6;arge){n=n+strlen(argv[i]),num[0]中统计字母A的个数;;
fclose(fp);%d&
s= 【14】 ,s=0,a，不属于软件需求规格说明书的作用的是A）作为确认测试和验收的依据B）反映出用户问题的结构,j=30;%d&quot,C) 0，实体与联系都可以表示成 A）属性
D）域 （12）下列定义变量的语句中错误的是
A）double int_,5, bbxxyy
B) acc,12345,d; func1（i）,b),for(m=0,char *s){while(*t,0,b= ‘B’.5
C）1234;(k*(k+1))、开发人员进行理解和交流D）便于开发人员进行需求分析 （7）下列叙述中错误的是
A）一个C语言程序只能实现一种算法 B）C程序可以由多个程序文件组成 C）C程序可以由一个或多个函数组成 D）一个函数可单独作为一个C程序文件存在（8）下列不属于软件工程的3个要素的是A）方法
D）环境 （9）下列叙述中正确的是A）数据库系统是一个独立的系统,i）;f(a; 给b赋值 *&#47,++*y; i&v2+=1!x）;a+b=%d\} A）hello
D）hlm （33）在下述程序中,&quot,
*p=*p1*(*p2),6;y1=%lf,c;),a+b); printf(&k&lt，必须将它关闭;i) {i=x[m],j=0,3,40
11: case ′i′; default!=0),6,*p2=&b2010年9月全国计算机等级考试二级C语言笔试试题及答案一; }
A）v0=11,b=0,aa); printf（&quot,1
D）2,n++),v2=11
B）v0=8;%d%d%d&}程序运行结果是A) accxyy ,42;） { i+=k;9; s^=32;）,c=20，必须先关闭该文件然后再打开;%d& 则 A）运行程序段后输出0
B）运行程序段后输出1 C）程序段中的控制表达式是非法的
D）程序段执行无限次 （24）有以下程序 fun(int x,I,t);%d,y,i=2.exe,m=I,--y);（47）若变量已正确定义,y++; z=x++,x[0][0],%2;printf(&（10）执行以下程序的输出结果是【11】______：
proc 123 45 67
则程序的输出结果是A）3 B）5
C）7 D）11（42）有以下程序
void fun2(n&quot,9,c=19
B）***a=26,m),b=40,b=%d.请填空;fprintf（###a=%d,c）; } 执行后的输出结果是
D）3（25）执行下面的程序段后;
C）char For
D）float USS（13）下列叙述中正确的是 A）C语言编译时不检查语法
B）C语言的子程序有过程和函数两种 C）C语言的函数可以嵌套定义 D）C语言的函数可以嵌套调用（14）以下选项中,循环队列是队列的【4】存储结构;,st［i］）;;。main(){,friend，s所指字符串为n&n&quot,b),k)
执行后的输出结果是A）5;i&lt,n=11; } 程序运行后的输出结果是 A）20 30
D）30 20 （40) 有以下程序
{int x[3][2]={0},p1、变量定义和赋值语句
struct STD
{char name[10];;%d
%d\for(i=0; B）int _printf(“%d”
其他类似问题
为您推荐：
您可能关注的推广
计算机二级考试的相关知识
等待您来回答
下载知道APP
随时随地咨询
出门在外也不愁matlab样条插值函数,已经有求插值函数的程序,想知道求一系列值的程序该怎么写function [sx]=cubic(x,y,dy0,dyn)n=length(x)-1;d=zeros(n+1,1);h=zeros(1,n-1);f1=zeros(1,n-1);f2=zeros(1,n-2);for i=1:nh(i)=x(i+1)-x(i);f1(i)=(y(i+1)-y(i))/h(i);endfor i=2:nf2(i)=(f1(i)-f1(i-1))/(x(i+1)-x(i-1));d(i)=6*f2(i);endd(1)=6*(f1(1)-dy0)/h(1);d(n+1)=6*(dyn-f1(n-1))/h(n-1);A=zeros(n+1,n+1);B=zeros(1,n-1);C=zeros(1,n-1);for i=1:n-1B(i)=h(i)/(h(i)+h(i+1));C(i)=1-B(i);endA(1,2)=1;A(n+1,n)=1;for i=1:n+1A(i,i)=2;endfor i=2:nA(i,i-1)=B(i-1);A(i,i+1)=C(i-1);endM=A\d;for i=1:nsx(i)=collect(y(i)+(f1(i)-(M(i)/3+M(i+1)/6)*h(i))*(t-x(i))+M(i)/2*(t-x(i))^2+(M(i+1)-M(i))/(6*h(i))*(t-x(i))^3);digits(4);sx(i)=vpa(sx(i));endfor i=1:ndisp('S(x)=');fprintf(' %s (%d,%d)\n',char(sx(i)),x(i),x(i+1));end这上面是求样条插值函数的程序,现在想求-1:1:100,这100个插值点函数的值,怎么运行,我编写的程序是S=zeros(1,99);for i=1:99t=-1+0.02*i;for j=1:nif (x(j)
报什么错啊?插值点时,你写的也不完整啊.你看看私信吧.
为您推荐：
其他类似问题
扫描下载二维码