polyfit 与polyder 问题求助

zhang598

各位大侠：
用polyfit可以对矩阵中某一列与对应列进行拟合，然后用polyder求导数；
如a=

5	-2110.88	-6477.71
4.5	-2344.45	-6250.34
4	-2573.46	-6026.88
3.5	-2796.84	-5808.39
3	-3013.1	-5596.38
2.5	-3229.21	-5383.86
2	-3445.52	-5170.49
1.5	-3661.92	-4956.38
1	-3878.31	-4741.62
0.5	-4094.58	-4526.34
0	-4310.62	-4310.62
-0.5	-4526.34	-4094.58
-1	-4741.62	-3878.31
-1.5	-4956.38	-3661.92
-2	-5170.49	-3445.52
-2.5	-5383.86	-3229.21
-3	-5596.38	-3013.1
-3.5	-5808.39	-2796.84
-4	-6026.88	-2573.46
-4.5	-6250.34	-2344.45
-5	-6477.71	-2110.88

上述矩阵a，可以用

q1=polyfit(a(:,1),a(:,2),1);
q2=polyfit(a(:,1),a(:,3),1);
d1=polyder(q1);
d2=polyder(q2);
求得第二、三列对第一列的拟合关系及一阶导数，
但是，若在a1=（-5，-2）、（-2，2）、（2，5）三段内分别求第二、三列对第一列的拟合关系及导数，怎么实现呢？

感谢指导！

ljelly

II   你可以按段分成三个小矩阵，再按你上面的代码执行即可
   既然知道，他们在矩阵中相应的位置，那就好办
   如果不知道位置，可以加个判断，找出相应范围内的数据
   再进行拟合和微分即可

1. 
   
2.      ad1=find(a(:,1)<=-2&a(:,1)>=-5)
   
3.      ad2=find(a(:,1)<=2&a(:,1)>=-2)
   
4.      ad3=find(a(:,1)<=5&a(:,1)>=2)
   
5.      a1=a(ad1,:)
   
6.      a2=a(ad2,:)
   
7.      a3=a(ad3,:)
   
8.      q1=polyfit(a1(:,1),a1(:,2),1);
   
9.      q2=polyfit(a1(:,1),a1(:,3),1);
   
10.      d1=polyder(q1);
    
11.      d2=polyder(q2);
    
12.      q3=polyfit(a2(:,1),a2(:,2),1);
    
13.      q4=polyfit(a2(:,1),a2(:,3),1);
    
14.      d3=polyder(q3);
    
15.      d4=polyder(q4);
    
16.      q5=polyfit(a3(:,1),a3(:,2),1);
    
17.      q6=polyfit(a3(:,1),a3(:,3),1);
    
18.      d5=polyder(q5);
    
19.      d6=polyder(q6);
    
20. 
    
21.      

复制代码

zhang598

斑竹非常到位！
我的数据非常庞大，上述仅是举了个例子，凭肉眼找其位置的话极其不便，斑竹指出了利用判断语句来查找并赋值的方法和命令，方法极其到位！
帮我解决了一个大问题。
十分感激！

nwpukang

恩，不错，很强大