lzth.net
当前位置:首页 >> mAtlAB菜鸟求助大神帮忙变出程序求出p(r,t)和F(t)的表达式,急用啊。。。穷鬼,只有20分... >>

mAtlAB菜鸟求助大神帮忙变出程序求出p(r,t)和F(t)的表达式,急用啊。。。穷鬼,只有20分...

function [p]=prt(r,t)if(t>=0 & t<= r) p=1110*exp(-6.185*t)else p=f(t-r)*exp(6.185*t)end上面是第一个函数,存为prt.m文件下面是第二个函数,存为myft.m文件function [f]=myft(t,r)f=mybeta(t)*int(1.0/768*(r-15)^4*exp((r-15)/2)*p(r,t),15,49)

f1=sym('t*(heaviside(t)-heaviside(t-1))'); f2=sym('(heaviside(t-1)-heaviside(t-2))'); f3=f1+f2; figure; subplot(311),ezplot(f3,[0,10]); title('f(t)'); f4=subs(f3,t,t-4); subplot(312),ezplot(f4,[0,10]); title('f(t-4)'); f5 = subs(f3,t,4-2*t); subplot(313),ezplot(f5,[0,10]); set(gca,'Ylim',[0,1.2]); title('f(-2t+4)');

n=28;syms p;a=n/2:n;b=n/2:-1:0;sum=0;for i=1:length(a) sum=sum+nchoosek(n,a(i))*p^(a(i))*(1-p)^(b(i));endp=double(solve(sum-0.9,p));disp(p(p==real(p)));它一个有28个复数域的解,存在结果p中,我最后显示了了两个实数域的解,如下: -0.2332 0.6013如果只需要正数,将最后一行改成下面的:disp(p(p(p==real(p))>0));希望采纳,设置为满意答案,谢谢.

上述程序不变,加上:R=subs(R,t,x)S=subs(S,t,x)I=subs(I,t,x)就可求出当t=x时的对应值.

syms x m t d e r z>> d=solve('cos(x)=(d^2+e^2-r^2)/(2*d*e)','d')d = (e^2*cos(x)^2 - e^2 + r^2)^(1/2) + e*cos(x) e*cos(x) - (e^2*cos(x)^2 - e^2 + r^2)^(1/2)>> z=solve('cos(z)=(r^2+d^2-e^2)/(2*r*d)','z')z = acos((d^2 - e^2 + r^2)/(2*d*r)) -acos((d^2 - e^2 + r^2

N=input('输入一个正偶数:');%%%输入一个偶数a=primes(N)%%%%%%求不大于N的质数;for i=1:length(a) %%%%length(a)表示质数的个数 循环从1到length if max(primes(N-a(i)))==N-a(i) %%N-a(i)表示N与第i个质数的差,max(primes(N-a(i)))最大的这个质数与N-a(i)相等 disp(['The two primes is ',num2str(a(i)),' and ',num2str(N-a(i))]) %%显示结果 num2str(a(i))表示数字转化成字符串 break; %%%%结束循环 end end

用matlab求解微分方程组的解析解,可以用dsolve()函数,如求其数值解,则可以用ode45()函数.解析解的格式[V1,V2,r]=dsolve('微分方程组表达式')数值解的格式[t,x]=ode45(@ode_fun,tspan,X0)由于题主给出的已知值和初始条件书写比较撩草,容易混淆.请重新上传清楚的表达式和条件.

n=1024;s1=zeros(1,200);s2=ones(1,200);ss=[s1,s2];fs=fft(ss,n)/nsubplot(1,2,1)x=3000/2*linspace(0,1,n/2+1);y=2*abs(fs(1:n/2+1))plot(x,y);grid;e=0.1;f=1000;wn=2*pi*1000;n1=wn^2;n2=[1 2*wn*e wn^2];g=tf(n1,n2);a=0:199;h=step(g,a/3000);fh=fft(h,n)/n;y=2*abs(fh(1:n/2+1))subplot(1,2,2)plot(x,y);grid;

明确回答:不能反解出t的显式解析式复杂的超越方程,即便计算机符号计算能力再强,也是望方程兴叹.

clear;clc;f=@(t)(t+1)/2.*t.*log(2./t);fplot(f,[0,5]);x=fminsearch(@(x)-f(x),2);hold on;plot(x,f(x),'r*');legend('f(t)',sprintf('(%.3g , %.3g)',x,f(x)))

相关文档
网站首页 | 网站地图
All rights reserved Powered by www.lzth.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com