谁有vbICA变分贝叶斯独立分量分析的matlab程序啊？哪位大神发一份啊！
答案:1 悬赏:70 手机版
解决时间 2021-03-11 00:05
- 提问者网友：疯孩纸
- 2021-03-10 02:49
谁有vbICA变分贝叶斯独立分量分析的matlab程序啊？哪位大神发一份啊！
最佳答案
- 五星知识达人网友：行路难
- 2021-03-10 03:28
用fsolve()函数求这个方程组的数值解比较好，具体代码如下：
1、建立方程组代码，保存为myfun.m文件
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
function F = myfun(x)
% 待求解的方程组，x为要求解的变量，是一个1×4的向量，前三个表示坐标，x(4)为辅助变量——因为n1和n2平行，x(4)即为二者的比值
A1=[0,0,0];
B=[2,1,1];
E=[1,1,0.5];
A2=x(1:3);
BA1=A1-B;
BA2=A2-B;
BE=E-B;
n1=cross(BA1,BA2);
n2 = cross(BA1,BE);
f = x(4)*n1-n2; % 因为n1和n2平行，则此式为零向量
F = [f(1);f(2);f(3);dot(BA1,BA2)/(norm(BA1)*norm(BA2))-cos(10/180*pi)];
2、在命令窗口直接求解：
>> x0=rand(1,4); % 任意假设一个初始解
>> fsolve(@myfun,x0)
Equation solved.
fsolve completed because the vector of function values is near zero
as measured by the default value of the function tolerance, and
the problem appears regular as measured by the gradient.
x =
0.1455 0.4794 0.0727 1.2296
这个题目的解好像不唯一。
1、建立方程组代码，保存为myfun.m文件
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
function F = myfun(x)
% 待求解的方程组，x为要求解的变量，是一个1×4的向量，前三个表示坐标，x(4)为辅助变量——因为n1和n2平行，x(4)即为二者的比值
A1=[0,0,0];
B=[2,1,1];
E=[1,1,0.5];
A2=x(1:3);
BA1=A1-B;
BA2=A2-B;
BE=E-B;
n1=cross(BA1,BA2);
n2 = cross(BA1,BE);
f = x(4)*n1-n2; % 因为n1和n2平行，则此式为零向量
F = [f(1);f(2);f(3);dot(BA1,BA2)/(norm(BA1)*norm(BA2))-cos(10/180*pi)];
2、在命令窗口直接求解：
>> x0=rand(1,4); % 任意假设一个初始解
>> fsolve(@myfun,x0)
Equation solved.
fsolve completed because the vector of function values is near zero
as measured by the default value of the function tolerance, and
the problem appears regular as measured by the gradient.
x =
0.1455 0.4794 0.0727 1.2296
这个题目的解好像不唯一。
我要举报
如以上回答内容为低俗、色情、不良、暴力、侵权、涉及违法等信息，可以点下面链接进行举报！
点此我要举报以上问答信息
大家都在看
推荐资讯