| 查看: 272 | 回复: 0 | ||
womenyijia金虫 (小有名气)
|
[求助]
是一个矩阵的求解矩阵项数的扩展想用一个for循环编写，麻烦大家帮我看下
|
|
RGB信号到XYZ直接变换可以表示为H=PM,基于最小二乘法转换矩阵M=inv(P'*P)*P'*H 我这里用到了24组数据。先做了个3*3的转换。然后想把每个矩阵A中的三个值扩充为8个值[1 rgb rg rb gb r g b]最终矩阵P变为24*8的矩阵，最终可得到 一个8*3的转换矩阵M。 编一个for循环不想把这24组数据一个一个的带进去算。求大神帮忙 clc,clear,close all %RGB值 A1=[14.65,11.29,2.25];A2=[48.38,36.37,8.51];A3=[17.11,16.98,10.74]; A4=[17.11,16.98,10.74];A5=[25.74,22.17,13.44];A6=[30.71,37.08,15.02]; A7=[51.58,35.8,2.2];A8=[12.14,10.64,12.34];A9=[40.05,23.87,4.27]; A10=[10.52,7.37,4.33];A11=[41.96,43.61,4.56];A12=[59.05,46.23,2.67]; A13=[5.73,5.21,8.68];A14=[15.96,20.69,3.76];A15=[30.71,16.81,1.79]; A16=[73.15,63.51,3.72];A17=[40.56,23.99,9.57];A18=[11.82,15.46,12.55]; A19=[95.77,87.25,30.27];A20=[61.71,56.5,20.08];A21=[38.4,35.28,12.64]; A22=[20.64,18.86,6.7];A23=[9.87,9.13,3.34];A24=[3.94,3.6,1.34]; P=[A1;A2;A3;A4;A5;A6;A7;A8;A9;A10;A11;A12;A13;A14;A15;A16;A17;A18;A19;A20;A21;A22;A23;A24]; %XYZ值 B1=[110,49,13];B2=[219,126,76];B3=[97,97,93];B4=[76,75,16];B5=[146,105,106]; B6=[125,170,125];B7=[231,109,6];B8=[81,65,106];B9=[223,63,34];B10=[92,33,34]; B11=[172,170,42];B12=[230,145,4];B13=[40,42,88];B14=[73,125,41];B15=[204,39,12]; B16=[236,183,24];B17=[224,70,79];B18=[51,111,117];B19=[245,211,165]; B20=[216,184,138];B21=[117,146,102];B22=[126,96,61];B23=[69,52,29];B24=[31,24,12]; H=[B1;B2;B3;B4;B5;B6;B7;B8;B9;B10;B11;B12;B13;B14;B15;B16;B17;B18;B19;B20;B21;B22;B23;B24]; %求转换矩阵ｍ M=inv(P'*P)*P'*H %将rgb拓展为8项 for i=1:1:24 C(i)=[1,A(i)(1)*A(i)(2)*A(i)(3),A(i)(1)*A(i)(2),A(i)(1)*A(i)(3),A(i)(2)*A(i)(3),A(i)(1),A(i)(2),A(i)(3)]; end %求系数矩阵m M8=inv(C(i)'*C(i))*C(i)'*H |
回复此楼» 猜你喜欢
filecode
已经有10人回复
-
听说今天filecode变了
已经有3人回复
-
filecode
已经有12人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有5人回复
-
一个有机合成实验室都需要哪些设备？
已经有12人回复
-
求各位大神看下
已经有33人回复
-
咨询面上基金
已经有7人回复
-
影响面上的因素
已经有11人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有11人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有30人回复
找到一些相关的精华帖子，希望有用哦~
请大家帮我看看真空泵是哪里出了问题。
已经有7人回复
- 关于Matlab中for循环的问题，求达人指点，谢谢了
已经有7人回复
- 菜鸟求助关于fortran if循环的问题，经验丰富的前辈肯定一眼能看出问题出在哪
已经有3人回复
- 求助大家帮我看看这句话是什么意思，审稿人要求修改图片
已经有19人回复
- 一份JMCA的评审意见，请大家帮我看看如何回复，谢谢！
已经有4人回复
- 求助，L502错误，这是哪里出了问题，该怎么解决呢？帮忙看看。非常感谢
已经有12人回复
- 我做的一个紫外，不知出了什么问题，大家帮忙看看
已经有9人回复
- 我该怎么回复呢？求求大家帮忙看看
已经有21人回复
- 大家能帮我看一下这个SEM图像是怎么回事？
已经有14人回复
- 大家爱帮我看看这段MATLAB程序为啥出错了
已经有5人回复
- 一个gjf文件，总是pending,各位帮我看看是什么问题？
已经有16人回复
- matlab程序问题，总是出错，麻烦大家看看问题在那里？
已经有6人回复
- C++编程，关于循环结构的，大家看看我这程序哪里错了？
已经有21人回复
- 大家帮我看一个fortran的程序，我总是计算不出正确的结果
已经有20人回复
- 请大家帮我看看怎么回复这样的proof,在线等。急！
已经有5人回复
- 大家帮我看一下，核磁，有几个峰搞的不是很清楚
已经有40人回复
- 论文一审回来了，有一些审稿意见，请大家帮我看看是什么意思。
已经有5人回复
- 【求助】大家帮我看看我的循环伏安曲线，给个解释呗
已经有8人回复
- 几乎要崩溃了！大家帮我看看是否有挽回的余地？谢谢了
已经有112人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
50