clearclcr=0.1;s=0.01;T=0:0.5:20;w=0;f0=(-tanh(r)^0./cosh(r))
来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/05/25 20:48:15
clear
clc
r=0.1;
s=0.01;
T=0:0.5:20;
w=0;
f0=(-tanh(r)^0./cosh(r));
f1=(-tanh(r)^1./cosh(r));
B0=(f0./4)*s*sqrt((0+1)*(0+2)).*(1-cos(2.*T/s));
B1=(f1./4)*s*sqrt((0+2)*(0+3)).*(1-cos(2.*T/s));
C0=(f0./2)*s*(0+1)*sin(2.*T/s);
C1=(f1./2)*s*(0+2)*sin(2.*T/s);
w=0.5*(0*(f0^2)+(0+1).*(C0^2)+(2*0+2).*(B0^2)+(0+1).*f0.*f1+C0.*C1.*(0+2)+2*sqrt((0+1)*(0+3)).*B0.*B1);
for n=1:1000;
f(n)=((-tanh(r))^n)./cosh(r);
f(n+1)=((-tanh(r))^(n+1))./cosh(r);
B(n)=(f(n)./4)*s*sqrt((n+1)*(n+2)).*(1-cos(2.*T/s));
B(n+1)=(f(n+1)./4)*s*sqrt((n+2)*(n+3)).*(1-cos(2.*T/s));
C(n)=(f(n)./2)*s*(n+1)*sin(2.*T/s);
C(n+1)=(f(n+1)./2)*s*(n+2)*sin(2.*T/s);
W(n)=0.5*(n.*(f(n)^2)+(n+1).*(C(n)^2)+(2*n+2).*(B(n)^2)+(n+1).*f(n).*f(n+1)+C(n).*C(n+1).*(n+2)+2*sqrt((n+1)*(n+3)).*B(n).*B(n+1));
w=w+W(n)
end
T=0:0.5:20;
plot(T,w)
这个程序的主要目的就时要做出w随T变化的图,w是一个含有T的求和公式,n从0到1000求和 ,0《=T《=20,步长0.5
clc
r=0.1;
s=0.01;
T=0:0.5:20;
w=0;
f0=(-tanh(r)^0./cosh(r));
f1=(-tanh(r)^1./cosh(r));
B0=(f0./4)*s*sqrt((0+1)*(0+2)).*(1-cos(2.*T/s));
B1=(f1./4)*s*sqrt((0+2)*(0+3)).*(1-cos(2.*T/s));
C0=(f0./2)*s*(0+1)*sin(2.*T/s);
C1=(f1./2)*s*(0+2)*sin(2.*T/s);
w=0.5*(0*(f0^2)+(0+1).*(C0^2)+(2*0+2).*(B0^2)+(0+1).*f0.*f1+C0.*C1.*(0+2)+2*sqrt((0+1)*(0+3)).*B0.*B1);
for n=1:1000;
f(n)=((-tanh(r))^n)./cosh(r);
f(n+1)=((-tanh(r))^(n+1))./cosh(r);
B(n)=(f(n)./4)*s*sqrt((n+1)*(n+2)).*(1-cos(2.*T/s));
B(n+1)=(f(n+1)./4)*s*sqrt((n+2)*(n+3)).*(1-cos(2.*T/s));
C(n)=(f(n)./2)*s*(n+1)*sin(2.*T/s);
C(n+1)=(f(n+1)./2)*s*(n+2)*sin(2.*T/s);
W(n)=0.5*(n.*(f(n)^2)+(n+1).*(C(n)^2)+(2*n+2).*(B(n)^2)+(n+1).*f(n).*f(n+1)+C(n).*C(n+1).*(n+2)+2*sqrt((n+1)*(n+3)).*B(n).*B(n+1));
w=w+W(n)
end
T=0:0.5:20;
plot(T,w)
这个程序的主要目的就时要做出w随T变化的图,w是一个含有T的求和公式,n从0到1000求和 ,0《=T《=20,步长0.5
1.
你的B0,B1,C0,C1都是行矩阵,即在你的程序中都是1行41列的矩阵,而只有矩阵是方阵时才能平方,因此
w=0.5*(0*(f0^2)+(0+1).*(C0^2)+(2*0+2).*(B0^2)+(0+1).*f0.*f1+C0.*C1.*(0+2)+2*sqrt((0+1)*(0+3)).*B0.*B1); 语句中的C0^2,B0^2是错误的,因为C0,B0都不是方阵,后面的有类似写法的也出错.
如果你要行矩阵乘以自身,那应该乘以自身的转置矩阵,写法是:a*a'
所以C0^2改为C0*CO',B0*B0',依此类推
2.B(n)=(f(n)./4)*s*sqrt((n+1)*(n+2)).*(1-cos(2.*T/s)); 有错
当n=1时,上式等号右端的
(f(n)./4)*s*sqrt((n+1)*(n+2)).*(1-cos(2.*T/s))因为T,所以有41个值,而
B(1)只能储存一个值,因此是不正确的,后面的语句有类似用法的也错了.
由于我不知道你要表达什么,所以不知道怎么改
你的B0,B1,C0,C1都是行矩阵,即在你的程序中都是1行41列的矩阵,而只有矩阵是方阵时才能平方,因此
w=0.5*(0*(f0^2)+(0+1).*(C0^2)+(2*0+2).*(B0^2)+(0+1).*f0.*f1+C0.*C1.*(0+2)+2*sqrt((0+1)*(0+3)).*B0.*B1); 语句中的C0^2,B0^2是错误的,因为C0,B0都不是方阵,后面的有类似写法的也出错.
如果你要行矩阵乘以自身,那应该乘以自身的转置矩阵,写法是:a*a'
所以C0^2改为C0*CO',B0*B0',依此类推
2.B(n)=(f(n)./4)*s*sqrt((n+1)*(n+2)).*(1-cos(2.*T/s)); 有错
当n=1时,上式等号右端的
(f(n)./4)*s*sqrt((n+1)*(n+2)).*(1-cos(2.*T/s))因为T,所以有41个值,而
B(1)只能储存一个值,因此是不正确的,后面的语句有类似用法的也错了.
由于我不知道你要表达什么,所以不知道怎么改
clearclcr=0.1;s=0.01;T=0:0.5:20;w=0;f0=(-tanh(r)^0./cosh(r))
请问如何在matlab中画函数图形 r=cosh(t)*exp(-t/100) 为什么会出现 Inner matrix
cos(x)*cosh(x)+1=0 solve('cos(t)*cosh(t)+1=0') ans = 1.87510
matlab采样问题t=0:0.01:2*pi;f0=5;fs=20; %采样频率y=cos(2*pi*f0*t); %
设R={,},试求r(R),s(R)和t(R)
设A={1,2,3},给定A上二元关系R={,,},求r(R),s(R)和t(R).
C语言T=R∩S T=R∪S T=RXS
char r[20]={'I','a','m',' ','s','t','u','d','e','n','t','\0'
【符号的意思】W=I^2*R*t
m=pv v=s/t p=m/vG=mg P=F/S W=Fs P=W/t I=U/R U=IR R=U/I W=UIt
有关微积分1)单说第三问的R上热量,0-4s在R上的热量的问题,R=1欧,Q=R*t*I^2,用三角形面积s=0.5*底
设R是集合A={a,b,c,d}上的二元关系,R={,,,}求r(R),s(R),t(R)