请matlab高手帮我看一个程序,能否解释一下这些代码的含义?
来自: 更新日期:早些时候
MATLAB程序 求逐步解释程序中代码的意思~
请matlab高手帮我看一个程序,能否解释一下这些代码的含义?视频
相关评论:18957339578:请matlab高手帮我看一个程序,能否解释一下这些代码的含义?
戴才些接下来是Q矩阵中元素:Q中2到9行,各元素与b矩阵元素的关系通过这个式子来得到:然后是St、Sf这两个矩阵的求解,基本同上面的解释是一样的,只不过在此是从1到m的元素一步一步迭代,另外mean(Y )表示的是矩阵Y元素的平均值。Sqrt(St\/Sf) ,St\/Sf是开方。abs是绝对值 ...
18957339578:matlab高手帮我看看怎么把几个相交的曲面弄透明,汇出交线~~~
戴才些说明 1、生成网格数据的方法建议用linspace函数,而不要直接指定间隔(原来的数据间隔取的不太好)。2、为了得到更好的效果,网格数据需要稍微密集一些。3、求交线的做法是,在对应的x、y坐标,看两个平面(或曲面)的z值是否相差足够小。这种做法思路比较直接,但要求两个曲面的x、y坐标必需对应(如果...
18957339578:利用matlab,做三次样条差值,用的是外推样条,程序如文件,希望高手帮忙...
戴才些functionmymianx=[0123];y=[00.521.5];p=myspline(x,y)x=0:0.3:3;y=ppval(p,x);plot(x,y)functionp=myspline(x,y)%构造三次样条差值函数n=length(x);M=zeros(1,n);fork=1:n-1h(k)=x(k... function mymianx=[0 1 2 3];y=[0 0.5 2 1.5];p=myspline(x,y)x=0:0.3:3;y=ppval...
18957339578:matlab高手帮我看看怎样把这2个曲线最左边都放在圆点处。。。
戴才些做法其实就是把第一条曲线整体下移了初始的偏移值,程序如下:原程序第四行plot(r,y);换成plot(r,y-y(1));当然这就造成整条曲线的函数值都比实际值小,如果在保证函数值不变的前提下左端移到同一点应该无法实现。
18957339578:请高手帮我看一下我的matlab程序错在哪里了? max(min{f1,f2,f3,f4,f...
戴才些(2)怎么知道最小的最大值就一定为-7呢?(另,我这里试过,加入第4项约束时,最小最大值的确还是-7)帮你改写了一下非线性约束条件:function [C,ceq] = ghun(x)C = []; ceq = zeros(14,1);for i=1:7 ceq(2*i-1) = sum(x(1:7).^(i+1)) - sum((1:7).^(i+...
18957339578:请matlab高手帮忙分析一下这个程序。详细解释每一步什么意思!答案满意...
戴才些你透露的信息太少了,起码应该告知一下这个程序是用来干什么的,大概用的是什么算法啊。看来只能猜了,从最后画出的图来看,感觉这个程序的目的用逼近的方法得到一条分界线的样子,具体是什么算法,看不出来。你看看我的注解,自己研究一下吧。在网页里可能显得比较乱,你复制粘贴到matlab的编辑器了再看...
18957339578:我想用MATLAB做一幅正态分布图,但是做出的图怎么这么怪啊?请高手帮我...
戴才些首先,var是求方差,normpdf里面第三个参数sigma要求用标准差。所以把var函数换成std函数试试。还有你的数据可能比较集中吧,那这样也不奇怪的。你可以用减去平均值的数据画图试试
18957339578:matlab对一列数据进行fft变换 得不出频谱图 麻烦高手帮我看一下程序问 ...
戴才些fs=5000;N=length(t);n=0:N-1;tx=n\/fs;a1=fft(y1,N);a2= fftshift(a1); mag=abs(a2);%求出变换后的振幅f=(0:length(a2)-1)'*fs\/length(a2);subplot (2,1,1),plot(t,y1);subplot(2,1,2),plot(f,mag);title('频谱图')改成如上:...
18957339578:matlab编程,计算多值的y=e^x,请哪位高手帮我看看
戴才些clc;clear;x=linspace(0,5,10);ep(1,:)=ones(size(x));for i=1:5 ep(i+1,:)=x.^(i)\/factorial(i);end y=sum(ep);disp(y);plot(x,y,'*-');xlabel('x');ylabel('y');title('the function of y=e^x');legend('y=e^x');你的ep应当是一个矩阵....
18957339578:MATLAB高手帮一下忙啊.
戴才些clear;clc x=linspace(-2,8,100);p=[1 -4 0 3 -2 6];y=polyval(p,x);plot(x,y,x,0) %效果不太好,看不清楚 yy='x.^5-4*x.^4+3*x.^2-2*x+6';figure ezplot(yy,[-2 4.5])%这下子看清楚了 hold on,plot([-2:0.1:4.5],0)【2】xx=roots(p)结果:xx = 3...
1.matlab不支持中文的名字。
m文件可以取任何英文名字
2.randn(1,length(t))是生成一个随机分布的1×(t的长度)的正态分布的数组。
3.最后一个是生成P是一个Y.×Y的共轭矩阵。
要搞清楚什么意思,首先要知道里面的函数的含义,以后再遇到什么问题去查查matlab函数库查询字典吧
参考资料:函数库
一步一步来:
[m,n]=size(data) m,n即你自己的数据data这个矩阵的行数和列数,
Y=data(:,9) 是data这个矩阵的第九列构成的数组
X=zeros(38,9) 是由38行、9列构成的一个初始的0矩阵
X(:,1)=1表示X的第一列元素均为1
Z=zeros(38,1)、t、Q的解释同上面X的解释一样
然后是对X这个矩阵个元素的重新制定:
矩阵X的行数i为1到m列数j为2到9的元素,为data的i行,j-1列的元素
A为X的平方
C为A的逆矩阵
b 是C与X和Y的矩阵的乘积(注意矩阵相乘有行和列数的匹配)
Z矩阵的1到m行,个元素与b矩阵中元素的关系就是由
这个式子来获得。
接下来是Q矩阵中元素:Q中2到9行,各元素与b矩阵元素的关系通过这个式子来得到:
然后是St、Sf这两个矩阵的求解,基本同上面的解释是一样的,只不过在此是从1到m的元素一步一步迭代,另外mean(Y )表示的是矩阵Y元素的平均值。Sqrt(St/Sf) ,St/Sf是开方。abs是绝对值
请matlab高手帮我看一个程序,能否解释一下这些代码的含义?视频
相关评论:
戴才些接下来是Q矩阵中元素:Q中2到9行,各元素与b矩阵元素的关系通过这个式子来得到:然后是St、Sf这两个矩阵的求解,基本同上面的解释是一样的,只不过在此是从1到m的元素一步一步迭代,另外mean(Y )表示的是矩阵Y元素的平均值。Sqrt(St\/Sf) ,St\/Sf是开方。abs是绝对值 ...
戴才些说明 1、生成网格数据的方法建议用linspace函数,而不要直接指定间隔(原来的数据间隔取的不太好)。2、为了得到更好的效果,网格数据需要稍微密集一些。3、求交线的做法是,在对应的x、y坐标,看两个平面(或曲面)的z值是否相差足够小。这种做法思路比较直接,但要求两个曲面的x、y坐标必需对应(如果...
戴才些functionmymianx=[0123];y=[00.521.5];p=myspline(x,y)x=0:0.3:3;y=ppval(p,x);plot(x,y)functionp=myspline(x,y)%构造三次样条差值函数n=length(x);M=zeros(1,n);fork=1:n-1h(k)=x(k... function mymianx=[0 1 2 3];y=[0 0.5 2 1.5];p=myspline(x,y)x=0:0.3:3;y=ppval...
戴才些做法其实就是把第一条曲线整体下移了初始的偏移值,程序如下:原程序第四行plot(r,y);换成plot(r,y-y(1));当然这就造成整条曲线的函数值都比实际值小,如果在保证函数值不变的前提下左端移到同一点应该无法实现。
戴才些(2)怎么知道最小的最大值就一定为-7呢?(另,我这里试过,加入第4项约束时,最小最大值的确还是-7)帮你改写了一下非线性约束条件:function [C,ceq] = ghun(x)C = []; ceq = zeros(14,1);for i=1:7 ceq(2*i-1) = sum(x(1:7).^(i+1)) - sum((1:7).^(i+...
戴才些你透露的信息太少了,起码应该告知一下这个程序是用来干什么的,大概用的是什么算法啊。看来只能猜了,从最后画出的图来看,感觉这个程序的目的用逼近的方法得到一条分界线的样子,具体是什么算法,看不出来。你看看我的注解,自己研究一下吧。在网页里可能显得比较乱,你复制粘贴到matlab的编辑器了再看...
戴才些首先,var是求方差,normpdf里面第三个参数sigma要求用标准差。所以把var函数换成std函数试试。还有你的数据可能比较集中吧,那这样也不奇怪的。你可以用减去平均值的数据画图试试
戴才些fs=5000;N=length(t);n=0:N-1;tx=n\/fs;a1=fft(y1,N);a2= fftshift(a1); mag=abs(a2);%求出变换后的振幅f=(0:length(a2)-1)'*fs\/length(a2);subplot (2,1,1),plot(t,y1);subplot(2,1,2),plot(f,mag);title('频谱图')改成如上:...
戴才些clc;clear;x=linspace(0,5,10);ep(1,:)=ones(size(x));for i=1:5 ep(i+1,:)=x.^(i)\/factorial(i);end y=sum(ep);disp(y);plot(x,y,'*-');xlabel('x');ylabel('y');title('the function of y=e^x');legend('y=e^x');你的ep应当是一个矩阵....
戴才些clear;clc x=linspace(-2,8,100);p=[1 -4 0 3 -2 6];y=polyval(p,x);plot(x,y,x,0) %效果不太好,看不清楚 yy='x.^5-4*x.^4+3*x.^2-2*x+6';figure ezplot(yy,[-2 4.5])%这下子看清楚了 hold on,plot([-2:0.1:4.5],0)【2】xx=roots(p)结果:xx = 3...