| 查看: 546 | 回复: 0 | ||
[求助]
如何matlab导出三次样条的方程？
|
|
p=[0 ,0.6283 , 1.2566 , 1.8850 , 2.5133 , 3.1416 , 3.7699 , 4.3982 5.0265, 5.6549 , 6.2832;0 , 0.1564 , 0.3090 , 0.4540 , 0.5878 , 0.7071 , 0.8090 , 0.8910 , 0.9511 , 0.9877 , 1.0000]; t=0:0.005:1; hold on for i=1:8 x=p(1,i)*(1/6)*(-t.^3+3*t.^2-3*t+1)+p(1,i+1)*(1/6)*(3*t.^3-6*t.^2+4)... +p(1,i+2)*(1/6)*(-3*t.^3+3*t.^2+3*t+1)+p(1,i+3)*(1/6)*t.^3; fprintf('x=%f*t.^3+%f*t.^2+%f*t+%f'); y=p(2,i)*(1/6)*(-t.^3+3*t.^2-3*t+1)+p(2,i+1)*(1/6)*(3*t.^3-6*t.^2+4)... +p(2,i+2)*(1/6)*(-3*t.^3+3*t.^2+3*t+1)+p(2,i+3)*(1/6)*t.^3; plot(x,y,'--r'); end 程序如上所示,但有错误，我想得到每次循环的x关于t的三次方程，手算比较麻烦，望各位大神指点。 |
回复此楼» 猜你喜欢
5个%2F，非固定段22个字母
已经有8人回复
-
面上的水很深
已经有15人回复
-
国社科通讯评审
已经有4人回复
-
生命口C13面上会评时间
已经有11人回复
-
内心匮乏
已经有16人回复
-
生命口面上会评结束了吗大家有知道结果吗？现在还没有结果正常不？
已经有9人回复
-
投英文期刊如何查重
已经有6人回复
-
不要再数国自然申请书的 filecode 的分隔符个数了
已经有21人回复
-
伯胺上甲基求助
已经有7人回复
-
跨出版社商投稿
已经有5人回复
10