题目
用下面代码建了一个function文件 function f=Language(x,y,x0) syms t; if(length(
用下面代码建了一个function文件
function f=Language(x,y,x0)
syms t;
if(length(x)=length(y))
n=length(x);
else
disp('x和y的维数不相等');
return;
end
f=0.0;
for(i=1:n)
l=y(i);
for(j=1:i-1)
l=l*(t-x(j))/(x(i)-x(j));
end;
for(j=i+1:n)
l=l*(t-x(j))/(x(i)-x(j));
end;
f=f+l;
simplify(f);
end
f0=subs(f,'t',x0);
然后在命令窗口输入如下三行
x=0:0.5:3;
>> y=[ 0 0.4794 0.8415 0.9975 0.9093 0.5985 0.1411];
>> [f,f0]=Language(x,y,1.6)
老是显示 Error:File:Language.m Line:3 Column:13
The expression to the left of the equals sign is not a valid target for an assignment.
请问各位大侠是哪里出错了!怎么修改?小弟新手!
用下面代码建了一个function文件
function f=Language(x,y,x0)
syms t;
if(length(x)=length(y))
n=length(x);
else
disp('x和y的维数不相等');
return;
end
f=0.0;
for(i=1:n)
l=y(i);
for(j=1:i-1)
l=l*(t-x(j))/(x(i)-x(j));
end;
for(j=i+1:n)
l=l*(t-x(j))/(x(i)-x(j));
end;
f=f+l;
simplify(f);
end
f0=subs(f,'t',x0);
然后在命令窗口输入如下三行
x=0:0.5:3;
>> y=[ 0 0.4794 0.8415 0.9975 0.9093 0.5985 0.1411];
>> [f,f0]=Language(x,y,1.6)
老是显示 Error:File:Language.m Line:3 Column:13
The expression to the left of the equals sign is not a valid target for an assignment.
请问各位大侠是哪里出错了!怎么修改?小弟新手!
提问时间:2021-03-03
答案
function [f,f0]=Language(x,y,x0)%%你要输出f和f0,那你的函数的输出参数个数必须是2个
syms t;
if(length(x)==length(y))%%这里是判断,不能用=,要用==
n=length(x);
else
disp('x和y的维数不相等');
return;
end
f=0.0;
for(i=1:n)
l=y(i);
for(j=1:i-1)
l=l*(t-x(j))/(x(i)-x(j));
end;
for(j=i+1:n)
l=l*(t-x(j))/(x(i)-x(j));
end;
f=f+l;
simplify(f);
end
f0=subs(f,'t',x0);
%%其它的不用改
>> x=0:0.5:3;
y=[ 0 0.4794 0.8415 0.9975 0.9093 0.5985 0.1411];
[f,f0]=Language(x,y,1.6)
f =
-799/3125*t*(t-1)*(t-3/2)*(t-2)*(t-5/2)*(t-3)+561/500*t*(t-1/2)*(t-3/2)*(t-2)*(t-5/2)*(t-3)-133/75*t*(t-1/2)*(t-1)*(t-2)*(t-5/2)*(t-3)+3031/2500*t*(t-1/2)*(t-1)*(t-3/2)*(t-5/2)*(t-3)-399/1250*t*(t-1/2)*(t-1)*(t-3/2)*(t-2)*(t-3)+1411/112500*t*(t-1/2)*(t-1)*(t-3/2)*(t-2)*(t-5/2)
f0 = 0.9996
%%%直接调用运行,上面是运行结果.
syms t;
if(length(x)==length(y))%%这里是判断,不能用=,要用==
n=length(x);
else
disp('x和y的维数不相等');
return;
end
f=0.0;
for(i=1:n)
l=y(i);
for(j=1:i-1)
l=l*(t-x(j))/(x(i)-x(j));
end;
for(j=i+1:n)
l=l*(t-x(j))/(x(i)-x(j));
end;
f=f+l;
simplify(f);
end
f0=subs(f,'t',x0);
%%其它的不用改
>> x=0:0.5:3;
y=[ 0 0.4794 0.8415 0.9975 0.9093 0.5985 0.1411];
[f,f0]=Language(x,y,1.6)
f =
-799/3125*t*(t-1)*(t-3/2)*(t-2)*(t-5/2)*(t-3)+561/500*t*(t-1/2)*(t-3/2)*(t-2)*(t-5/2)*(t-3)-133/75*t*(t-1/2)*(t-1)*(t-2)*(t-5/2)*(t-3)+3031/2500*t*(t-1/2)*(t-1)*(t-3/2)*(t-5/2)*(t-3)-399/1250*t*(t-1/2)*(t-1)*(t-3/2)*(t-2)*(t-3)+1411/112500*t*(t-1/2)*(t-1)*(t-3/2)*(t-2)*(t-5/2)
f0 = 0.9996
%%%直接调用运行,上面是运行结果.
举一反三
我想写一篇关于奥巴马的演讲的文章,写哪一篇好呢?为什么好
奥巴马演讲不用看稿子.为什么中国领导演讲要看?
想找英语初三上学期的首字母填空练习……
英语翻译
1,人们染上烟瘾,最终因吸烟使自己丧命.
最新试题
- 1要一个新的钟表怎么样?那些是旧的.这个是新的 英语翻译
- 2调查方式有哪几种?
- 3f(x)是以T为周期的函数,求f(x)+f(2x)+f(3x)+f(4x)的周期函数,
- 46666乘3333的计算方法
- 5有个数学题,算成本价的,
- 6沈石溪写的作品有哪些?
- 7二氧化碳__(难,易)溶于水
- 8园内收藏着无数的珍宝.圆明园被称为当时世界上最大的博物馆,艺术馆(用关联词把两句话和一句)
- 9一个长、宽、高分别为21、15、12(厘米)的长方体.现在从它的上面尽可能大的切下一个正方体.然后从剩余部分再尽可能大的切下一个正方体.最后再从第二次剩余的部分尽可能大的切下
- 101520千克的盐水中 含盐率为25% 要是这些盐水变为含盐率50%的盐水 需要蒸发多少千克的水 / 如果不蒸发水 而采用加盐的方法是它变为含盐率为50%的盐水 要加几千克的盐?(是两问 分别作答 要有
热门考点