不拟合怎么求用矩阵画出的图像的交点？

purple_paradise

不拟合怎么求用矩阵画出的图像的交点？
比如：
>> a=[1 3];b=[2 0];plot(a,'r');
>> grid on;hold on;plot(b,'g');

现在想求两条线的交点，不用拟合的方法，怎么求？
跪求啊

ch_j1985

plot(a,'r');
plot可以这么用吗？
如果可以，画出来的是什么图形？麻烦LZ贴出来

purple_paradise

可以，画矩阵啊，
但是我不知道怎么贴出来，找不到附件添加的

ch_j1985

发表新回复－>上传附件中的【浏览】－>选中图像，点击【插入】即可

以前没有用过这种方法画直线，故有点儿好奇。
电脑上没有装Matlab，还要烦请LZ贴图！

[ 本帖最后由 ChaChing 于 2010-4-3 07:24 编辑 ]

purple_paradise

贴不上来呀，要求：
文件尺寸: 小于 195.3125 kb
我怎么截都大于 195.3125 kb

friendchj

原帖由 ch_j1985 于 2008-11-26 21:42 发表

                               
登录/注册后可看大图

plot(a,'r');
plot可以这么用吗？
如果可以，画出来的是什么图形？麻烦LZ贴出来

plot(a,'r');应相当于plot(1:length(a),a,'r');
我记得如果a是矩阵的话，plot(a);应是对a的每一列画一条线
假设两个直线方程分别为：y=k1*x+b1;y=k2*x+b2，很容易求解方程组得到交点。但有多个点求解直线方程时，可能就需要拟合了，要不然就任取两点求直线方程参数，有时效果可能会差很多。

[ 本帖最后由 friendchj 于 2008-11-26 22:39 编辑 ]

ch_j1985

用附件中的画图工具或QQ截图，保存时注意格式（.jpg格式的图片比较小）

wwbeyondww

不明白为啥会想到用拟合来求交点。。。。。
当然用插值求，方便快捷啊
a=[1 3];b=[2 0];plot(a,'r');
grid on;hold on;plot(b,'g');
x=interp1(a-b,1:length(a),0);
y=interp1(1:length(a),a,x);

purple_paradise

有两个矩阵a_lim、n，都是aXb的，现在画了一个图，图上的点是（a_lim(i,j)，n(i,j)）;
for j=1:b
  plot(n(:,j),a_lim(:,j)); hold on;grid on;
end
这样每对应一个j,就有一条曲线，也就是有j条线，现在想求线的交点。怎么求？用插值怎么做？ 但是每一条都是独立的。
我现在就是要最下面的包络线，不要上面的，所以要求交点。但是不能通过拟合来求。

这里有7条线，最底下一个叶瓣对应一条，但是条条之间有很多的交叉

for i=1:b-1
    g(i)=interp1(qa_lim(:,i)-qa_lim(:,i+1),1:a,0)
    f(i)=interp1(1:a,qa_lim(:,i),g(i))
end
其中a、b为qa_lim的大小

??? Error using ==> interp1 at 261
The values of X should be distinct.

Error in ==> yeban_plot00014paper at 70
    g(i)=interp1(qa_lim(:,i)-qa_lim(:,i+1),1:a,0)

[ 本帖最后由 ChaChing 于 2010-4-3 07:30 编辑 ]