我今天正面对的关于线性代数的问题

PassByYou

前天策划区有人说我自大。非也。我一点都不自大。

因为我想策划的数学都是可以的，所以我抓了幅图，来咨询关于线性代数的算法。我不怕麻烦的，就怕这个问题解决不了，或则没找到最佳的算法公式。

我的最终目的是，根据4个来自角的母数，平滑一个区域。区域大小有3种情况，32*32和64*64，128*128，均为2的次方。为了分析简易可行，我将demo的数值用手写的，填在表格中。我的目的是根据4个来自角的母数，平滑整个区域。

该问题是我昨天遇到的，因为我用photoshop的柔边计算的线性代数有问题，结果出来的数字全是爆炸值，没起到平滑作用。

该问题就是对普通的线性代数，进行2次平滑插值，不知道该问题是否让各位策划感到有点为难？

GeeGee

这不用线性代数不一样能做么
先把4个边求出来
再一行一行的求就完事了……

GeeGee

#include <iostream>
using namespace std;

#define HOR 16
#define VER 16


//  a    i     b
//      
//  j  
//
//  d          c
void main(void)
{
        int a,b,c,d;
        int m[HOR][VER];
        int i,j;

        cout<<"please input a:"<<endl;
        cin>>a;
        cout<<"please input b:"<<endl;
        cin>>b;
        cout<<"please input c:"<<endl;
        cin>>c;
        cout<<"please input d:"<<endl;
        cin>>d;

       
        for(i = 0; i < HOR; i++)
        {
                m[0] = a + (b - a) * i / HOR;
                m[HOR - 1] = d + (c - d) * i / HOR;
        }
        for(j = 0; j < VER; j++)
        {
                m[0][j] = a + (d - a) * j / VER;
                m[VER - 1][j] = b + (c - b) * j / VER;
        }
        for(i = 1; i < HOR - 1; i++)
        {
                for(j = 1; j < VER - 1; j++)
                {
                        m[j] = m[0] + (m[VER - 1] - m[0]) * j / (VER - 2);
                }
        }
        for(i = 0; i < HOR; i++)
        {
                for(j = 0; j < VER; j++)
                {
                        cout<<m[j]<<" ";
                }
                cout<<endl;
        }       
}

PassByYou

geegee能不能解释一下这一句的计算思路？
m[j] = m[0] + (m[VER - 1] - m[0]) * j / (VER - 2);

PassByYou

如果当前计算坐标不是在4个边的话，会发生4个数交替，

geegee，请不能介绍一下，如果4个数交替，你的处理思路？

neelkey

仿佛跟线性代数关系不大

PassByYou

你只算出了单线，没有算出临接数
不过还是谢谢你，至少给了我不少启发。

关于临接的数值，应该是4个数中的最大数减掉，或则4个数的最小数加上，4个数的两个差的最大差的除2。geegee如果你把算法加以改进，考虑一下交叉值就对了。

看你是临时随手写的算法，思路很清晰，表示佩服。虽然你忘了必须round取整。在math库中。这个算法我已经明白了。多谢启发。

zfscnu

记得matlab中有类似的函数

以前做毕业设计的时候，处理图像的时候，需要对图像进行锐化，自己写过一个算法，但效率和效果都不太好，查看帮助，看到一个函数，使用后效率和效果都大大好过自己写的算法

matlab我还记得带有转为C的工具，也就是说在matlab写好的代码可以转为C代码

zfscnu

忘记说了，matlab的算法大多是基于矩阵的，也是楼主希望的“线性代数”算法^…^

GeeGee

不是啊，我是一行一行算的，因为9 x 7那行的开头和结尾本来就是算出来的，相邻2行之间的关系早在算开头和结尾的时候就已经确定了

所以根本不用看6 x 10那列对x的影响

另外int也不用取整，他自动就把尾数忽略掉了，这里面没必要四舍五入吧