没落贵族 发表于 2010-12-9 13:38

同心圆的问题

用matlab做一个100个圆的同心圆,各个圆的颜色由给定的行列式的值来赋予,怎么做?
现在有行列式A,根据A的各个值从外圆开始给各个圆付上颜色,用来观察同心圆从内到外的光变化。请问有什么好方法吗?
A=[ 1.33E-41 2.14E-41 1.83E-41 8.16E-42 7.66E-42 2.26E-41 3.80E-41 3.78E-41 2.45E-41 1.14E-41 3.73E-42 1.95E-43 3.72E-42 1.54E-41 2.35E-41 1.72E-41 9.64E-42 2.60E-41 6.45E-41 9.07E-41 8.02E-41 4.80E-41 2.22E-41 9.34E-42 2.88E-42 6.31E-42 2.17E-41 3.08E-41 2.24E-41 2.84E-41 9.03E-41 1.86E-40 2.35E-40 1.99E-40 1.27E-40 8.37E-41 7.17E-41 6.28E-41 6.24E-41 8.37E-41 9.63E-41 8.39E-41 1.31E-40 3.45E-40 6.67E-40 8.79E-40 8.65E-40 7.65E-40 7.76E-40 9.01E-40 1.00E-39 1.05E-39 1.13E-39 1.18E-39 1.11E-39 1.07E-39 1.54E-39 2.72E-39 4.17E-39 5.21E-39 5.78E-39 6.51E-39 7.88E-39 9.63E-39 1.13E-38 1.27E-38 1.40E-38 1.48E-38 1.45E-38 1.28E-38 1.07E-38 9.57E-39 9.60E-39 9.91E-39 9.51E-39 8.27E-39 6.78E-39 5.60E-39 4.84E-39 4.23E-39 3.62E-39 3.19E-39 3.47E-39 4.74E-39 6.47E-39 7.45E-39 6.80E-39 4.86E-39 2.80E-39 1.55E-39 1.10E-39 1.03E-39 1.11E-39 1.60E-39 3.08E-39 6.45E-39 1.24E-38 2.03E-38 2.73E-38 3.02E-38];

如图,横轴代表圆的个数,纵轴代表颜色的变化。x=100的时候是圆心,请大家帮帮忙,谢谢

qibbxxt 发表于 2010-12-9 13:57

做个循环画,不就可以了吗
不同的圆传入不同的半径和颜色不就可以了吗?
另外,不明白A这个行列式是什么意思,行列式就是一个数,不知道该怎么用?

没落贵族 发表于 2010-12-9 14:15

A行列里的每一个数都代表颜色的变化。再问一下,怎么把不同的颜色由明到暗带入各个圆里??

leeking30 发表于 2010-12-9 20:18

本帖最后由 leeking30 于 2010-12-9 20:19 编辑

抛砖引玉 最简单的办法
clear;clc
A=[ 1.33E-41 2.14E-41 1.83E-41 8.16E-42 7.66E-42 2.26E-41 3.80E-41 3.78E-41 2.45E-41 1.14E-41 3.73E-42 1.95E-43 3.72E-42 1.54E-41 2.35E-41 1.72E-41 9.64E-42 2.60E-41 6.45E-41 9.07E-41 8.02E-41 4.80E-41 2.22E-41 9.34E-42 2.88E-42 6.31E-42 2.17E-41 3.08E-41 2.24E-41 2.84E-41 9.03E-41 1.86E-40 2.35E-40 1.99E-40 1.27E-40 8.37E-41 7.17E-41 6.28E-41 6.24E-41 8.37E-41 9.63E-41 8.39E-41 1.31E-40 3.45E-40 6.67E-40 8.79E-40 8.65E-40 7.65E-40 7.76E-40 9.01E-40 1.00E-39 1.05E-39 1.13E-39 1.18E-39 1.11E-39 1.07E-39 1.54E-39 2.72E-39 4.17E-39 5.21E-39 5.78E-39 6.51E-39 7.88E-39 9.63E-39 1.13E-38 1.27E-38 1.40E-38 1.48E-38 1.45E-38 1.28E-38 1.07E-38 9.57E-39 9.60E-39 9.91E-39 9.51E-39 8.27E-39 6.78E-39 5.60E-39 4.84E-39 4.23E-39 3.62E-39 3.19E-39 3.47E-39 4.74E-39 6.47E-39 7.45E-39 6.80E-39 4.86E-39 2.80E-39 1.55E-39 1.10E-39 1.03E-39 1.11E-39 1.60E-39 3.08E-39 6.45E-39 1.24E-38 2.03E-38 2.73E-38 3.02E-38];
A=A*1E38;
E=A/max(A);
for i=1:1:100
    radius=i;
    THETA=linspace(0,2*pi,1000);
    RHO=ones(1,1000)*radius;
    = pol2cart(THETA,RHO);
    Z=ones(1,1000)*A(i);
    plot3(X,Y,Z,'linewidth',1.5,'color',.* E(i) );hold on
end
view()

代码贴出来 不用下载费体能了

没落贵族 发表于 2010-12-10 15:46

谢谢leeking30了,非常感谢
页: [1]
查看完整版本: 同心圆的问题